THE screvesn 
BELL SYSTEM 
TECHNICAL JOURNAL 


J. A. Arnaud Pulse Broadening in Multimode 
Optical Fibers 


J. Shah Effects of Environmental Nuclear 
Radiation on Optical Fibers 


C. M. Miller Loose Tube Splices for Optical 
Fibers 
G. Miller The Effect of Longitudinal 


Imbalance on Crosstalk 


T. J. Aprille, Jr. Wideband Amplifier Design Using. 
Major Multiloop Feedback Techniques 


A. Descloux Variance of Load Measurements in 
Markovian Service Systems 


D. R. Means Tl Carbon Transmitter Model for Use 
in Computer-Aided Analysis of Tele- 
phone Set Transmission Characteristics 


M. J. Gans and Some Far-Field Studies of an Offset 


R. A. Semplak Launcher 
A. D. Wyner Upper Bound on Error Probability 


for Detection With Unbounded 
Intersymbol Interference 


Contributors to This Issue 


1179 


1207 


1215 


1227 


1253 


1277 


1301 


1319 


1341 


1353 


THE BELL SYSTEM TECHNICAL JOURNAL 


ADVISORY BOARD 


D. E. PROCKNOW, President, 
Western Electric Company, Incorporated 


W. O. BAKER, President, 
Bell Telephone Laboratories, Incorporated 


W. L. LINDHOLM, Vice Chairman of the Board, 
American Telephone and Telegraph Company 


EDITORIAL COMMITTEE 


W. E. DANIELSON, Chairman 


F. T. ANDREWS, JR. J. M. NEMECEK 

S. J. BUCHSBAUM Cc. B. SHARP 

i. DORROS B. E. STRASSER 

D. GILLETTE D. G. THOMAS 
W. ULRICH 


EDITORIAL STAFF 


L. A. HOWARD, JR., Editor 
P. WHEELER, Associate Editor 


J. B. FRY, Art and Production Editor 
F. J. SCHWETJE, Circulation 


THE BELL SYSTEM TECHNICAL JOURNAL is published ten times a year by 
the American Telephone and Telegraph Company, J. D. deButts, Chairman 
and Chief Executive Officer, R. D. Lilley, President, C. L. Brown, Executive 
Vice President and Chief Financial Officer, F. A. Hutson, Jr., Secretary. Checks 
for subscriptions should be made payable to American Telephone and Telegraph 
Company and should be addressed to the Treasury Department, Room 1038, 195 
Broadway, New York, N. Y. 10007. Subscriptions $15.00 per year; single copies 
$1.75 each. Foreign postage $1.00 per year; 15 cents per copy. Printed in U.S.A. 


THE BELL SYSTEM 
TECHNICAL JOURNAL 


DEVOTED TO THE SCIENTIFIC AND ENGINEERING 
ASPECTS OF ELECTRICAL COMMUNICATION 


Volume 54 September 1975 Number 7 


Copyright © 1975, American Telephone and Telegraph Company. Printed in U.S.A. 


Pulse Broadening in Multimode Optical Fibers 


By J. A. ARNAUD 
(Manuscript received February 21, 1975) 


Closed-form expresstons are obtained for the impulse response of 
graded-index fibers whose relative permittivity is a homogeneous function 
of the two transverse coordinates x, y, and for the impulse width in graded- 
index fibers whose profile departs slightly, but otherwise arbitrarily, from 
a square law. The inhomogeneous dispersion of the material is taken into 
account. Pulse broadening can be reduced by a factor of 12 from the value 
obtained for square-law fibers. Simple expressions are found for the 
acceptance of highly oversized fibers. 


1. INTRODUCTION 


Light-emitting diodes supply their optical power in a time and space 
incoherent form. The line width is typically of the order of 200 A, 
and the radiation is approximately lambertian with an emissive area 
of the order of 50 X 50 wm. Time and space incoherent optical pulses 
can be transmitted by oversized optical fibers. However, optical pulses 
propagating in such fibers tend to broaden as they travel. This is in 
part due to the nonzero line width of the source and the dispersion 
(d?k/dw*) of the fiber material. The other cause of pulse broadening 
is associated with the fact that the time of flight of a pulse along a 
ray depends on the ray trajectory. Pulses traveling along axial rays 
usually go faster than pulses traveling along rays of large amplitude. 
Because both types of rays are excited by spatially incoherent sources, 
the difference in axial group velocity causes a broadening of the input 
pulse. In the main text of this paper, we assume that the carrier is 
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monochromatic and that the spatial distribution of the rays is time- 
invariant. This is the case, for instance, when the source is an injection 
laser that oscillates simultaneously on many transverse modes. The 
difference in frequency between these various transverse modes can 
usually be neglected. 

It was first pointed out by Kompfner! that pulse broadening in 
step-index fibers could be drastically reduced by introducing ray 
equalizers at various locations along the fiber. The role of ray equalizers 
is to exchange fast and slow rays. A possible implementation of this 
idea is shown in Fig. 1 together with the calculated impulse response 
for uncorrected and corrected step-index fibers.2, Because natural 
mode mixing appears to be very small in the most recently made 
optical fibers, ray converters may be practical. They have not been 
experimented with, however, and we shall therefore restrict ourselves 
to uniform, uncorrected fibers. 

Important results concerning the broadening of spatially incoherent 
optical pulses in graded-index fibers have already been reported. In 
Refs. 3 to 9, the difference in group velocity between the various 
modes that can propagate in step-index and graded-index fibers has 
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Fig. 1—Ray converter that minimizes pulse broadening in step-index fibers. (a) 
Angular spread of a step-index fiber. (b) Optical arrangement with confocal lenses. 
The first and last lenses are unconventional. (c) Calculated impulse response for 
uncorrected [P(é) = 1 within the pulse ] and corrected step-index fibers (from Ref. 2). 


1180 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1975 


been evaluated. The impulse response is obtained by adding the con- 
tribution of each mode, under the assumption that all modes are 
equally excited by the source. The calculation of the group velocities 
can be simplified with the help of the W.IX.B. approximation (see, 
in particular, Ref. 4). 

Let us now describe an alternative ray-optics method. The time of 
flight of a pulse along a ray is first evaluated according to the laws of 
geometrical optics. A ray can be defined by the point x, y where it 
intersects the input plane of the fiber (plane z = 0), and by the 
transverse components, k,, k, of the wave vector k. k is, by definition, 
directed along the ray and has magnitude (27/\o)n, where \» denotes 
the free-space wavelength and n the refractive index of the fiber 
material, usually a function of « and y. Thus, the time of flight of a 
pulse (at a fixed carrier frequency) is, in general, a function of the 
four parameters 2, y, kz, ky. These four parameters can be considered 
the components of a four-vector p, in the so-called phase space. The 
impulse response is subsequently obtained by assuming that the 
density of rays is equal to (27)-? in the phase space. In other words, we 
assume that the number of rays whose points of intersection with 
the input plane are between x, x + dx and y, y + dy, and whose 
direction is defined by values of k,, ky between kz, kz + dk, and ky, 
ky + dky, is equal to dxdy dk.dk,/(27)?. The total power transmitted 
is the acceptance (or number of modes) of the fiber. This is the power 
transmitted for a source of luminance unity (see, for example, Ref. 10). 

The approach used in Refs. 11 to 13 is based on the conventional 
ray equations. We have shown in Refs. 14 and 15 that it brings a 
considerable simplification to write the ray equations in the Hamil- 
tonian form. The relationship between the ray-optics method and the 
W.Ix.B. method becomes more obvious with the Hamiltonian form. 
It can be shown that the W.Ix.B. method and the ray-optics method 
are essentially identical.' 

An important difference, however, should be noted. In the W.Is.B. 
method, modes whose axial wave number k, is less than the free wave 
number k, in the surrounding medium (or cladding) are assumed to 
leak out so rapidly that they can be ignored. On that basis, the accep- 
tance of a step-index round fiber with radius a, for example, is found 
to be 


N? = (k? — k2)a2/2 = V2/2. 


The radiation loss of leaky modes can be small in the case of highly 
oversized fibers, however, as was pointed out by Snyder.!® The ray- 
optics condition is distinctly different: Only those rays are ignored 
whose tangential component of the wave vector at the core-cladding 


PULSE BROADENING 1181 


interface [ (k2 + k%)?, where k, denotes the azimuthal wave number ] 
is less than the free wave number k, in the surrounding medium. 
According to ray optics, the acceptance of a step-index fiber is N? = V? 
instead of V?/2. The influence of the slightly leaky rays on the im- 
pulse response of fibers has not been observed. This is perhaps because 
high-order modes are more sensitive to irregularities than low-order 
modes. Slightly leaky rays may become important when highly over- 
sized fibers of good quality are fabricated. This is even more so for 
graded-index (e.g., near-square-law) fibers, because the field decays 
exponentially beyond the caustic line, which bounds the ray 
trajectories. 

In most previous works, the effect of inhomogeneous dispersion* on 
quasi-monochromatic pulse broadening was neglected. This effect, 
however, was taken into account for square-law and linear-law graded- 
index fibers in Appendix B of Ref. 14, and by Gambling and Matsuhara® 
for circularly symmetric modes in square-law fibers perturbed by an 
r term. The result for arbitrary small deviations from square-law was 
given by Arnaud in Ref. 15. Olshansky and Keck? first pointed out 
that inhomogeneous dispersion is of great practical importance, at 
least for fibers doped with TiO:. Dispersion for the promising GeO» 
doped fibers is not known at the time of this writing. The variation of 
the loss of that material as a function of doping is likewise unknown. 
If we consider further that the sources used in pulse broadening experi- 
ments are not fully characterized in terms of their distribution in 
phase space, it appears that a precise comparison between theory and 
experiment is difficult at the moment. We shall therefore restrict our- 
selves to the theoretical evaluation of pulse broadening. 


Il. GENERAL RESULTS 


The derivations of the general results given in this section appear in 
Appendix A. They follow in a straightforward manner from the Hamil- 
ton equations for pulse trajectories in space-time. 

Fibers are most often characterized by a refractive-index profile: 
n(x, y, w). However, the quantity that enters directly in the expressions 
for pulse broadening is the square of the wave number k’(z, y, w) 
= (27/)o)?n? (a, y, w), Where \» denotes the wavelength in free space. 
We shall therefore deal directly with k?(a, y, w). 

Let x(z), y(z) denote a ray trajectory. Assuming that the fiber is 
time-invariant and uniform and that the material is isotropic, we ob- 


* Inhomogeneous dispersion refers to the spatial variations of the ratio of the local 
phase to group velocities in the material. This parameter should not be confused with 
the parameter d?k/dw?, usually refered to as ‘‘material dispersion.” The latter is 
important only for broadband carriers. 
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tain from the ray equations the following differential equation (see 
Appendix A): 


Lee (X + Y)/det = hk? — k2 + Xok?/aX + Yak?/ay, (1) 


where we have set, for convenience, X = x?(z), Y = y?(z). The quantity 
k. in (1) denotes the axial (z) component of the wave vector k and 
is a constant of motion. In other words, kz remains the same along 
any given ray. In a wave theory, k, corresponds to the propagation 
constant of a mode (sometimes denoted 8). Note that, in spite of the 
fact that we are using the language of wave optics, the theory given 
in this paper is based strictly on ray optics, except when we impose 
the condition k, > k, to make contact with previous results. 

It follows from the space-time ray equations that the ratio of the 
time of flight of a pulse along a ray to the corresponding time on axis 
is (see Appendix A) 


t = (ko/kz) (0k?/dw*)/ (dki/duw*), (2a) 


where ky = k(0O, 0, w). The sign ( ) denotes an average over a ray 
period. For any function a(z, y, w), we have defined 


(a(x, , 0) = 2" f alee), yO), ole (2b) 


where 2x(z), y(z) denotes a particular ray trajectory and Z the ray 
period. If the ray trajectory is not periodic, (a) should be understood 
as the limit of the right-hand side of (2b) when Z — o. In the special 
case where the inhomogeneous dispersion of the material can be 
neglected, & is proportional to w and, consequently, dk?/dw? = k?/w?. 
Thus, (2a) reduces to 


t = (k*)/ (kok). (2c) 


Finally, if the source of rays has a distribution f(p) in the phase 
space p = {2, y, kz, ky}, the response of the fiber to an input P’(#) is 
(see Appendix A) 


PW) = f PT t@@)7@) @). (3) 


The quantity 7(p) is the transmission of a ray (usually 7 < 1), and 
(dp) = dxdy dk,dk,. In the special case of a uniform lambertian source 
of luminance unity, we have f(p) = 1/(2z7)*. For simplicity, we can 
assume that 7'(p) is unity when the point zx, y is within the core cross 
section and the components k;, ky of p are within some area to be 
defined later for specific examples and zero outside that area. All the 
subsequent results follow from (1), (2), and (8). 
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Ill. IMPULSE RESPONSE WHEN k°(x, y) — k3 IS A HOMOGENEOUS 
FUNCTION OF x AND y 
Let the differential equation (1) be averaged over a ray period. 
The left-hand side of (1) vanishes because d(X + Y)/dz assumes the 
same values at the ends of the integration interval. (In this integra- 
tion, k, can be considered a constant.) Thus, we have 


(k2 — k2 + Xoak?/aX + Yak?/aY) = 0. (4) 


Let us further assume that h(X, Y) = k®(X, Y) — kG is a homogeneous 
function of degree x in X = 2? and Y = 7’. This means that, for any X, 


hQX, AY) = Mh(X, Y). (5) 
If we differentiate (5) with respect to \ and set \ = 1, we obtain 
XOh/AX + Yah/aY = xh(X, Y). (6) 
Thus, going back to k?(z, y, w), 
XOk?/0X + Yok?/aY = x(k? — ko). (7) 


In that case, a simple and general expression for the relative delay 
in the absence of material dispersion is readily obtained from (2c), 
(4), and (7), 


t = L(k./ko) + «(ko/kz) ]/( 4+ x). (8) 


The relative delay ¢ is plotted in Fig. 2 as a function of k./ko with x 
as a parameter. This result is applicable, for example, to the index 
profile 


(x,y) = ko —alx| — Blyl, (9) 


where a and 8 denote constants. In that example, x = 3. Note that 
the fiber described by (9) is not circularly symmetric, even if a = 8. 
Examples of circularly symmetric fibers that satisfy (5) will be given 
in the next section. 

In almost any z-invariant focusing system, any initial distribution 
eventually reaches a steady state. This steady state in general differs 
from the initial distribution. A lambertian distribution f = constant, 
however, remains lambertian because it is a (trivial) solution of the 
Liouville equation (see Appendix A). Note that the distribution f in 
(3) represents a ray (or pulse) density. If the medium introduces a 
nonuniform attenuation on the rays, the power density 7'(p)f(p) in 
phase space needs to be distinguished from the distribution f. 

A fiber is usually surrounded by a homogeneous material, called 
the cladding, with wave number &,. For fibers that are not highly 
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Fig. 2—Relative time of flight in a fiber where k?(z, y) — k? is a homogeneous 
function of degree x in z and y [ky = &(0, 0)]. For most fibers, k,/k» is close to unity. 


overmoded, the transmission law 


T (x, y, kz, ky) = : (10) 
is often acceptable. Equation (10) says that rays whose axial wave 
number is less than the free wave number in the surrounding medium 
are leaking sufficiently rapidly to be ignored. The distribution f of 
the lambertian source is set equal to 1/(27)? so that the luminance is 
unity. In that case, the total power transmitted is the acceptance of 
the fiber. The relative time of flight is, within the present assumptions, 
solely a function of k,. The upper and lower bounds on k, are k(a, y) 
and k,, respectively. It remains to express the volume element 
dk,dk,dxdy in (3) as a function of dk, dx, dy. For given x, y, a constant 
value of k, corresponds to a circle of radius squared k?(z, y) — k3 in 
the kz, ky space because kz + k}? = k? (ax, y) — k2. Thus, 


dkdky = mdk?. (11) 


Let us evaluate the acceptance of the fiber. The light acceptance of 
any optical system is, as we have seen, the volume in phase space of 
the accepted rays divided by (27). It is also equal to the effective 
number of modes that the system can transmit. If we integrate P(t) 
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fromt = —o© tol = +. in (3), the integral over P’ in the integrand 
is unity, and we obtain 


N* = (1/4x) ff CeCe, v) — bdzay, (12) 


where we have used (10) and (11). Thus, 47N? is the volume enclosed 
by the profile: k?(x, y). For a step-index fiber of any shape with cross- 
section area A, for example, we have from (12) 


N? = (A/4r) (k? — k2). (13) 
This expression should be multiplied by 2 to take into account the 


two states of polarization. 
The pulse transformation in (3) becomes, using (11), 


P(t) = (\/4n) / angi is ee Dice — 4(b.) dk. (14) 


Let the input pulse P’ (¢) be a symbolic 6 function (e.g., a rectangular 
pulse of width At and height At in the limit At 0). The output 
pulse in (14) becomes 


P(t) = (1/4m) |dkz/dt|A(k.), ka > Is, (15) 


where |dk2/dt| denotes the absolute value of dk?/dt and A (k.) denotes 
the cross-section area that satisfies k(x, y) > kz. kz can be expressed 
as a function of the delay ¢ by inverting the relation ¢(k,) given earlier. 
We obtain, from (8), 


dk?/dt = 2(1 + «)kL1 — (x/k,?)], (16) 
where 


k, = k./ko = (1+ «)t(/2& {[0 + #t/2P — x}. (17) 


If x > 1, there is only one value of k; between k, = k;/ko and 1, for 
any k,. If 
koe, (18) 


there are two values of k; that need be considered. Their contributions 
to P should be added. If « < k?, there is again only one relevant value 
of kz. 
Let us consider as an example a (noncircularly symmetric) square- 
law medium 
kK? (a, y) = ke (1 — 0? — ay’), (19) 


where 2,, Q, denote arbitrary constants. 27/Q, and 27/Q,, for small 
x, y, are the periods of ray oscillation in the xz and yz planes, respec- 
tively. The area A (k.) defined earlier is the interior of an ellipse 


A(kz) = r(1 — k2)/0,0,. (20) 
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The impulse response is obtained from (15) and (16) with « = 1, 
and (20), 
kle. | (995), keh 


PO = ae os 


(21) 
where, from (17), k; = t — (# — 1)}. Because, in most fibers, k, re- 
mains close to ko, the variation of k; can be neglected, and the pulse 
response is almost rectangular. 

For a step-index fiber, the area A is the area of the core cross sec- 
tion, and t = k)/k,. Thus, the impulse response of a step-index fiber 
with cross-section area A is simply 


PQ) = kjA/2r 6, 1<t< ko/k,. (22) 


Because, in most fibers, t remains close to unity, the pulse response is 
almost rectangular.? The pulse width, however, is considerably larger 
than for square-law fibers, as we shall see in more detail later. 


IV. CIRCULARLY SYMMETRIC FIBERS WITH k’— k? A POWER 
OF THE RADIUS 


Let the wave-number profile be of the form 
| (BR, w) = ko(w) — kilo) Rs, (23) 


where R = X + Y = 7? denotes the square of the radius. The relative 
time of flight is, substituting (23) in (2a), 


t = (ko/kz) (Ok2/de2)/ (dk3/dw?) 


= (ko/Kz)(1 — &Dx(*)), (24) 

where we have defined 
& = ki/ko (25) 
Dg = ko (dhit/dw*) /kee (dko/ de’). (26) 


D, is a dispersion factor equal to unity in the absence of dispersion. 
Thus, we need to evaluate (R*). It is interesting that we can do that 
without solving the ray equations. The quantity (R*) is, of course, 
independent of dispersion, so we may omit the w arguments. 

lor circularly symmetric fibers, (1) can be written 


LPR /de = d(k2R)/dR — k?. (27) 
Averaging (27) over a ray period, we obtain 
k2 = (d(k’R)/dR). (28) 
We have also, directly from (23), 
(k?) = ko — ki(R*) (29a) 
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and, from (28) and (28), 


(kh?) = (kz + xko)/(1 + «). (29b) 
quating the two expressions (29a) and (29b) for (k?), we obtain 
(Rt) = (1 — k?)/(1 + x), (30) 


where k; = k./ko. Thus, substituting (R*) from (30) into (24), the 
relative time of flight is 


t = k,~* — Dy(k,-* — k)/(1 + «). (31) 


In applications, we need kz as a function of t. Solving (31) for k; and 
setting D, = D,/(1 + «), we obtain 


k, = (t/2D,) & CG/2D,)? + 1 — Di}. (32) 
By differentiating (32), we further obtain 
dk,?/dt = 2k,[.D, — (1 — Dp)/k?}". (33) 


To obtain explicitly the impulse response in (15), we need the area 
A(k,) defined by k, < k(R). For k(R) in (28), this area is 


A(kz) = rR(kz) = r[(1 — k?)/e |". (34) 


If «, were kept a constant as the parameter « varies, the core radius a, 
defined by k(a) = k,, would vary. Thus, it is preferable to express « 
as a function of the core radius a. We have 


le = (1 = beat, (35) 


where k, = k,/ko. The impulse response is finally obtained from (15), 
(33), (84), and (35) ; 


P(t) = (koa?/2)k.L (1 — he?)/(1 — ke?) "/LDe — (1 — DOke~*]. (36) 


The possibly doubled value k; is expressed as a function of t by (82). 
Thus, a closed-form expression has been obtained for the impulse 
response of a fiber with k? — k§ a power of r, that takes inhomogeneous 
dispersion into account. 

In the absence of dispersion, we have D, = 1/(1 + «), and (86) 
reduces to 


P(t) = (kga?/2)k.0 (1 — k”)/(L — hs) "(1 + «)/(. — «hz). (87) 


As indicated in the previous section, there are in general two values 
of k, that contribute to P. Note that the shape of the impulse response 
does not depend on the core radius a. 

The impulse response P(t) in (87) is shown in Figs. 3 and 4 for 
various values of the parameter x. These curves are essentially the 
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0.990 0.992 0.994 0.996 0.998 1.000 1.002 1.004 1.006 1.008 1.010 
t 


Fig. 3—Impulse response of a circularly symmetric fiber with k?(r) = kj —_kir*« 
for a lambertian source and various values of x. The optimum impulse response is for 
ek = 0.9. 
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k', =ks/kg = 0.9 


co (STEP—INDEX FIBER) 





~ (2—k’s) 
Fig. 4—Continuation of Fig. 3 for larger values of x. 
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same as those shown in Ref. 4. Figures 3 and 4, however, are much 
more detailed. We have assumed that k,/ky) =0.9, that is, An/n=10%, 
a rather large value. For «x = 1 (square-law fiber), the pulse width + 
is 0.0054. For example, if n = 1.45 and the fiber length is 1 km, the 
pulse width is 26 ns. For x = 0.9, however, the corresponding pulse 
width is only 7 ns. We find, in agreement with Ref. 4, that the mini- 
mum pulse width occurs when « = k;. For a step-index fiber (x >), 
the pulse width would be as large as 630 ns. Note the following detailed 
features on the curves in Figs. 3 and 4. For (0.9)? < « < 1, the re- 
sponse starts from infinity because of the minimum in the ¢(k,) curve. 
For « = 0.85, P drops suddenly for t } 0.998. This is because, at that 
time, the smaller of the two k,; becomes less than 0.9, and is rejected. 
For x = 0.95, the response crosses the ¢ = 1 axis. 

Tigure 4 is applicable to larger values of x. We note that, for a 
very large « (step-index fiber), the response is almost rectangular. The 
slow decay in power shown in I’ig. 4 would be almost negligible for 
small An/n. 

The effect of inhomogeneous dispersion is shown in Fig. 5. The 
parameter x is kept equal to 0.9 (this is the optimum value in the 
absence of inhomogeneous dispersion), but D, is made to vary in the 


30 


20 


P IN DECIBELS 





0 
0.990 0.995 1.000 1.005 1.010 
t 
Fig. 5—Impulse response for a fiber with k?(r) = k2 + k2(r?)°° for various values 
of the parameter D that expresses the inhomogeneous dispersion of the material. 


D = 1 corresponds to the absence of dispersion. D ~ 1 merely introduces a shift 
in the optimum value of x. 
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neighborhood of unity. These curves have a striking resemblance to 
those in Fig. 3. This means that the effect of inhomogeneous dispersion 
merely consists in shifting the optimum value of «. The impulse re- 
sponse remains essentially the same, at least for x 1. 

The total pulse power is the acceptance of the fiber, a function of x. 
The acceptance is, in the present case, 


Nee [> Poa 24 i; *Tee(R) — k2 Jak 


= (4) (ko — ks)a? — kia? (a?)*/(x + 1) ] 
= [«/4(« + 1) ](ko — ks)a’. (38) 


The coefficient in the last expression in (38) is + for step-index fibers 
(x +0) and # for square-law fibers. The acceptance given in (38) 
should be multiplied by 2 to account for the two states of polarization. 
The same rule applies to all the expressions given in this paper. It is 
more difficult to obtain the ray-optics acceptance of fibers. The result 
is derived in Appendix B. 

In the next section, we consider fibers whose profile departs slightly, 
but otherwise arbitrarily, from a square law. 


V. NEAR-SQUARE-LAW FIBERS 


Let us rewrite the differential equation (1) for circularly symmetric 
fibers 


4h2PR/dz = d(k?R)/dRk — k?. (39) 
For square-law fibers 
k2(R) = ki — kik, (40) 
the solution of (39) is 
R(z) = Ro + (Ro — 12/kj)? cos (202/k:), (41) 
where 
Ro = (ko — k2)/ki = 301 — ke?) /@ (42) 


and Q = ky/ko. We have introduced in (41) the axial component of the 
angular momentum (or Bouguer invariant) 


be = thy — oka; (43) 
which is the second constant of motion. Let us set 
6 = (1./kiR a)’, (44) 


where Ry denotes the maximum radius squared. Note that, for merid- 
ional rays, 6 = 0 and, for helical rays, 6 = 1. Equation (41) can be 
written in the convenient form 


R = sRy(1 + 6) + $Rar(1 — 8) cos (292/k:). (45) 
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Tor later use let us evaluate (R”), the average of R” over a ray 
period. Using the binomial expansion and the result 


{cos”) = m!2-"[ (m/2)!}? (46) 


for m even and 0 for m odd, we obtain 





ieee, Sain n (1 — G)?-7 =e 0)7 
oe Ric 2... 2(n — m)'[ (m/2) PP (47) 
In particular, 
(R?) = Rir(8@ + 26 + 3)/8 (48a) 
(R) = R3,(1 + 0) (502 — 26 + 5)/16 (48b) 


(R*) = Rj, (3564 + 2063 + 186? + 206 + 35)/128. (48c) 


Let us now show that a closed-form expression can be obtained for 
the times of flight in fibers whose permittivity profiles depart slightly 
from a square law. Inhomogeneous dispersion is taken into account. 
Let the profile be of the form 


N 
ie(R) = kj — MR + Do BAR”. (49) 


We assume that ¢,f", n = 2, is of the order e << 1, where €, = k2/kj. 
Substituting (49) in (2a), we obtain (with Q? = e«; = ki/k@) 


t= Wi (1 — D,OX(R) + > é,D,(R*)), (50) 


where we have defined inhomogeneous dispersion factors 
D, = (kedki/ dw?) / (kadk§/dw*). (51) 


The D,, are unity in the absence of inhomogeneous dispersion. Because 
the perturbation is small, (R”) in the sum (50) can be replaced by 
the expression (47) applicable to square-law fibers. This approximation 
is not permissible, however, for the term (Ff) in (50) because this term 
is not small. We need an exact expression for (R). We proceed as in 
the previous section. We first observe that, for k? in (49), 


N 
d(k?R)/dR = 2k? ~ kk + XS (n — 1)kZR*. (52) 
n=2 
Integrating (89) over a ray period, the left-hand side vanishes and, 


using (52), we obtain an expression for (k?) that does not involve 


(R) 
(l8) = FQ3+ WB) + 4D — myaRer), (53) 
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We also have, directly from (44), 
N 
(kh?) = ko — ki(R) + Do kn(R™). (54) 
n=2 
Thus, by comparing (53) and (54), 


ER) = 408 — WD) +4 (+ DAR. (55) 


Substituting this expression for (R) in (50), we obtain the relative 
time of flight for circularly symmetric near-square-law fibers 


t= RML 40 — BPD. + Y [Da $+ WDilelR)}. (66) 


Alternatively, ¢ can be expressed in terms of the azimuthal and radial 
mode numbers. The result is given in Appendix C. 
In the absence of inhomogeneous dispersion, (56) reduces to 


t= POP +14 E 1 — we(R)] (57) 


Limiting ourselves to an r* correction to the square-law profile, 
és; = eg =---=0, and setting e. = e, (57) becomes, using (48), 


é=4[1 — pw(1 + 6) S32 — pul + 6) — ep (80? + 26 + 3)/8] 
1+ pul — 86) + 4 — 266+ (2 — 36)0)/16 + O(pix) 
pu = (Ry. (58) 


The first two terms in (58) give sufficient accuracy when px < 0.01, 
that is, when the total relative change in refractive index An/n  pa/2 
is less than 0.005 (p. = ?a?). 

The total pulse width 7 is the maximum variation of ¢for0 < @< 1 
and 0 < p < pa. For the square-law fiber [e = 0 in (58) ], we obtain 


r = p2/2 (ray optics). (59) 


It should be noted that, in defining 7 in (59), we have specified that 
the maximum radius of the ray be less than a for any 6. This condition 
is different from the condition used earlier that k, be larger than k,. 
The ray-optics condition py < pa is applicable to highly oversized 
fibers. 

If we now consider the expression in (58) with a correction term in 
m, we find that ¢ = 1 for meridional rays (@ = 0) when e. = ? in 
agreement with Ref. 17, where it is shown that all the rays have 
exactly the same optical length when k?(*) = [cosh (2) |? 21 — 2 
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+ (3)a4 + ---. We also find that ¢ = 1 for helical rays (@ = 1) when 
€. = 1, in agreement with Ref. 3, where it is shown that helical rays 
have exactly the same optical length when k?(r) = (1+ r)73 1 
—r+7r+-+ ---, By considering all rays whose maximum radius is 
less than a, we find that the minimum 7 is obtained for eg = 0.91. In 
that case, 7 = 0.04692. The improvement compared with square-law 
media is therefore as large as 11. If we had imposed instead the wave- 
optics condition k, > k,, the vertical scale in Fig. 6 would be divided 
by (1 + 6)?. For eg = 0, for example, the wave-optics pulse width is 
p2/8 instead of p2/2 as in (59). With the wave-optics limit, the opti- 
mum value of ¢, turns out to be 2 instead of 0.91. The improvement 
over the square-law case is only 4, instead of 11. 

Let us now consider the effect of 7° terms. Figure 7 shows the varia- 
tion of the pulse width 7, defined as the maximum variation of ¢ for 
any 0 < 6< landanyO < par < 0.002, as a function of 2 for various 
values of e3. The effect of ¢3 is essentially to shift the optimum value 
of «2 to lower values. The reduction in pulse width is rather modest. 
Nevertheless, a small improvement is obtained, compared to the case 


< 
Ry & 0.01 


2 
M 


(t—1)/R 





0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
0= (£,/Rm) 
Fig. 6—Variation of the normalized time of flight for a fiber with k?(r) = kj — Kir? 
tea (8/13) i in the absence of material dispersion for various values of the parameter 


€2. 6 = 0 corresponds to meridional rays and @ = 1 to helical rays. e: has been 
redefined to be dimensionless. 
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Fig. 7—Variation of the width of the impulse response with e2 for various values 
of e3 for a fiber with k?(r) = kh? — kir? + e2(ki/k2)r* + €3(k8/ki)r8. 


where e; = 0 when 


ko2k?(r) = 1 — p +. 0.615p2 + 709, pa = 0.002 


p = Or’, 600) 


We give only the result when the departures from a square-law 
profile are not circularly symmetric. The free wave number in the fiber 
is now in the form 


(x, yo) = BG) — MOR + LY kw)X'V!, G1) 


where, as before, X = v?, Y=y,R=2+7? = 1°. The ratio, t, of 
the time of flight along a ray to the corresponding time on-axis is found 
to be 


N on 
bom Teg Slee Re a a Xd [Dur — a(n + D1] 
x Eni (X! Yartyt (62) 


where €n, = k?,/kg and D,, is defined as D, in (51) with k, replaced 
by kx. Let us assume that it is permissible to use the sinusoidal rays 
of the square-law medium to evaluate the quantity (X'Y"~'). Because 
the average over one cycle of the product of powers of sinusoidal func- 
tions is known, the relative delay ¢ can be written in closed form. 


PULSE BROADENING 1195 


Let the ray trajectory be written 
x(z) = 2 cos (az + oz) (63a) 
y(z) = yo cos (az + $y). (63b) 


The coefficient a does not enter in the final result and is henceforth 
omitted. We evaluate 


(Xt sh) = (arya) 
= 1348 (Leos (¢ + $2) F*Leos (2 + 41) FO). (64) 


It can be shown that}8 


(Leos (2 + $2) Leos (2 + oy) P-?) 
nl — 
-20{aF ( l) )(,2 ,) eos [2s (dz a dy) | 


sol \2 —l— 8 
+O PVE oo 


(5) = a aye (66) 


Thus, given a ray trajectory, defined by the parameters Xo, Yo, oz, 
and ¢, (or, equivalently, by the values of 2, y, kz, and ky at the input 
of the fiber), we can evaluate in closed form the quantity (X'Y"~') 
that enters in formula (62) for the relative time of flight, from (64) to 
(66). 

The above calculation is incomplete for the following reasons. When 
the power law n?(r) of a fiber departs from the exact square law, pro- 
jected ray trajectories in the (ry) transverse plane are precessing 
ellipses.* That is, the principal axes of the near-elliptical trajectories 
slowly rotate as a function of z. This precession is unimportant for 
circularly symmetric fibers. For noncircularly symmetric fibers, how- 
ever, the ellipse precession introduces an averaging effect. I'urther- 
more, the noncircularly symmetric components of n?(r, ¢) change the 
eccentricity of the precessing ellipse. The axial component kz of the 
wave vector remains a constant, but the axial component /, of the 
angular momentum varies. Finally, in real fibers, slow (adiabatic) 
changes of the refractive index law along the fiber axis are likely to 
occur that must be taken into account. The twists of the fiber axis 
must also be taken into account. Thus, a realistic assessment of the 
effect of small noncircularly symmetric departures of the index law 


where 


*Tt is well known in mechanics that the only r* potentials (potential U ~ n?) 
that give closed trajectories are the harmonic potential U(r) ~ n?(r) = 1 — 7? and 
the Newton potential U(r) ~ n?(r) = 1/r. 
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from square law on pulse broadening requires a deeper and more 
intricate analysis than the one given in the present section. However, 
the result in (62) and (65) can be used as a basis for more complete 
analyses. 


VI. CONCLUSION 


From a rather straightforward application of the Hamilton ray 
equations, we have obtained closed-form expressions for the pulse 
width in graded-index fibers when k?(2, y) — ki is a homogeneous 
function of x and y, and for fibers whose profile departs slightly, but 
otherwise arbitrarily, from a square law. Inhomogeneous dispersion 
was taken into account. The expressions obtained are exact. The small 
angle (or weakly guiding) approximation need not be made. We have 
also given simple expressions for the wave optics and ray optics 
acceptance of weakly guiding graded-index fibers. 

The algebraic results given should prove more accurate and require 
much less computer time than the straightforward numerical integra- 
tion of time along ray trajectories. We have carried the perturbation 
only to first order in the small parameter e«. To obtain more accurate 
results, up to order &, we need a more accurate expression of the ray 
trajectory. This expression can be obtained, for example, by the 
method of strained coordinates.!® These more accurate expressions are 
probably not needed, however, in most practical cases. 


Vil. ACKNOWLEDGMENT 


The author expresses his thanks to EH. A. J. Marcatili for useful 
comments. 


APPENDIX A 
The Hamilton Equations of Ray Optics 


The Hamiltonian form of the ray equations are well known in 
mechanics and wave dynamics,” and they have also been used fre- 
quently in optics (e.g., Refs. 14, 15, 17, 21, and 22). However, their 
simplicity and power is not always appreciated. The physical difficulty 
is that it is not always recognized that ray momenta and wave vectors 
(or photon momenta) are identical concepts. On the other hand, ray 
momenta (proportional to the wave vectors) need be carefully dis- 
tinguished from mass-carrying momenta (proportional to the group 
velocities). On the mathematical side, we need distinguish a function 
such as k,(x, y) and the value kz assumed by that function. We must 
also be aware that da/dz denotes a total derivative, that is, in the 
present context, the variation of the quantity a along some given ray. 
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If a is a known function of x and y, and x = x(z), y = y(z) denote a 
known ray trajectory, then da/dz = (da/dx)(dx/dz) + (da/dy) (dy/dz) 
can be evaluated explicitly as a function of z. Here again, an arbitrary 
point in space x, y should not be confused with a specific ray trajec- 
tory « = z(z), y = y(z). Unfortunately, it is not possible to go into 
more details here. An excellent reference is Lighthill’s paper.” A com- 
parison between the W.K.B. method and the Hamilton equations is 
given in Ref. 14. 

Let X = (2, y, 2, ict) denote a point in space-time (¢ is time) and 
K = (kz, ky, kz, w/c) denote the four-wave vector, with w the angular 
frequency. An arbitrary medium is characterized by a function of K 
and X that we denote 


H(K, X) = 0. (67) 

The Hamilton equations for light pulses X(c), K(c) are 
dX/do = aH/aK (68a) 
dK/do = —0H/0X, (68b) 


where o denotes an arbitrary parameter. 

Equations (68a) and (68b) can be considered the basic postulates 
of geometrical optics. From a wave-optics point of view, (68a) follows 
from the requirement that the wave lengths and periods of the waves 
that constitute a wave packet be the same in the direction of a ray. 
Equation (68b) follows from (67), (68a), and the fact that K is the 
gradient of an eikonal function. Thus, in wave optics, the Hamilton 
equations (68) are derived from first principles and need not be 
postulated. 

Let — denote a point in phase space (kz, ky, w, x, y, t) at the input 
plane, and &’ a point in phase space at the output plane. The optical 
system maps the input phase space into the output phase space, that is, 


Sele). (69) 


It follows from (67) and (68) that the Jacobian of the transformation 
(69) is unity, a result often used in photometry. Equivalently, we can 
say that the determinant of paraxial ray matrices is unity or that the 
ray density in phase space is a constant of motion (Liouville theorem). 
These three statements are obviously equivalent, provided the rays 
are not reflected. 

A source of light that is time and space incoherent is described by 
a distribution S(£) in phase space. Each small volume in phase space 
can be pictured as an optical pulse, provided the dimensions of the 
volume are larger than unity. More precisely, this picture requires 
that AwAt>>1, Ak,Ax >> 1, and Ak,Ay > 1. The detailed structure 
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of the pulse is ignored in ray optics. Only the motion of the center is 
considered. 

The transmission 7, of an optical pulse through the optical system 
is a presumably known function of ~— that we denote as 


T, = T1(&). (70) 


For lossy optical systems, 7, < 1. Because the Jacobian of the trans- 
formation — — ~’ is unity, the output distribution is simply 


S'(é') = S(T). (71) 


The power emitted by the source and the power that can be collected 
at the output of the optical system are obtained by integrating S 
(or S’) over all variables, except ¢ (or t’). Thus, 


Pm = [ S@@® (72a) 


P(t’) = / S’(#) (de’), (72b) 


where € = (kz, ky, w, x, y) and @ is similarly defined. The terms (d@) 
and (dt’) denote elementary volumes in ¢ and @’ spaces, respectively. 
The response of the detector could be described by a function D(é’). 
For simplicity, we do not take the detector response into considera- 
tion. All subsequent results follow in a rather straightforward manner 
from the above results, through a succession of approximations. 

Let us assume that the properties of the fiber do not vary with time. 
This means that the Hamiltonian in (67), the transmission 7, the 
mapping ~— &’, and the pulse delay do not depend on time. In 
particular, 


t’ =t+ h(¢). (73) 
Sources that are t-separable, on the other hand, have the property that 
S(&) = P@F(%). (74) 


That is, the distribution in ¢-space does not vary with time. For a hot 
tungsten wire whose temperature varies with time, the spatial phase- 
space distribution is almost lambertian at all times, but the frequency 
spectrum (approximately given by the Plank law of radiation) varies 
with time. Thus, (74) is not applicable to that source. For consistency 
with (72), we assume that F'() is normalized to unity. 

For most sources, we can further assume that 


F(@) = Q(@)f(p), (75) 
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where p = (kz, ky, 2, y) denotes a point in spatial phase-space. That 
is, we assume that the spatial distribution does not depend on what 
part of the frequency spectrum we are considering. Both 0 and f are 
assumed normalized to unity. This ensures that F is normalized to 
unity. When the spectral width of the source is small (e.g., less than 
1 percent, as is the case for light-emitting diodes) and the fiber material 
absorption does not exhibit sharp resonances in that band, we can 
assume that 

T1(%) = To()T() (76) 
and 


t1(¢) = tow) + t(p). (77) 


For definiteness, we assume that the maximum value of J(w) is 

unity, and we define to(w) as the delay experienced by axial pulses. 

We evaluate in the main text ¢(p)/to at a fixed angular frequency. 
The pulse response is obtained from (71) to (77), 


Prt) = f PL — be) — Up)]2@) Tow) f(0)T ©) (dp) de 
= f PPL ~ b() 12 ©) Tole) de, (78) 
where 
pr) = f PL’ — t@)¥@)T@) @). (79) 


In writing (78) we have used the fact that the Jacobian of the trans- 
formation — — &’ is unity and that dw = dw’. The pulse response is 
the convolution of the pulse response in (79), which we may call the 
quasi-monochromatic pulse response, and the spectral width of the 
source. In most cases, 7’p(w) is a constant. For injection lasers, the quasi- 
monochromatic pulse response is the most important contribution. 

In what follows, we assume that the fiber is uniform and long com- 
pared with the period of ray oscillation and therefore approximately 
z-invariant. Let the Hamiltonian in (67) be written 


H =k, — ki(kz, ky, w, x, y) = 0. (80) 
The Hamilton equations (68) are 
dx/dz = —0k,/okz (81a) 
dy/dz = —0k,/dk, (81b) 
dt/dz = dk,/dw (81c) 
dk,/dz = dk,/dx (81d) 
dk,/dz = 0k./dy. (81le) 
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Let us assume further that the medium is isotropic, that is, 


k2 = k?(w, 2, y) — k2 — Ke. (82) 
Thus, (Sla) to (Sle) are 
dx/dz = k./k, (83a) 
dy/dz = k,/kz (83b) 
dt/dz = (0k?/dw)/2k, (83c) 
dk,/dz = (dk?/dx)/2k, (83d) 
dky/dz = (ak?/dy)/2k,. (83e) 


According to (838c), the time of flight of a pulse along a ray for a period 
(period = Z) is obtained by integrating (0k?/dw)/2k, from z = 0 to 
z= Z. If ko(w) = k(w, 0, 0) denotes the wave number on axis, the 
time of flight of a pulse along the z axis is similarly obtained by inte- 
grating (0kj/dw)/2k,. Thus, the ratio of the time of flight of a pulse 
along a ray to the corresponding time on axis is 


t = (ko/kz) (0k? /dw*)/ (dkj/dw?), (84) 


where ( ) denotes an average over a ray period. If the trajectory is not 
periodic, { ) is understood as a limit for zg — «©. When k is proportional 
to w, (84) reduces to 

t = (k?)/kok:z. (85) 


Let us now observe that, from (88a), (83b), (83d), and (838e), 
ghia? (X + Y)/d22 = kh? — k2 + Xok?/0X + Yok?/aY, (86) 


where X = 2, Y = y’. This is easily verified by carrying out the differ- 
entiations. Equations (86), (84), and (79) (with a slightly different 
notation) are those used in the main text. 


APPENDIX B 
Acceptance of Highly Oversized Fibers 


The acceptance, or effective number of modes transmitted by the 
optical system, is the volume of the accepted rays in phase space 
divided by (27)?. We have said earlier that, if the fiber is very long, 
all leaky rays are eliminated and the acceptance is simply the volume 
enclosed by the profile k?(x, y) divided by 4r. If the fiber is highly 
oversized, however, many leaky rays (kz. < k;) are not significantly 
attenuated.1® We need then consider the ray-opties condition that the 
tangential (rather than axial) component of the wave vector be larger 
than k, at the core-cladding interface. The ray-optics acceptance is 
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now evaluated for circularly symmetric fibers. We specify that 
k2 + k2, > ki, a 7 Se (87a) 


where k, denotes the azimuthal wave number at the interface. We 
also have the condition 
kz > 0, (87b) 


which is not implied by (87a). In this appendix, we restrict ourselves 
to small differences in refractive index, in which case condition (87b) 
can be ignored. Because of the conservation of 1, (the axial component 
of the angular momentum), we have 


rk, = ak, (88) 


fora ray with x = r, y = 0, kz, ky, at the input plane, that can reach 
the interface r = a. Thus, condition (87a) is 


k(r) — k2 — k2 + (7?/a®)k2 > k?. (89) 


Equation (89) defines an area in the k,, k, plane bounded by an ellipse. 
We have to make sure, however, that rays outside that area do in 
fact reach the interface. This is not necessarily the case. The maxi- 
mum ray radius rj; is defined implicitly by 


ke + (1 — /ridky = W(r) — (rm), (90) 


where rar is the largest real number that satisfies (90). (The initial 
radius r is considered a constant in the present discussion.) Equation 
(90) shows that the k,, ky that correspond to rj, are contained in an 
ellipse with semi-axes squared k2) = k?(r) — k?(rar) and key = [k?(r) 
—k* (rar) ]/(1 — 1?/ris), respectively. If k?(r) is never increasing, we 
are sure that k.) keeps increasing as rj, increases from 7 to a. We do 
not have any such assurance for kyo, however. When ra, reaches a, 
there may be acceptable values of k,, ky that are located outside the 
ellipse defined above. For each profile, we need therefore verify that 
kio(rm) never exceeds k7y(a). We easily verify that this is the case for 
square-law fibers, because 


ko = ki (rig — 1°)/Q — 9?/ru) = Kir (91) 


increases with ry, for any r. 
Thus, for square-law fibers at least, we can proceed with the calcu- 
lation of the area of the ellipse defined by (89). This area is 


mLk'(r) — ks]. — 1°/a?)-4. (92) 


Substituting this result in the general expression for the acceptance 
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factor, we obtain 
N=) f° Te@) — HIG = /at)-tar (93) 
0 


This expression simplifies if we introduce the variable wu = (1 — 1?/a?)}. 
Equation (93) becomes 


N? = (a2/2) i eC sean. (94) 


Thus, the ray-optics acceptance of most circularly symmetric fibers is 
half the area enclosed by the curve k?(u)a?. For a step-index fiber, we 
obtain from (94) 


N? = (k2 — k2)a?/2  (step-index, ray optics). (95) 


This is twice the wave-optics acceptance. Thus, for step-index fibers, 
the slightly leaky rays carry half the power. Our result agrees with 
that in Ref. 16 for weakly guiding fibers. For a square-law fiber, with 
k(a) = k,, we obtain 


N? = (kj — k2)a?/6  (square-law, ray optics). (96) 


In square-law fibers, 25 percent of the total power is carried by slightly 
leaky rays.* 


APPENDIX C 


Impulse response width of near-square law fibers 


When the source distribution is lambertian, all propagating modes 
are equally excited. It is convenient in that case to express the relative 
time of flight ¢ for near-square-law fibers given in (56) as a function of 
the mode numbers (azimuthal number: p =---—2, —1, 0, 1, 2 --- 
and radial number: a = 0, 1, 2 ---) rather than k, and l,. This can be 
done by quantizing the ray trajectories. [If the W.K.B. method is 
used, it is essential to first remove the singularity of the Helmholtz 
equation at r= 0. This is achieved by changing the independent 
variable from r to log (r)._] One easily finds that the axial component 
of the ray angular momentum lI, is equal to uw. Furthermore, we can 
use for k, the well-known expression applicable to square-law media. 
The result (56) is written below as a function of a, yw, for the reader’s 
convenience. We have 


an) = BL gBD + LPN), (07) 


* This is in agreement with a recent result by D. N. Payne. 
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where 


B = 29K¥?/K, 
=2a+ |u| +1; |u| = abs. val. (x) 
y!2-7LD, — $(v + 1) DiJK,/(KoKY”) 


Ny = (29) SE {2™(y — m) IL (m/2) PILL — (u/g)* (98) 


ea 
I 


The parameters K, = k?, y = 0,1.--- and D,, y = 1, 2 --+ are ob- 
tained from the square of the wave number: K(R) = k?(R) = (w/c)? 
n?(R) of the fiber as a function of R = r?, measured at the nominal 
wavelength > and at a slightly different wavelength, \o, expanded in 
power series of R as follows 


K(R) = Ko — Kih + BR? +--+ (Ao) (99) 
K'(R) = Ky — KiR + KjR? +--+ (Xo). 


The D, are obtained from (99) 
D, = K.(K, — K,)/K,(Ko — Ky). (100) 


If we can neglect the power in the leaky modes, the mode numbers 
a, w are restricted by the condition k, > ks, that is, 


B<1— K./Ko & 2An/n, (101) 


where K, = k} is the square of the cladding wave number. The root- 
mean-square impulse response width is defined as 


o = 5,000L(#) — (¢)?]} ns/km, (102) 


where ( ) denotes an average taken over all the modes permitted by 
(101). Thus, it is a straightforward matter to evaluate from our ex- 
pression in (56) the root-mean-square width of the impulse response of 
any circularly symmetric near-square law fiber, provided the wave- 
number profile can be measured with sufficient accuracy at two 
wavelengths. 
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Effects of Environmental Nuclear 
Radiation on Optical Fibers 


By J. SHAH 
(Manuscript received March 11, 1975) 


A comparison of the available data on environmental radiation and 
on the radiation damage in fibers and glasses under controlled laboratory 
conditions leads us to the conclusion that long-term exposure to gamma rays 
and neutrons in the environment does not pose a serious problem for the 
optical fibers. 


Resistance of optical fibers to damage from long-term exposure to 
environmental nuclear radiation is an important factor to be considered 
when planning a communication system using these fibers. In this 
report, we first summarize the nature and intensity of the natural 
environmental radiation to which the fibers will be exposed, and then 
discuss the available data on radiation damage in optical fibers and 
glasses under controlled laboratory conditions. A comparison of these 
data leads us to conclude that long-term exposure to gamma rays and 
neutrons in the natural environment does not pose a serious problem 
for the optical fibers. 


l. NATURAL ENVIRONMENTAL RADIATION 


The total background radiation at sea level is divided approximately 
equally between extraterrestrial and terrestrial components. The 
extraterrestrial component results from the secondary radiations in- 
duced by cosmic rays, solar radiation, and Van Allen belt radiation. 
The terrestrial component is due to the radiation from naturally occur- 
ring radionuclides in the earth. Gamma rays (y rays) and neutrons (n) 
are important constituents of this radiation? and we will concentrate 
on them for the purposes of this report. 


1.1 Gamma rays 

A number of measurements of the intensity of the environmental 
radiation as a function of location, altitude, and latitude have been 
made. According to Hollaender,? the worldwide average exposure is 
approximately 0.5 R/year (R stands for roentgen, a unit of exposure 


1207 


dose. A brief discussion of units relevant to this report is given in the 
appendix). More recent studies (e.g., Ref. 4) indicate that values 
range from 0.1 to 0.2 rad/year in normal regions (rad is a unit of 
absorbed dose—see the appendix). An average value of 0.13 rad/year 
appears to be generally acceptable’ for normal regions.* However, for 
the purposes of our discussion of radiation damage, we will deliberately 
overestimate the y-ray dose and assume a value of 0.5 rad/year. 


1.2 Neutrons 


Hess et al® have measured the extraterrestrial neutron flux as a 
function of neutron energy and found that the total neutron flux 
[So? o(E)dE] at sea level is % 1.5 X 10° n/cm?-year.t Measurements 
by Herbst’ indicate that the additional neutron flux from terrestrial 
sources is 10° n/cm?-year in open air. However, in tunnels or above 
rocks containing a high density of radioactive nuclides, or in regions 
with high radioactivity, Herbst obtained a flux of up to 4 X 10’ 
n/cm?-year. For the purposes of estimating neutron-induced damage, 
we will assume a rather high value of 1 X 108 n/cm?-year to provide 
us with an extra margin of safety. 


Il. RADIATION DAMAGE IN FIBERS AND GLASSES 
2.1 Gamma rays 


The y rays interact with glasses principally by forcing the electrons 
to leave their normal positions and move through the glass network. 
The primary consequence of this is an increase in the absorption 
coefficient in the vuv-visible-near-Ir range. A detailed study of 
y-induced damage in fibers has been made by G. H. Sigel and co- 
workers’? at the Naval Research Laboratory. They find that the 
y-induced change in the refractive index is small (<10-*) at doses as 
high as 10° rads. They also find that the y-induced losses in optical 
fibers depend strongly on the fiber composition and vary from 1074 
dB/km-rad for bulk Suprasil SiO, to 5 dB/km-rad for Corning fiber 
No. 5010 at 8000 A. Thus, pure fused silica is extremely resistant to 
radiation, while the Corning 5010 is quite susceptible to it. 

A 20-year exposure to natural environmental y radiation (assumed 
to be 0.5 rad/year) would lead to an increase of 50 dB/km for Corning 


* There are regions with exceptionally high level of natural background radiation. 
In some special areas such as Kerala in India or the Santo Spirito province in Brazil, 
values of up to 14 R/year have been reported (Ref. 2). Certain regions of the Black 
Forest (Germany) have shown dose rates up to 1.8 R/year. However, these regions 
are rare and will not concern us in this report. 

t Later reports (see Ref. 7) indicate that the values reported by Hess et al may 
be too high by a factor of 2 to 4. 
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5010 fiber. However, even in this case, the normal bleaching of the 
damage® would probably reduce the total (20-year) loss to $15 dB/km. 
Although this number appears large, it is only a small fraction of the 
loss of 1000 dB/km present in the Corning 5010 fiber before exposure 
to any radiation. Furthermore, this is the worst case reported by 
Sigel et al.° The y-induced losses are generally smaller in fibers with 
smaller initial losses.* For example, the Corning low-loss fiber (type B), 
having germanium-doped silica core and pure silica cladding, has an 
initial loss of 10 dB/km and a y-induced loss of 0.01 dB/km-rad be- 
tween 8000 A and 12,000 A (1.2 um). Thus, even if we neglect bleach- 
ing, the y-induced loss in 20 years would amount to only 0.1 dB/km. 
Since fibers with small initial losses are precisely the ones that will be 
used in communication systems, it seems reasonable to conclude that 
long-term exposure to environmental y radiation will not seriously 
affect the fiber performance. 


2.2 Neutrons 


Neutrons interact principally with the nuclei rather than electrons 
in solids. Neutron radiation, therefore, results not only in increased 
absorption losses but also in structural changes that lead to changes in 
density, refractive index, rotary power, birefringence, thermal con- 
ductivity etc. Since small differences in refractive indices of the core 
and the cladding are essential to fiber performance, we will pay par- 
ticular attention to refractive index changes as well as to increased 
losses caused by n-irradiation. 

To our knowledge, the only study of n-induced losses in optical 
fibers is by Maurer et al.!° They irradiated high-silica-glass multimode 
fiber waveguides with 14-MeV neutrons, using doses of as high as 
1.4 X 10” n/cm?. They concluded that the n-induced loss varies 
roughly linearly with the total dose and is less than 1.5 X 10-! 
(dB/km)/(n/em?) in the 8000-A to 12,000-A region. This number, 
which is obtained from the figure given by Maurer et al,’ is in fact an 
overestimate of n-induced damage, because we have disregarded the 
fact that the n-irradiated samples also received a simultaneous dose 
= 1000 rads of y radiation. However, even if we assume this to be the 
true value, a 20-year exposure to environmental n-irradiation (2 X 10° 
n/cm?) would increase the loss by only about 3 X 10-? dB/km. It 
should be emphasized that this extrapolation is only approximately 
valid because neutrons in the environment have a wide range of 
energies (from 0.01 eV to 10° eV), whereas the neutrons in the con- 


* While there is no evidence that the correlation between low radiation damage 
and low initial losses is universally valid, such a correlation definitely exists in the 
presently available data. 
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trolled experiment were monoenergetic (14 MeV). However, even 
after a thirtyfold increase, the n-induced losses would still be less than 
1 dB/km. Therefore, it seems reasonable to conclude that absorption 
losses induced by long-term exposure to environmental n-radiation 
will not seriously affect fiber performance. 

Neutron-induced changes in the refractive index of the fibers can 
be a potential source of problems. We know of no measurements on 
fibers which can shed light on this problem. However, an extensive 
literature exists on the effects of n-irradiation on various forms of 
silica and other commonly used glasses (a good summary is given in 
Ref. 11). The refractive index of vitreous silica changes by 0.67 percent 
under a flux of 2 K 10” n/cm? of thermal (<0.1 eV) neutrons.” From 
the measurement by Primak,!? we deduce that the rate of increase of 
the refractive index of vitreous silica is approximately 5 X 10-” per 
(n/em?) for doses less than 1 X 10 n/cm?. This suggests that the 
changes in refractive index induced by environmental neutrons 
(2 X 10° n/cm? in twenty years) will be less than 1 X 107”, a truly 
negligible effect when we consider the fact that the difference in the 
refractive index of the core and the cladding is typically larger than 
107°. 

No data are available on the n-induced changes in refractive indices 
of other glasses. However, density changes have been investigated 
for many glasses."' For vitreous silica,” the density increases approxi- 
mately linearly (10- percent per n/cm?) up to 2.5 X 10" n/em? and 
then saturates. Other glasses (except borosilicate glasses) are also 
quite resistant to neutrons and show very few changes up to about 
10!7-10'8 n/cm?.* The borosilicate glasses are more susceptible because 
boron, like other light elements, has high neutron cross section. How- 
ever, even these glasses show damage only when flux levels exceed 
104 n/em?,* which is some five orders of magnitude larger than the 
accumulated (20 years) flux of 2 X 10° n/cm? encountered in the 
environment. 


Ill. CONCLUSIONS 


We have summarized the available data on environmental nuclear 
radiation and also the data on radiation damage in glasses under con- 
trolled laboratory conditions. Unfortunately, the laboratory experi- 
ments were not performed with the exact y ray or neutron energy 
distributions that one encounters in environmental radiation. It is 
difficult, therefore, to make accurate predictions about the radiation 
damage in fibers caused by environmental radiation. However, we 


*See Table 6.12 in Ref. 10. 
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have made some approximate estimates from which it is reasonable 
to conclude that damage induced by environmental y or neutron 
radiation should not pose any serious problems to optical fibers so far 
as their optical loss or refractive index are concerned. More recent 
experimental works by Evans and Sigel!’ and Mattern et al. do not 
affect this conclusion. 

Some general comments seem to be appropriate in conclusion. Pure 
fused silica seems to be extremely resistant to radiation damage. It is 
also useful to remember that the addition of small amounts (0.1 to 
0.2 percent) of Cerium?! makes most glasses more resistant to radia- 
tion. We have not discussed damage by ea particles, but it is appro- 
priate to mention here that a particles have very short ranges in air 
as well as in most other materials. Therefore, it seems unlikely that a 
particles will pose any problems for the optical fibers if the fibers are 
enclosed in a conduit. Finally, the background luminescence induced 
by environmental ionizing radiation has been considered by Cohen 
and Lanzerotti!® and found to be not significant for fiber optic com- 
munications systems. 


IV. ACKNOWLEDGMENTS 


We acknowledge helpful discussions with G. H. Sigel, Jr. and M. 
Weiss. 


APPENDIX 
Units of Dosimetry 


The most useful units in the study of radiation damage in solids 
are the particle or photon fluxes as a function of their energy. Thus, 
¢(£)dH, expressing the number of particles/em?-sec in the energy 
range LH to # + dH, completely specifies the incident radiation field. 
However, many special units are frequently used in specifying the 
radiation. Roentgen (R) is a unit of exposure dose used for X rays 
and y rays and is defined as follows. Roentgen is that exposure of X 
or y radiation which gives a dose of 87.7 ergs/g to air. 

A special unit of absorbed dose is called a “rad.’”’ One rad = 100 
ergs/g. 

Unlike the roentgen, the rad is independent of the irradiated ma- 
terial. This means that a given beam of radiation acting for the same 
time will deliver different doses, expressed in rads, according to whether 
it is absorbed in air, tissue, or other materials. The rad in Section I 
refers to air as the reference material. The rad as used here is indirectly 
a measure of the radiation field rather than the absorbed dose in the 
sample because it refers to energy absorbed by air rather than the 
sample under study. Under these conditions, rad and roentgen are 
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numerically equal within about 20 percent and can be used 
interchangeably. 

In the work reported by Sigel et al.,° Si is used as the reference 
absorbing material. The differences in using air or Si as the reference 
material are small (less than a factor of two) and are inconsequential 
for the purposes of this report. 

The conversion between rad and n/cm? and photons/cm? as a func- 
tion of energy are given by H. Stern.!® (See also the report by J. 
Moteff.!7) For example, for 1 MeV y-ray photons, 1 rad} 2 X 10° 
photons/em?. For 1 MeV neutrons, 1 rad } 2.6 X 108 n/em?. For 
y rays with energy E between 0.07 and 2 MeV, 1 rad (air) } 2 X 10°/E 
photon/cm?. Conversion factors at other energies may be obtained 
from the above references. See also the American Institute of Physics 
Handbook.!8 
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Loose Tube Splices for Optical Fibers 


By C. M. MILLER 
(Manuscript received March 14, 1975) 


A technique for splicing optical fibers has been developed that uses a 
self-aligning square cross-section tube, with winner dimensions slightly 
larger than the optical fiber. A total loss of 0.58 dB was obtained for eight 
splices connected in series using a graded-index fiber with a 68-pm core 
diameter. The splices were made one at a time without the use of micro- 
scopes or micromantpulators; however, the fabrication process could be 
mechanized and extended to groups of fibers. A holding fixture could be 
added to adapt this technique to a connect-disconnect type splice. The 
size of the splice is presently 0.012 in. square, making zt surtable for use 
within cables. Measurement set refinements that were needed to measure 
endwidual splice losses as low as 0.05 dB include an improved detector 
and means for better control of launching conditions. 


]. INTRODUCTION 


The basic requirements for low-loss splices are (2) accurate align- 
ment, (72) good fiber ends, and (777) accurate diameter control. Trans- 
verse alignment accuracy of approximately -++0.1-fiber-core radius 
(typically, +0.0001 in.) is required to achieve a splice loss of 0.1 dB. 
Good fiber ends may be produced by scoring and breaking,! grinding 
and polishing, or dise sawing. Accurate fiber diameter control is also 
needed; however, significant progress is being made in this area. Of 
these three requirements, accurate transverse alignment may be the 
most difficult problem to solve, especially when the field environment 
and variability of craftsmen’s skill are considered. 

Single-fiber splicing has been accomplished by Bisbee? and Dyott 
et al? using heat fusion. Someda‘ suggested using embossed plastic to 
obtain transverse alignment. This paper describes a splicing technique 
that uses a loose-fitting, square, cross-section tube to align the fibers. 
The splices produced are small, exhibit very low losses, and are simple 
and inexpensive. 

Previously, snug-fitting sleeves have been suggested,® but three 
problems are usually encountered. 
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(2) 


(1) 


(217) 


The 


If the glass sleeve is to support the fiber with the required 
alignment accuracy, it must be less than typically 0.0001 in. 
larger than the fiber. Both fiber and sleeve must be highly circu- 
lar, and the fiber diameter must be controlled to at least the 
same tolerances. These tolerance requirements have discouraged 
efforts to use a snug tube as an alignment mechanism. 

Given a snug tube of the proper dimensions, the initial insertion 
of a fiber into that tube is difficult. Pinnow® has described a 
method of flaring the inner diameter of capillary tubes, which 
reduces the initial insertion difficulty. 

Contaminants that are scraped off the inside wall of snug- 
fitting sleeves during fiber insertion are trapped between the 
fiber ends where the effect of contamination is worst. 


“‘loose’’-fitting square-tube splice described below reduces these 


difficulties substantially and appears to have potential application in 
several places in a fiber-optic communication system. 


Il. SPLICE CONFIGURATION AND ASSEMBLY 


The 


loose-tube splice combines the alignment accuracy obtainable 


by using a groove for alignment with the small size and simplicity of 
glass sleeves.> The fiber ends are biased to one corner of the square 
cross section by bending the fiber outside the tube. Figure 1 is a 
pictorial layout of the square tube with two fiber ends in position 


within 
radius 


the tube. The tube has nearly flat interior walls and a small 
in the interior corners, as shown in the cross section in Fig. 2. 


One corner of the square is used as a groove for aligning the fibers. 
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Fig. 1—Splice configuration. 
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Fig. 2—Square tube cross section. 


Epoxy is forced into the square tube prior to insertion of the fibers 
and serves several useful functions. 


(t) The epoxy serves as an adhesive after curing to hold the splice 
together. 
(21) The epoxy also serves as an index-matching material with good 
glass-wetting characteristics. 
(227) Contamination on the fiber ends is washed away by the flow 
of epoxy around the fiber ends during insertion of the fibers. 


Assembly of a splice involves inserting two fibers with good ends 
approximately halfway into each end of a square cross-section tube 
filled with uncured epoxy. No particular orientation of the square- 
tube cross section is required. The fibers are placed on a flat surface 
and bent in a curved pattern. This causes forces to be generated that 
rotate the tube so that a diagonal of the square cross section is in the 
same plane as the bent fibers. The tube is therefore self-aligned and 
the fibers biased to one corner by action of the fiber stiffness. After 
the bends are made, the fibers are taped to a flat surface in the bent 
configuration and the fibers pushed into the tube until they touch each 
other. Figure 3 is a cross-section photograph of a splice showing a 
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Fig. 3—Splice cross section showing position of fiber in vicinity of joint. 


fiber in a corner of the square tube. Figure 4 is a magnified view of one 
splice, and a longitudinal section is shown in Fig. 5. In spite of the 
small angle between the fiber ends caused by one end not being broken 





Fig. 4—Single loose tube splice, with tube approximately 0.5-inch long and 0.012- 
inch wide. 
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Fig. 5—Longitudinal section, 250X. 


perpendicular to the fiber axis, the splice loss for the splice in Tig. 5 
was only 0.07 dB. 


Ill. END PREPARATION 


Good fiber ends are necessary for the fabrication of low-loss splices. 
As mentioned earlier, several techniques exist for preparing suitable 
ends. A score-and-break technique! was used for end preparation on 
all splices reported in this paper. A single fiber is clamped in the 
apparatus shown in Tig. 6 with approximately 100-g load applied to the 
fiber by a spring. The fiber rests in a groove along a 2-in. radius are 
and is scored lightly with a hand-held diamond knife edge. Fiber ends 
prepared by this method are nearly perfect, as shown in Fig. 7. A very 
small amount of edge chipping is present where the fiber was scored. 


IV. SPLICE-LOSS MEASUREMENTS 


After constructing just a few square-tube splices, it became evident 
that the measurement set-up being used was not adequate for losses 
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Fig. 7—Good fiber end. 
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Fig. 8—Single fiber detector. 


below 0.1 dB. A new detector was built, which repeated to within 
+0.015 dB. Accurate positional repeatability was accomplished by 
using a 0.05-in. inner-diameter capillary tube which was tapered to 0.006 
in. inner diameter by heating and pulling. The fiber is easy to insert in 
this detector, and the positional repeatability is excellent. The solar cell 
sensor was Immersed in index-matching liquid to reduce reflections, 
and a microscope cover slip was used to protect the cell from damage 
by the fiber being measured. A photograph of the detector is shown 
in Fig. 8. At the input end of the measurement set, a vacuum chuck 
was added to ensure repeatability of launching conditions. This chuck 
positions the fiber accurately along the optical axis of the 30X launch- 
ing objective lens and the laser. The overall repeatability of splice- 
loss measurements is within +0.03 dB. 

Care was taken to place the fiber in the same coiled configuration 
after splicing so that bending losses before and after splicing would 
be similar. Fiber loss of approximately 0.01 dB/m was subtracted from 
the total loss measurements so that losses stated are for the splices 
only. All loss measurements were made at a wavelength of 0.6328 um. 


V. RESULTS 


Initially, a fixture was used to hold the square tube and control the 
fiber bending. Although the maximum loss measured on 25 consecutive 
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splices fabricated with this fixture was 0.21 dB, it was found that the 
tube would align itself if allowed to rotate. Losses were lower and the 
assembly of the splice, as previously mentioned, was much easier. 

Eight epoxied square-tube splices fabricated in series using no 
fixtures except the end-making apparatus produced a total splice loss 
of 0.58 dB or 0.073 dB per splice. The splices had approximately 14 m 
of fiber between each splice and on each end and were put together 
in series to increase the total loss to an accurately measurable quantity. 

It has been found that splices measured in series have higher losses 
than when measured individually. Ten earlier splices were fabricated 
in series and measured 1.37 dB. The ten-series splices were measured 
separately and gave the distribution shown in Fig. 9. The average loss 
per splice was 0.077 dB compared to 0.137 dB in series, or a nearly 
2-to-1 increase for splices in series. The process of peaking up the 
power through a single splice probably selects the launching conditions 
and therefore the mode structure best suited to the particular imperfec- 
tions of that splice. A loss measurement made in this way gives a value 
that is too low. That is, the loss of a splice with long lengths on either 
side or with other splices nearby is apt to be considerably higher than 
when measured separately with short fibers on each side of the splice. 
The eight-series splices mentioned earlier were not measured separately 
because losses as low as 0.03 dB, which would be expected based on the 
series loss, could not be accurately measured. 

A slight longitudinal separation of fiber ends within the tube occurred 
during epoxy cure for the 10 splices in series. Loss measured before 
epoxy cure was 0.69 dB and 1.37 dB after epoxy cure as stated above. 
Iibers were bent in a 90-degree are and taped to an optical table while 
the epoxy cured. This configuration applied very little, if any, force 
component to hold the fibers in place during cure. The eight-series 
splices were bent through an are of approximately 45 degrees and 
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Fig. 9—Splice loss histogram. 
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taped to the table with a definite force applied to keep the fiber ends in 
contact. Although this force (a result of the fiber stiffness) is small, end 
separation did not occur and a loss of 0.58 dB was measured before 
and after epoxy cure. 

These extremely encouraging results have stimulated thinking as to 
how loose square-tube splices could be applied to other types of 
splices, e.g., connect-disconnect configuration. 


VI. CONNECT-DISCONNECT SPLICE FOR SINGLE FIBERS 


Several configurations based on the square tube can be envisioned 
for a connect-disconnect splice, that is, a splice that can be reassembled 
many times and used as a connector. Figure 10 is a photograph of a 
simple fixture that supports a single square-tube splice by the fibers 
on each side of the splice. The splice itself is suspended in air. This 
fixture is not intended to be a finished design, but it does produce 
losses of 0.1 dB or less. The clamps are lined with a thick, soft EVA 
layer that grips the fibers and holds them in position. More practical 
designs are sketched in Fig. 11. An index-matching material is neces- 
sary to achieve 0.1 dB, although a liquid index-matching material 
may be feasible. 

Polymethylmethacrylate (PMMA) was suggested as an alternative 
index-matching material by Pinnow.® This thermoplastic could be 
drawn into a fiber and inserted into the square tube. Heat could then 
be applied to melt the PMMA, the fibers inserted, and the PMMA 
allowed to cool. The splice could be disassembled by again heating the 
tube and removing the fibers. These connect-disconnect splices may 
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Fig. 10—Quick-connect holding fixture. 
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Fig. 11—Connect-disconnect holders for single-fiber splices. 


be useful as methods for connecting sources, detectors, and line 
regenerators in a fiber-optic communication system. 
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The Effect of Longitudinal Imbalance 
on Crosstalk 


By G. MILLER 


(Manuscript received December 30, 1974) 


Multipair cables are carefully constructed to control the coupling be- 
tween wire pairs to prevent excessive crosstalk. Several types of coupling 
modes exist between pairs, but the principal effort is devoted to controlling 
the ‘‘metallic-to-metallic’”’ coupling mode because the coupling loss in this 
path zs the most important in ordinary telephonic use of the cable. Over a 
half-century ago, Campbell showed that crosstalk behavior of this mode 
could be characterized by measuring a simple function called capacitance 
unbalance. This paper shows that at voice frequencies the crosstalk charac- 
teristics of the “longttudinal-to-metallic’ mode can be predicted by 
measuring a second similar function of the same parameters that are con- 
tained in the capacitance unbalance. With the aid of these two functions, 
it ts shown how the longitudinal balance of terminal equipment connected 
to a cable pair affects crosstalk. It is further shown that a longitudinal 
balance of approximately 40 dB or more ts necessary for any station or 
terminal equipment used in the telephone network so that it will not 
significantly increase the small amounts of crosstalk inherent in the careful 
cable design. Also, a limitation is established for the maximum longitudinal 
voltages at voice frequency that can be applied without noticeably increasing 
crosstalk and noise in other cable pairs. This limitation 1s approximately 
40 dB more restrictive than the tariff limitations for metallic voltages. 


Il. INTRODUCTION 


A multipair cable consists of many insulated but unshielded con- 
ductors within a protective conducting sheath. The individual con- 
ductors are used to form circuits. In one configuration, called a metallic 
circuit, two conductors are paired and form the circuit. Signals are 
applied between them. This is called metallic excitation of the circuit, 
and the signal is said to propagate in the metallic mode. In another 
configuration, called a longitudinal circuit, two conductors are paral- 
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leled and these, plus the conducting cable sheath, form the circuit. 
Signals are applied between the paralleled conductors and the conduct- 
ing cable sheath, which is grounded. This is called longitudinal excita- 
tion of the circuit, and the signal is said to propagate in the longitudinal 
mode. It is also possible for one wire pair to be used for both circuits 
simultaneously and, consequently, for one wire pair to be simultan- 
eously excited in the metallic and longitudinal modes of propagation. 
This happens when the terminal equipment is longitudinally un- 
balanced, as will be explained. 

Because the conductors are not shielded and are in close proximity 
to each other, electromagnetic fields generated by current flowing 
through the conductors cause energy to be coupled from one circuit 
to another. This is called crosstalk and is undesirable, since it may 
cause noise in other circuits that can impair the performance of digital 
and analog systems, or even be intelligible speech that is overheard 
and leads to loss of privacy. 

Crosstalk cannot be eliminated, but several things can be done to 
reduce it, that is, to increase the crosstalk loss between circuits. First, 
metallic circuits are used rather than longitudinal circuits, because it 
was found by experience that the crosstalk loss between two metallic 
circuits is generally greater than the loss between two longitudinal 
circuits or between a longitudinal and a metallic circuit. Second, ad- 
jacent conductors are paired and often twisted and are used for the 
metallic circuits because they are less susceptible to inductive noise 
and the crosstalk loss between twisted pairs is generally greater than 
between nontwisted pairs. Twisting reduces crosstalk by assuring that 
each pair of the cable is exposed to opposing couplings by transposing 
its conductors relative to the disturbing pair. Third, the terminal 
equipment at both ends of a pair should be longitudinally balanced, 
i.e., have impedance symmetry with respect to ground, because longi- 
tudinal imbalance has the effect of producing longitudinal excitation 
which consequently can increase crosstalk. Finally, the cable pairs are 
also constructed to have longitudinal impedance symmetry for the 
same reason. 

Since metallic circuits are usually used, and both cables and terminal 
equipment are usually constructed to be longitudinally balanced, most 
crosstalk studies to date have concentrated on what is called metallic- 
to-metallic crosstalk, i.e., crosstalk between balanced metallic circuits. 
Much less is known analytically about the crosstalk loss between 
longitudinal circuits or between a longitudinal and a metallic circuit. 
For example, to explain crosstalk between balanced metallic circuits, 
Campbell! assumed that all circuits within a cable were longitudinally 
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symmetrical, that the pairs were excited metallically, and that, con- 
sequently, the applied metallic signal would not excite any longitudinal 
voltage in the disturbing or disturbed pairs. Thus, crosstalk would be 
due to cable characteristics alone. Campbell was then able to show 
that crosstalk at low or voice frequencies, where inductive coupling is 
negligible, was very nearly proportional to the capacitance unbalance, 
which is a function of the four interwire capacitances between two 
cable pairs, and is now used as a measure of quality of a cable with 
regard to crosstalk performance. 

In another study, Foschini? developed an accurate transmission 
model of cable systems for computing crosstalk which is an extension 
of Campbell’s work. He too assumed longitudinal symmetry and 
showed that crosstalk coupling losses between metallic circuits can be 
predicted quite accurately from Campbell’s capacitance unbalance. AI- 
though his results are valuable in the study of crosstalk for the metallic 
mode of propagation, they, as well as Campbell’s results, do not con- 
sider the effects of terminal imbalance on crosstalk loss. 

The objectives of this paper are to extend the results of Campbell 
and Foschini by first removing the constraints of metallic circuits, 
terminal balance, and pair symmetry; and to construct a model to 
permit calculating the crosstalk loss between pairs as a function of 
terminal balance and pair symmetry. The model is used to show why 
terminal imbalance can greatly increase crosstalk by causing longi- 
tudinal excitation of a cable pair and, consequently, why limitations 
must be imposed on the longitudinal balance of terminal equipment 
and on the direct application of longitudinal signals. These objectives 
are accomplished by showing, through numerical solutions and ex- 
perimental results, that longitudinal excitation couples energy into 
adjacent wire pairs with much less loss than does metallic excitation. 

The paper is divided into four sections. First, the important results 
on longitudinal balance and longitudinal voltage restrictions are given. 
Next, the model of crosstalk between two wire pairs in a cable is 
analyzed using transmission line equations. This model is used to 
derive a new set of crosstalk coupling coefficients that can be used to 
relate the crosstalk loss between two metallic circuits, a longitudinal 
and a metallic circuit, and two longitudinal circuits. Third, average 
values for these coupling coefficients for a typical cable are obtained, 
derived from measured characteristics. Using the coupling coefficients 
the predicted increase in crosstalk resulting from longitudinal excita- 
tion is compared with direct measurements of the increase made on 
another cable. Finally, restrictions on longitudinal balance and longi- 
tudinal voltages are established. 
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ll. RESULTS 
2.1 Requirement on longitudinal balance 


For application to crosstalk performance, it is appropriate to define 
a longitudinal balance* of terminal equipment as 


em(f) 
ex(f) 


where ¢é, is the longitudinal voltage produced when a metallic voltage 
em is applied at any frequency f. The subscript ‘““mM — L”’ means the 
conversion from a metallically applied voltage to a longitudinal voltage. 
This paper shows that a balance of approximately 40 dB or more in 
the voice frequency region is required for any terminal device to ensure 
that the level of crosstalk that already exists in the network will not 
be significantly increased. This requirement is based on measurements 
of the near-end crosstalk at 1000 Hz of cable with a balanced and 
unbalanced termination. It is assumed that any metallic signal applied 
to the telephone network does not exceed the power level specified in 
Ref. 3. Longitudinal and metallic voltages are defined in Section 3.3. 





BALm_1(f) = 20 log io 





2.2 Restriction of longitudinally applied voltages 


Crosstalk coupling losses decrease with increasing frequency and 
hence voltage restrictions are frequency dependent. Figure 1 shows 
the limitations on applied longitudinal voltages established so as to 
increase the crosstalk energy already present in the telephone network 
by no more than about 1 dB. 


2.3 Derivation of crosstalk coupling coefficients 


Three new capacitive coupling coefficients have been derived that 
can be used with a simple but reliable computation method to predict 
the degradation in crosstalk performance for a particular cable when 
any of its terminations are unbalanced. These coefficients are given in 
Table I. The coefficients are defined in eqs. (10) to (13), and the 
interpair capacitances given in the formulas are the capacitances be- 
tween the pairs shown in Fig. 2. 


lt. COUPLING BETWEEN TWO WIRE PAIRS IN A CABLE 
3.1 Transmission line model 


Figure 2 models two wire pairs within a cable of length 1. The follow- 
ing assumptions about a cable are made to construct this model: 


* A second type of balance for noise immunity purposes is a separate but important 
consideration for good telephone network performance. It is defined in Section 4.2. 
However, crosstalk does not enter into establishing its restrictions. 
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(c) The impedance and admittance per unit length of each wire 
pair, the admittances to ground per unit length, etc., are 
constant. 

(<7) The conductance to ground and between wire pairs is negligible, 
i.e., the admittances are purely capacitive. 

(222) The inductive coupling between pairs is negligible at voice 
frequencies. 

(iv) The impedances per unit length of all wire pairs are equal. 


Let the admittances per unit length between the two circuits be Y43, 
Yo3, Yos, and Yi, connected between conductors | and 3, 3 and 2, 2 and 
4, and 4 and 1, respectively, where conductors 1-2 form one twisted 
wire pair and 3-4 form the other pair. The impedances per unit length 
of the four wires are Z1, Z2, Z3, and Z, and the admittances of wires 
to ground are Yy,, Yo,, Ys,, and Y4,. The admittances per unit length 
of the wire pairs are Yy. and Y3,, and the voltages and currents are 
labeled in the figure. 

Consider a differential section of the model of length Az. It is readily 
seen that the following eight current-voltage relationships hold for 
this differential section: 


Vila + Ax) = Vila) — Lh(x)Z Ax (1a) 
Vo(a + Az) = Vo(x) — I2(x)Z2Ax (1b) 
V3(a2 + Ax) = V3(a) — I3(x)Z3Ax (1c) 
Vila + Ax) = Vala) — Iu(2)ZsAx (1d) 


Ii(a + Ax) = Iy(a) — {Vi(a) — Vo2(2)} Yi2dx 
— V,(z) Y,,Ax = {Vi(2) = V3(a)} YisAx 
= {Vi(2) = Va(x)} VigAx (le) 


In(a + Ax) = Io(a) — {Veo(x) — Vilx)} Virda 
— V2(2) Y2,Ax = { V2(2) = V3(x)} YosAx 


I3(a + Ax) => I3(2) = {V3(x) = Va(x)} Vasa 
—V3(2) Y;,Ax = {V3(2) = Vo(x)} YosAx 
—{V3(2) — Vi(x)}VYy3Ax (1g) 
I(x + Ax) = In(x) — {Va(x) — Va(x)} YasAr 
= V4(x) Y4,Ax = {Va(z) = V2(a) } Yo,Ax 
Dividing through by Az, taking the limit as Az approaches zero and 


recognizing the definition of the derivative, using assumption 3, and 
writing the resulting eight equations in matrix form, we obtain eq. (2). 
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ceel 


SZ6l YAEWALdSS “IWNYNOP WOINHOSL WALSAS 1139 SHL 


dV ,/dx 
dV 2/dx 
dV3/dx 
dV 4/dx 
dI,/dx 


dI2/dx 


dI3/dx 


dl s/dx 


0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 
= (V1, + Yio Yie Yi3 Yu4 
+¥ia+ Vu) 

Yie oi (Yo, + Yie Yos Yo 

+ Yost Yau) | 
Yis Yo3 = (Y3, =e V4 Y 34 
+ Yi3 + Y 93) 
Yu You Y34 "Ty (Yas + Y34 
+¥ut Yo) 


Vi 
Vo 
V3 
Va 
I, 


I. 


I; 


Is 
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Fig. 1—Longitudinal signal limitations. 
3.2 Transmission line equations 
Equation (2) above can be written in matrix notation as 
dV 
— = —ZI 3a 
Te (3a) 
dl 
—= —YV. 3b 
dx ? ( ) 
where 
Vi(x) Zi 0 0 0 I, (x) 
Va(a) |’ 0 O 4; Of’ I3(x) 
V(x) 0 0 O Zs I4(x) 


Table |— Average crosstalk coupling coefficients 
for a multipair cable 





Average 
Coefficient Magnitude Formula 
(picofarads) 
Cum 7.5 Cis — Cys — Coz + Cog 
Cran 65.2 2{Cis + Cre — Crs — Cra} 
Crum 64.0 2{C13 - Cis + Co3 — Cx4} 
Crilh 9076.0 A{Cis + Cig + Coz + Cog} 
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Fig. 2—Model of wire pairs and terminations used to analyze crosstalk coupling. 


are the voltage, impedance, and current matrices, respectively, and, 
using assumption (22), 


C1 + Ce —C pe —Cy3 —C 4 
+C13 ie Cis 
—C ye Cog + Cie —C 3 —Cr 
Y = jw +C 23 + Co 
—C; —C 23 C's, + C34 —C34 
+C43 i C23 
—C 14 —Co4 C4, a lg C34 
+C i i C24 


is the admittance matrix. Equations (3a) and (3b) are basic transmis- 
sion line equations describing the voltage-current relationships be- 
tween wire pairs. We will use them to calculate crosstalk coupling 
between wire pairs. They are more conveniently written as a matrix 
differential equation: 


alr}-~[y ollr} n 


where O is a 4 X 4 null matrix. The solution of eq. (4) is straight- 
forward and is discussed in Appendix A. 


3.3 Crosstalk coefficients for various coupling modes 


Since there is negligible inductive coupling at voice frequencies, the 
insight to crosstalk coupling can be obtained from the transmission 
line equations involving the admittance matrix only. Rewriting eq. 
(8b) explicitly, we have 


dI;/dx Cig + Cr —Che —C 13 —Cu Vi 
+C13 + Cus 
dI./dz —Cip Cog + Cre —Co3 —C x V2 
; +C23 + Cr 
== — jw © 
dI3/dzx —C3 —C3 C39 + C34 —Cx V; 
+C13 + C23 
al ,/dz —Cy —C4 —Cu C19 + C's4 V4 
+Cu + Co 
(5) 


It is clear from eq. (5) that coupling between wire pairs 1-2 and 3-4 
could not possibly occur if the four interpair capacitances Cy3, Cy4, 
C23, and C's, were all zero regardless of the longitudinal imbalance at 
the terminations. Furthermore, inspection of eq. (5) shows that the 
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coupling between the wire pairs is some function of these four inter- 
wire capacitances. This is seen by observing the contributions to dJ,/dx 
and dI,/dx from V3 and V4, and the contributions to dI3/dx and 
dI,/dx from V, and V2. 

The insight needed to understand why longitudinal excitation and 
longitudinally unbalanced terminations increase crosstalk is obtained 
when eq. (5) is transformed and expressed in terms of the longitudinal 
and metallic voltages and currents, rather than in terms of the con- 
ductor currents and conductor-to-ground voltages. This transforma- 
tion is easily made because the longitudinal and metallic voltages and 
currents are linearly related to the conductor voltages and currents. 
If wire pair 1-2 is now denoted as circuit one and wire pair 3-4 is 
denoted as circuit two, then the metallic voltages and currents on the 
two circuits are defined to be 

I, a I, 


Vim = Vi—- Vs, Li ge 


and 
Is a Ia 


Vom = Vs —_ V4, Lom = 5 


The longitudinal voltages and currents on the two circuits are 


Vp = Vat Vs 


5 j Tin = I, + L2 
and 
Vv 
Var = — Ion = Is + Ih. 


Expressed in matrix form, these eight equations become 


Vi % 1 O O}l Via 
Ve} _ |-% 1 0 O||Viz 6) 
V3 0 0 5 1 Vom 
Va 0 0 —4 1 Vor 
and 
Tu 2-3 0 O} lL; 
Tin} _ {1 1 0 Of |de] 
Ieu| |0 O 4 —4| [Zs iy) 
Ton 0 01 41)U, 


Now, by using eqs. (6) and (7), eq. (5) can be expressed in terms of 
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the longitudinal and metallic voltages and currents as 





dT ya¢/dx 3-2 0 0 
al i,/dx 1 1 0 O 
= — jw 
dT oa¢/dx 0 O } -3 
dIy1/dx 0 01 1 
Cig t Cie —Cpe —Cis —Ci4 z 1 O O) {Viar 
+O 3+C 4 
—Che Cog t Cre —Co3 —C 4 —} 1 0 0 Vit 
+Co3+Co4 
x ) 
—C13 —C2x C3g C4 —C34 0 0 —} 1 Vom 
+CistCos 
—Ci4 —Co —Cs4 Cig tC 00 43 1 Vor 
+CyutCau J 
(8) 
or, performing the matrix multiplications, 
dln /dx Qiao iz Qs 4 Vim 
dIy,/dx - _ jw Gor G22 Qes Goa} | Viz (9) 
dT on /dx 4 Q31 «332 33) 34 Vom 
dIo1/dx G41 G42 Gaz 4a Vor 


where the elements of the 4 X 4 symmetric matrix Q in eq. (9) are 
given in Appendix B. 

Much useful information can be obtained by simple inspection of 
some elements of Q. First, the derivative of the metallic current in 
circuit one due to the metallic voltage in circuit two is proportional 
to quis, 1.e., the crosstalk coupling loss between the metallic circuits is 
directly related to qi3. Thus, the coupling between two metallic circuits, 
Le., the metallic-to-metallic coupling, is proportional to 


Cum, = —Q13 = Cis — Cr — Cos + Cas. (10) 


This is the capacitance unbalance term first derived by Campbell! and 
used today as one measure of cable quality. Referring again to eq. (9), 
we see that the derivative of the metallic current in circuit one due 
to the longitudinal voltage in circuit two is proportional to qi, and that 
the derivative of the metallic current in circuit two due to the longi- 
tudinal voltage in circuit one is proportional to q32. In other words, the 
crosstalk coupling from a longitudinal to a metallic circuit is pro- 
portional to 

Crim, = —Qu = 2(Cis + Cis — Coz — Co) (11) 
or 


Crim, == 932) = 2(Cis — Cy + Co3 — C24). (12) 
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The subscript Z.M/, means “from the longitudinal mode in circuit 
two to the metallic mode in circuit one.’’ Also, we can readily see that 
the derivative of longitudinal current in circuit one resulting from the 
longitudinal voltage in circuit two is proportional to qz4. In other words, 
the crosstalk coupling between two longitudinal circuits is proportional 
to 


Crt, = — Gea = 4(Ciz + Cis + Cas + Cay). (13) 


Using these four coupling coefficients, it is now possible to compare 
the difference in crosstalk loss between two metallic circuits, a longi- 
tudinal and a metallic circuit, and two longitudinal circuits. This com- 
parison was made for one cable and the results are discussed in the 
next section. 


3.4 Comparison of crosstalk using the coupling coefficients 


One good feature of the four coupling coefficients given in eqs. (10) 
through (13) is that they are easily measured. Hence, they provide a 
simple method for comparing the difference in crosstalk between two 
metallic circuits, a longitudinal and a metallic circuit, and two longi- 
tudinal circuits. To make such a comparison, it is necessary to have 
data on the interwire capacitances, C13, C14, Ce3, and Cas, for real cable. 
Such measurements were made in 1968 on a 22-gauge, pulp-insulated 
cable manufactured by Western Electric. These measurements were 
made on many different 50-pair binder groups.* The data on inter- 
Wire capacitances were taken for random samples out of the 1225 
possible sets' of interwire pair combinations within each binder group. 

Using these data, the average value of the four coupling coefficients 
were calculated and are given in Table I. These show that, on the 
average, the coupling between two metallic circuits is significantly 
less than the coupling between a longitudinal and a metallic circuit, 
and that the coupling between two longitudinal circuits is by far the 
greatest. Hence, the fundamental reason why terminal longitudinal 
imbalance increases crosstalk is that longitudinal imbalance causes excita- 
tion of the longitudinal circutt. 

Comparison of the values of the coupling coefficients made so far 
does not provide any quantitative estimate of the amount of the 
differences in crosstalk losses to be expected. Such an estimate can be 
obtained by using the coupling coefficients for individual wire-pair 
combinations to construct distributions of 1000-Hz near-end crosstalk 


“A binder group is a unit of 12, 16, 20, 25, 50, or 100 twisted wire pairs bound 
together within a cable. 

T For a 50-wire pair cable there are n(n — 1)/2 = 50(49)/2 = 1225 possible two- 
wire pair combinations. The sample sizes ranged from 200 to 600 pair combinations. 
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loss. This was done for the metallic-to-metallic, longitudinal-to- 
metallic, and longitudinal-to-longitudinal crosstalk loss distributions 
by using the formula given in Ref. 4, 


Ni = 20 logio | ; 


where C,, is the capacitance unbalance, 1e., Car,ar,, Cr,ar,, OF C1, 1,, @ 18 
the radian frequency in Hertz, and Zp is the characteristic impedance 
of 22-gauge pulp. The inductive contribution is neglected. The dis- 
tributions are shown in Fig. 3. The rms crosstalk loss corresponds to 
that loss which would result in the average crosstalk power in watts. 
Consequently, crosstalk power transferred between two circuits with 
crosstalk loss equal to the rms value would be the average crosstalk 
power. The rms values are 105.2-dB loss between metallic circuits and 
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Fig. 3—Computed 1000-Hz near-end crosstalk based on the interwire capacitances 
of 416 pair combinations of 1319 ft of 22-gauge pulp cable. 
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only 83.9-dB loss between a longitudinal and a metallic circuit. Hence, 
the new coefficient Cz,.1, predicts that longitudinal excitation of a 
wire pair in the cable measured causes 21.3 dB (105.2 — 83.9) more 
average crosstalk power in a metallic circuit than metallic excitation. 
Similarly, the rms loss between two longitudinal circuits, predicted by 
Cr,1,, 18 49.2 dB, which is 56 dB (105.2—49.2) less than the metallic- 
to-metallic loss. These results are compared in Section 3.6 to crosstalk 
loss measurements made on another cable. In Section 3.5, these re- 
sults are compared to results obtained from computer solution of the 
transmission line equations (9), i.e., by simulation of the cable. 


3.5 Numerical solutions 


A second, more difficult method of calculating the crosstalk between 
the various modes is direct solution of the transmission line equations 
on a computer with an appropriate set of boundary conditions. 

A computer program has been written to solve these equations that 
simulates a cable of the same length and identical characteristics of 
the Western Electric cable used to obtain the coupling coefficients. 
Two conditions of interest were simulated on the computer. First, 
metallic excitation by a balanced 1000-Hz signal generator in series 
with a 600-ohm resistance was applied to a pair, denoted the disturbing 
pair, and all other pairs were terminated metallically with 600-ohm 
resistors from tip to ring. Second, the same conditions applied except 
one wire of the disturbing pair was grounded. This resulted in a deg- 
radation of 14.5 dB in the rms value of the balanced near-end cross- 
talk loss distribution. 

To compare the degradation in crosstalk obtained by the two 
methods, i.e., coupling coefficients versus numerical solutions, it is 
necessary to note that grounding a wire connected to a signal generator 
produces a longitudinal voltage that is one-half the value of the applied 
metallic voltage. This follows directly from the definitions of longi- 
tudinal and metallic voltages in terms of the voltage from each wire 
pair conductor to ground, eq. (6) with V2. = 0. Hence, an approximate 
6-dB adjustment must be made when using the longitudinal-to-metallic 
coupling coefficients Cz,1, and Cz,u,, which predict a 15.3-dB deg- 
radation in the near-end rms crosstalk loss at 1000 Hz due to ground- 
ing, as compared to 14.5 dB predicted by the numerical computation. 
This good agreement suggests that the new capacitive coupling coeffi- 
cients do provide a simple but reliable method of predicting the 
degradation in crosstalk performance for a particular cable when its 
terminations are unbalanced. 
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3.6 Comparison with measured data 


In 1962, measurements of the degradation of near-end crosstalk loss 
caused by grounding one conductor of the disturbing or disturbed pairs 
on randomly selected pair combinations in a 7500-ft length of a 22- 
gauge multipair pulp-insulated exchange grade trunk cable were made 
by Bell Laboratories. This was a real working cable between Oceanside 
and Vista, California. All load coils in the section under test here were 
first removed and the cable ends spliced. The results of these measure- 
ments are shown in Fig. 4. The results reveal, at 940 Hz,* a degradation 
of about 19.4 dB in the near-end rms crosstalk loss when a ground 
was applied to one wire of either the disturbing or the disturbed pairs. 

Since grounding one conductor causes a longitudinal voltage excita- 
tion that is one-half the metallic voltage, a 6-dB numerical adjust- 
ment was made on the measurements to predict that the rms cross- 
talk loss between a longitudinal and a metallic circuit is 25.4 dB worse 
than the rms loss between metallic circuits. This is compared to 21.3 dB 
obtained using the coupling coefficients for the cable discussed in the 
previous section. This 4.1-dB difference may be due to the fact that 
the two cables were not the same, each having different value parame- 
ters characterizing them as well as different lengths. 

When one conductor of both disturbing and disturbed pairs were 
grounded, the measured rms crosstalk loss was 32.6 dB, as shown in 
Fig. 4. This is the loss between the two longitudinal circuits' and, as 
can be seen, it is 61.5 dB less than this rms loss between the metallic 
circuits. This measured difference compares favorably to the calculated 
difference of 56 dB as shown in Fig. 3. The 5.5-dB difference may be 
due to cable differences. In conclusion, direct crosstalk measurements 
on another cable substantially support the analytical method for 
calculating crosstalk using the coupling coefficients or computer 
simulations. 


3.7 Metallic-to-longitudinal conversion because of wire pair imbalance 


So far, we have analyzed the effect of direct longitudinal excitation 
of wire pairs on crosstalk between pairs. This excitation results when 
longitudinally unbalanced terminations are used. However, now we 
discuss how longitudinal excitation can also result because of “pair 
longitudinal imbalance,”’ which is defined as any lack of symmetry 


*940 Hz is close enough to 1000 Hz to permit direct comparison with calculated 
results. 

t Referring to the definitions of longitudinal and metallic voltages, it is simple to 
show that the coupling loss for the longitudinal-to-longitudinal mode is the same as 
for both pairs grounded. 
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in the wire pairs with respect to ground or with respect to each other. 
Such asymmetry can cause part of the metallic signal to be converted 
to a longitudinal excitation even when there is perfect longitudinal 
balance at the terminations. 

To understand the causes of wire pair longitudinal imbalance, 
again refer to eq. (9). Perfect pair balance is the condition that exists 
whenever a metallic signal does not excite the longitudinal modes in 
either the disturbing or disturbed wire pair. This requirement can be 
met if and only if 

G21 = G23 = Gar = Qazs = O. 


These last four conditions are satisfied if 


Ci, = C2, (14a) 
Cig = Cay (14b) 
Cu = Cos (14¢) 
Cis = Cas. (14d) 


Equations (14a) and (14b) are necessary since, for example, if C1, 
were not equal to C2,, there would be a lack of longitudinal symmetry 
in wire pair one even if the terminations were all perfectly balanced. 
Kquations (14c) and (14d) imply that equal and opposite currents are 
coupled (metallic-to-metallic crosstalk) from each of the wires in the 
disturbing pair to the disturbed pair preserving the pair symmetry. 

In other words, if the conditions of eqs. (14) are met and all the 
terminations are balanced, then all the currents are strictly confined 
to the metallic circuits. This is not to say that crosstalk cannot occur. 
It means that only one of the three kinds of coupling can occur, i.e., 
from metallic circuit to metallic circuit. In fact, the crosstalk will then 
be proportional to Campbell’s capacitance unbalance expression which 
simplifies to 


Cuym, = Cis + Cos — Crs — Co3 = 2(Cis — Cu). (15) 


Cable data reveal that the capacitances to ground for wire pairs are 
nearly equal, their differences on the average being less than 2 percent 
of their magnitude. The percent differences in the interwire capaci- 
tances are larger (e.g., 10 percent), but they are much smaller than 
the capacitances to ground. This suggests that metallic-to-longitudinal 
conversion of signals due to the cable characteristics alone is small. 
Computer simulation of wire pairs, using eq. (9) and assuming balanced 
terminations, supports this suggestion. To put it another way, the 
high quality of manufactured multipair cable used in the Bell System 
ensures excellent pair longitudinal balance. The small imbalance in 
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Fig. 4—940 cycles near-end crosstalk coupling loss measured on 60 pair combina- 
tions within a 51-pair unit of a DSAC-202 cable (cable length 7500 ft). All grounds 
are applied at measure end. 


the pairs is rarely a significant factor in the contribution to longi- 
tudinal voltages that degrade crosstalk. Substantial conversion to 
longitudinal modes does occur when there is imbalance at the termina- 
tions, as revealed by our analysis and, consequently, it is necessary to 
place limits on permissible terminal longitudinal imbalance. 


IV. LONGITUDINAL BALANCE REQUIREMENTS 


Crosstalk energy can reach the metallic mode in the disturbed circuit, 
circuit two, from an applied metallic signal in the disturbing pair, 
circuit one, in three different ways. 


(t) Direct coupling from a metallic signal in circuit one to a 
metallic disturbance in circuit two. 

(zt) Conversion of the metallic signal of circuit one to a longi- 
tudinal signal in circuit one because of an unbalanced 
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termination* in that circuit, then coupling of the longitudinal 
signal in circuit one to a metallic disturbance in circuit two. 

(222) Conversion of the metallic signal in circuit one to a longi- 
tudinal signal in circuit one because of an unbalanced termina- 
tion, coupling of the longitudinal signal in circuit one to a 
longitudinal signal in circuit two, and, finally, conversion of 
the longitudinal signal of circuit two back to a metallic dis- 
turbance in circuit two due to an unbalanced termination on 
circuit two. 


The crosstalk described in (7) above is independent of the imbalance 
at the terminations. It is the result of the capacitance unbalance C'y,m, 
between the individual wire pairs and there is little more that can 
practically be done to circuits to reduce it. The important thing ts to 
make sure that any equipment that 1s connected at the terminations of the 
cable does not degrade the low levels of crosstalk that currently exist by 
introducing longitudinal excitations. 


4.1 Longitudinal balance requirement 


The data on the vulnerability of cable to longitudinal imbalance 
have been obtained by measurements made on two different cables.' 
Measurements on the cable in California, with the 6-dB numerical 
adjustment, revealed that longitudinal signals, on the average, cross- 
talk into adjacent wire pairs with 25.4-dB less coupling loss for that 
cable than do metallic signals. The data on the Western Electric reel 
of cable, used in the newly derived capacitance unbalance formulas, 
showed 21.3-dB less coupling loss for longitudinal signals. 

The definition of longitudinal balance, for application to crosstalk 
performance, is repeated: 


em(f) 
ex(f) |’ 


where e, is the longitudinal voltage produced when a metallic voltage 
€m 1S applied at any frequency f. The measurements made on the 
cable in California establish the more stringent longitudinal balance 
requirement, and it shall therefore be assumed that rms longitudinal- 
to-metallic crosstalk loss is 25 dB less than metallic-to-metallic cross- 


BALm_1(f) = 20 logio 








* Conversion because of imbalance in the cable itself can be neglected, as discussed 
in Section 3.7. 

Tt Subsequent to the beginning of this investigation, measurements on the vulner- 
ability of one other cable to crosstalk because of longitudinal imbalance have been 
made. These measurements do not alter the conclusions reached by using the data 
on the first: two cables only. 
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talk loss. This implies that, if a network which is sending a metallic 
signal has a balance of about 25 dB, the longitudinal signal developed 
because of this imbalance may contribute the same amount of cross- 
talk power to nearby cable pairs as the direct metallic signal applied 
to it. It is not known exactly how the two components of the crosstalk 
produced by metallic and longitudinal signals will add, that is, on a 
voltage basis, a power basis, or somewhere in between. However, a 
longitudinal signal developed because of imbalance is likely to be 
correlated to the metallic signal causing it. Hence, it will be assumed 
that the signals add approximately on a voltage basis. 

What is needed is a balance such that the contribution to crosstalk 
power because of imbalance is small compared to the crosstalk that 
exists when a metallic signal is applied. For illustrative purposes, it is 
assumed that an increase of 1.0 dB is not too noticeable and is thus a 
permissible contribution. In Fig. 5, which shows how two voltages 
expressed in dB are added, it is seen that, in order for the power in a 
signal to be increased by no more than 1.0 dB because of the presence 
of a second signal, the voltage difference must be over 17 dB. Thus, a 
longitudinal balance of approximately 42 dB is required (we will use 
40 dB) to ensure that crosstalk is increased by no more than this 
amount, due to the longitudinal-to-metallic coupling path, type (72), 
described at the beginning of Section IV. 

We now show that the crosstalk resulting from the coupling path 
described as type (772) is less severe and has no bearing in determining 
the balance requirement. To do this requires discussing a second 
measure of balance, that known as longitudinal-to-metallic balance. 


SUM OF TWO VOLTAGES 
ABOVE LARGER IN DECIBELS 


0 2 4 6 8 10 12 14 16 18 20 
DIFFERENCE BETWEEN TWO VOLTAGES IN DECIBELS 


Fig. 5—Sum of two voltages expressed in decibels. 
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4.2 Longitudinal-to-metallic balance 


In addition to the possibility of severely degrading the crosstalk 
levels already occurring in the telephone network, there is a second 
important reason why high longitudinal balance at the terminations 
must be maintained. This is to reduce the metallic noise resulting from 
power line induction. Here, the disturbing signals are longitudinal in 
nature and, to assure good performance of the user’s circuit, the con- 
version loss from a longitudinal noise signal on his circuit to a metallic 
signal on his circuit must be large. A measure of this conversion loss, 
defined as longitudinal-to-metallic balance, is: 


ex(f) 
em(f) 


where ¢, is the applied longitudinal voltage source and e€,, is the result- 
ing metallic signal. The balance subscript ‘‘L — m’’? means a conver- 
sion from an applied longitudinal signal to a metallic signal. It is im- 
portant to note that the two measures of balance, BALM, and 
BAL,z_M, are not equal, i.e., reciprocity does not necessarily apply. 
Moreover, they are not necessarily correlated. 

Generally, to assure good performance, the minimum balance 
BAL z_m of a termination is well in excess of 40 dB over the voiceband. 
Consequently, we use 40 dB as a lower bound on the BALy_m of 
terminations, keeping in mind that this in no way establishes 40 dB 
as the necessary performance minimum. Discussion of that topic is 
outside the scope of this paper. 

Using the assumed minimum bound on BALy_m = 40 dB of termi- 
nations on the disturbed circuit, it is now shown that the crosstalk 
coupling loss path (777) is at least 20 dB less than the metallic-to- 
metallic path and, hence, is not a factor. It is also assumed that the 
balance BALm_z of the terminations on the disturbing circuit is 40 dB, 
determined from the crosstalk requirement because of the coupling 
path of type (zz). Now, since there is approximately 60 dB less cross- 
talk loss between two longitudinal circuits than between two metallic 
circuits, the difference in the losses between type (¢) and type (772) 
paths is BALm_y (disturbing circuit) —60 dB + BALy_™ (disturbed 
circuit), or at least 20 dB. 





BALyi_m(f) = 20 log io 


) 








V. REQUIREMENTS ON LONGITUDINALLY APPLIED SIGNALS 


So far, we have considered how longitudinally unbalanced termina- 
tions can cause increased crosstalk, because they longitudinally excite 
a wire pair. We have recognized that it is this longitudinal excitation 
that is the fundamental cause of the increased crosstalk, and we have 
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recommended a metallic-to-longitudinal balance, BALy_—xz, limit for 
termination of 40-dB minimum to restrict the amount of longitudinal 
excitation produced. However, it is also possible to longitudinally ex- 
cite wire pairs directly from a voltage source connected between the 
tip and ring of a wire pair and the cable sheath, or ground. Such direct 
excitation must also be limited, because it too causes crosstalk in a 
disturbed metallic circuit in two ways: 


(¢) Direct coupling from the longitudinal mode in circuit one to the 
metallic mode in circuit two. 

(iz) Coupling of the longitudinal mode in circuit one to the longi- 
tudinal mode in circuit two and conversion of energy in the 
longitudinal mode of circuit two to the metallic mode in circuit 
two because of an unbalanced termination in circuit two. 


Since the effect of directly applying longitudinal signals is the same 
as longitudinal signals arising from metallic-to-longitudinal imbalance, 
and since the rms crosstalk loss for this type of signal is on the average 
25 dB less than metallic signals, longitudinal voltage limits should be 
40 dB more restrictive than metallic voltage limits. Figure 1 shows 
the restriction on longitudinally applied voltages as a function of 
frequency. It is based on the restrictions already placed on metallic 
voltages determined by a previous study at Bell Laboratories* and 
the 40-dB restriction determined here. 


VI. SUMMARY 
The following has been accomplished in this paper: 


(7) Three new capacitive coupling coefficients have been derived 
that provide a simple but reliable method of predicting the 
degradation in crosstalk performance for a particular cable 
when its terminations are unbalanced. 

(27) It has been demonstrated that a metallic-to-longitudinal 
balance requirement of 40 dB or more for any terminations 
connected into network should not noticeably increase the low 
levels of crosstalk that are already present. 

(227) A requirement has been established on longitudinally applied 
signals that if met should not degrade crosstalk performance. 
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APPENDIX A 
In eq. (4), let T(x) be an 8 X 1 column matrix where 


_ | Va) 
r= li | — 
and let A be an 8 X 8 constant matrix (not a function of 2), 
O Z 
A= E ol: (17) 
Then 
dT(x) 
le AT. (18) 


The solution to this matrix differential equation’ is known to be 
T(x) = exp (— Az) T(0). 


Since the parameters that characterize the line are independent of z, 
it is readily seen by solving eq. (4) that 


VFE| _ _[O Z 1 VNE (19) 
CFE} “Pi ly o] |Lene|’ 
where the far-end voltages and currents where x = / are | 
Vid) I,(1) 
V2(1) I,(2) 
ead 0) ea AL 
Vi) I4(1) 


and the near-end voltages and currents are 


Vi(0) 1,(0) 
_ |¥2(0) _ | 12(0) 
VNE = V3(0) |? CNE = T3(0) . 

V4(0) 1,(0) 


With the eight equations given in (19) and a knowledge of the 
boundary relations at the terminations, we can characterize the model 
of the system at each point in space (x) by a vector pair of voltages 
V(x) and pair currents I(<). 

The matrix exponential, 


wo {-[2 ZH. 
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may be evaluated in closed form. With a closed-form representation, 
the voltages and currents can be expressed in closed form and more 
complex structures such as spliced cable systems can be simulated. 
The way to explicitly determine the matrix exponential is to use the 


fact that 
O Z 
ie “\9 0 


satisfies its own characteristic equation (Cayley-Hamilton Theorem). 
Then 


7 
expQ = 2) a:Q*. 
k=0 


Replacing Q by a diagonal matrix consisting of the eight eigenvalues 
of Q enables us to solve for a; However, for fairly short unspliced 
cable systems we may use the first few terms in a power series, i.e., 


exp[-QI] =1-1Q+5Q—-, 


where I is the identity matrix, to accurately approximate the matrix 
exponential as was done for the cable in the numerical solutions section. 


A.1 Boundary conditions 


For any two wire pairs within a cable, four sets of current-voltage 
relationships exist at the wire terminations. Referring back to Fig. 2, 
we define the near end to be the subscriber side of the loop with its 
termination where xz = 0, and the far end, where x = l, to be the 
other termination, possibly a central office. The disturbing pair will 
always be designated wire pair 1-2 with a generator of some kind at the 
near end, and the disturbed pair will be designated 3-4. Suppose the 
generator is two voltage sources each grounded at one end and in 
series with an impedance and the remaining terminations consist each 
of two complex impedances to ground shown in Fig. 2. Then we have 
the eight relations at the boundaries 


Vi(0) = VS(1) — ZNE(1, 1)14(0) (20a) 
V2(0) = VS(2) — ZNE(2, 2)12(0) (20b) 
V;(0) = —ZNE(3, 3)1;(0) (20¢) 
V.(0) = —ZNE(4, 4)1;(0) (20d) 
Vill) = ZFE(A, 1 (1) (20e) 
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V.(l) = ZFE(2, 2)I2(2) (208) 
V3(l) = ZFE(3, 3)13(1) (20g) 
Vil) = ZFE(4, DAD, (20h) 


which may be written in matrix form as 


VNE = —ZNECNE + VS (21a) 
VFE = ZFE CFE, (21b) 
where 
ZN&K(I, 1) 0 0 0 
_ 0 ZNE(2, 2) 0 0 
enti 0 0 ZNE(3, 3) 0 
0 0 0 ZNE(4, 4) 
ZFE(1,1) 0 0 0 
_ 0 ZFE(2, 2) 0 0 
as 0 0 ZF E(3, 3) 0 
0 0 0 ZF E(4, 4) 
and* 


VST = [VS(1), VS(2), 0, 07. 


We can solve the 16 equations (19) and (21) and determine VNE and 
CNE. Now we have the model completely characterized by the vector 
pair of voltages V(x) and I(x) via the equation 


V O Z VNE 
F |= 2 Ly oF tlene]| 

It should be pointed out that the terminations are not always simple 
impedances to ground. For instance, for a second type of termination, 
where an ordinary telephone set is connected to a wire pair, tip, and 
ring, there is no direct conducting path to ground. If the impedance to 
ground from the tip and ring is assumed to be infinite, then we cannot 
write a simple impedance matrix relating the current to voltage as in 
eqs. (21a) or (21b). As a result, there is a rather tedious but straight- 
forward rearrangement of eqs. (19) and (21). Finally, a third type of 
termination could be a central office that will also require modification 
of the impedance matrix. All three of these types of terminations have 
been simulated in computer programs. 


* The superscript 7 means transpose. 
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APPENDIX B 


The elements of the admittance matrix in eq. (13) are given below. 


Qu = {Cig + Cog + 4Ci2 + Cis + Cig + Coz + Cos} 

Qi2 = Qa = 2{(Cig — Cag) + (Cis + Cis — Cos — C22)} 
Qis = qs1 = —(Ciz3 — Cig — Co3 + C24) 

Qua = Qasr = —2(C3 + Crs — Cos — Cos) 

Qa2 = 4{ (Cig + Cag) + (Cis + Cig + C23 + C21)} 

Qos = G32 = —2(C13 — Cra + Co3 — Cor) 

Qos = Gao = —4(Cy3 + Cig + Coz + Cr1) 

Q33 {C39 + Cag + 4034 + Cis + Cig + Cos + Cog} 
Qsa = Gas = 2{(C3g — Cag) + (Cis — Cis + Cos — Co4)} 
qa4 = 4{ (C3, 2 C'4q) + (C13 + Cy, + Cos + C'o4) }. 
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Wideband Amplifier Design Using Major 
Multiloop Feedback Techniques 


By T. J. APRILLE, JR. 
(Manuscript received February 11, 1974) 


Multiloop feedback has heretofore been ignored as a means of obtaining 
shaped gain amplifiers. In this paper, a theoretical basis is developed for 
using dual major loop feedback amplifiers to obtain shaped power gain 
with input and output reflection coefficient constraints. From the theoretical 
results, practical design procedures can easily be developed and one such 
procedure 1s discussed. 

The aim of this study was an alternative to the “‘brute force” termina- 
tion technique of realizing input and output impedance matches. The 
development is otherwise unique in that it uses no hybrid transformers for 
beta circuit coupling or for realization of the reflection coefficient 
constraints. 


Il. INTRODUCTION 


Wideband feedback amplifier design has heretofore mainly been 
accomplished by the use of single major loop feedback techniques.’ 
Major loop feedback implies that the current or voltage on the input 
to the basic amplifying element is manipulated by the current or volt- 
age that appears on the output of the basic amplifying element. The 
design concept follows the classical feedback design procedure of 
assuming a unilateral forward amplifying element of voltage gain u and 
a feedback path with voltage gain 8. Existing multiloop feedback 
techniques have been primarily concerned with stability considerations 
of “tandem” * and minor multiloop?? feedback arrangements. 

In many applications, input and output impedance matching of the 
amplifier is necessary. The communications amplifier is one such 
example, since it requires very low levels of signal interference due to 
input or output impedance mismatch. The classical single-loop feed- 
back techniques offer little help in designing for the impedance match- 
ing constraint. This is due to the fact that the more loop gain in a 
single-loop feedback circuit, the more extreme (zero or infinite) the 
input and output impedance becomes.'~? Two techniques that are 
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used to circumnavigate this problem are “brute force’ terminations 
and bridge couplings. 

The brute force approach obtains the impedance match by placing 
a resistor in series with the input (or output) for a feedback amplifier 
with zero input (or output) impedance. A parallel resistor is used for 
the infinite input or output impedance case. 

The use of a balanced resistor bridge is also useful in obtaining an 
impedance match. This is accomplished by balancing the bridge com- 
ponents with respect to the input (or output) impedance of the feed- 
back amplifier. The use of a resistive bridge is limited, though, due to 
the excessive resistive losses associated with such a bridge. A useful 
four-port device, which exhibits the same qualities as a resistive 
bridge but with much less through loss, is the hybrid transformer.® 
The impedance match with this device is obtained by manipulation 
of the two unused port impedances.’ 

Since the hybrid transformer is similar to a bridge, one of the two 
unused ports can be used for the @ return path. This technique is 
theoretically the best alternative mentioned since a property of such 
a connection is that the impedance match is improved with the 
amount of loop gain.? This technique has been used to advantage on 
several communications amplifiers.®:? 

The limitations of the above alternatives of obtaining an impedance 
match become evident when other design constraints are investigated. 
For example, the noise figure of an amplifier is degraded by any loss 
that exists on the input to the amplifier.!° Thus, the use of brute force 
or hybrid transformer coupling causes an increase in noise figure. On 
the output side, a loss increases the power requirement on the last 
stage of the amplifier. Even if this is no problem, the resultant in- 
crease in the distortion may be. This is due to the fact that second- 
order distortion power increases twice as fast as fundamental power 
and third-order distortion power three times as fast.!! Thus, the losses 
associated with the matching techniques will increase the power re- 
quirement and reduce the linearity of the overall amplifier. 

The use of the hybrid transformer in the 6 path may also cause a 
stability problem. Since the transformer introduces phase shift, due 
both to the physical length and techniques of construction, their use 
is limited at very high frequencies. 

Investigation into alternative methods of design is therefore de- 
sirable. To this end, this paper presents fundamental concepts on the 
techniques of using major multiloop feedback in amplifier design. The 
objective is the design of wideband-frequency-dependent gain ampli- 
fiers with input and output match constraints. The design procedure 
does not use hybrid transformers and attempts to minimize brute 
force termination techniques. 
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In Section II, the basic amplifier element is introduced. The analysis 
that follows is applicable to configurations of active devices that can 
be modeled by this basic amplifier. A circuit form using this basic 
amplifier is then introduced. Matched impedances and gain relation- 
ships are developed for this circuit form in such a way as to make the 
open loop gain characteristics evident. This paves the way for an 
initial design approach that is independent of the loop gain 
characteristics. 

In Section III, a second circuit form of shunt-series feedback using 
the same basic amplifier is introduced. Matched impedances and gain 
relationships are again developed. The derivations in this section 
exactly parallel those of Section II. 

In Section IV is given the results of the two previous sections to 
demonstrate the procedure used to obtain an initial circuit design for 
a practical amplifier configuration. The configuration treated is that 
of a cascade of N common emitter transistor stages. It is shown that 
for N odd, the results of Section II can be used, and for N even, the 
results of Section III apply. One numerical example is supplied for 
each case. Two appendices provide the calculations used to derive the 
results in Sections II and III. 


Il. SHUNT TRANSADMITTANCE: SERIES TRANSIMPEDANCE FEEDBACK 


Each dual-loop feedback amplifier discussed in this paper contains 
three major components: two feedback networks and one amplifying 
element. Each major component is assumed to be made up of any 
number of passive and active elements. Characteristics of importance 
for the amplifying element component are given in Fig. 1; this ab- 
breviated model is designated a basic amplifier. In this figure, z, is 
the input impedance and J, is a current-controlled current source. J; 
is given by the product of a frequency-dependent variable k and the 
current through 2,. 

In Fig. 2, the first multiloop feedback circuit form is given. Series 
feedback voltage source al, sums up the most important characteristic 


INPUT 2x | OUTPUT 


Fig. 1—Basic amplifier. 
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I, = kix | ZL 





Zins Pin 


Fig. 2—First feedback form. 


of one feedback network. Shunt current source —bV, likewise is the 
important characteristic of the second feedback network. Since the 
series feedback voltage source is dependent upon output current, it 
represents a transimpedance feedback. Similarly, the shunt current 
source is dependent on the output voltage yielding a transadmittance 
feedback. Source and load impedances, Z, and Z,, summarize amplifier 
interaction with the driving circuitry and the loading circuitry, 
respectively. 


2.1 Input and output impedance 

Zin and Zour, the input and output impedance, are desired to be 
matched to Z, and Zz, respectively. Thus, Zin and Zout are needed and 
are given by 
22 + ka 


Aha eg (1) 
_ &e + Z, tka 
celal ae 


If the amplifier gain k is large, then Zin and Zour become 


ak a 

Zin = bkZ BZ’ (3) 
ak a 

Lisi = okZ, = bZ. (4) 


For the matched condition, Zin = Zand Zou: = Zz. Using these condi- 
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tions in eqs. (3) and (4) yields 


* a 
Ao Bia og (5) 


Zi =. Zivat = 7k (6) 


Substituting (Zz)* = Zz from eq. (6) into (5) yields 


a 
= SGP)’ 7) 


a* a\* a 
mao (E) HE (8) 
Thus, conjugate matching yields the requirement that the ratio of 


a to b (or more generally ka to kb) must be real. Given this fact, eqs. 
(5) and (6) are identical, i.e., 
a 


Z2Z,= 2,0. = 7 (9) 


Z; 


which implies 


The imaginary part of Zz Z} is therefore constrained by 
Im {Zz} Re {Z,} — Re {Zz} Im {Z,} = 0. (10) 


A necessary condition for an amplifier to be absolutely stable is that 
Zin and Zour be passive.’ This is satisfied when the real parts of Z, 
and Zz, are positive. Thus, the imaginary part of Z, and Zr, have 
the same sign, implying that if the matched load impedance is capaci- 
tive (inductive), then the matched source impedance is capacitive 
(inductive). 

If the reflection coefficient [reflection coefficient p is defined as 
(Z — Zret)/(Z + Zret) 4 at the input is evaluated (assuming 
Z, = a/bZ;), the following is obtained: 


eee 2ab Re (Zr) i 
Pin = pin(1 ae ein ta” F 


In eq. (11) pin, is the input reflection coefficient when k = 0. 
Evaluating the return ratio 7’ (Ref. 2) of the circuit in Fig. 2 with 
respect to the output dependent current source gives 
—2ab Re (Zz) k 
202, +a — 


Return difference F (Ref. 2) is defined as 1 — 7; thus, eq. (11) can be 
rewritten as 


(11) 


T= (12) 


1 
Pin = Piny F . (13) 
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Therefore, the return loss (20 log 1/|p|) is improved with increasing 
return difference F. Since the output reflection coefficient is given by 


js 
Pout = Pout, FP ’ (14) 


it also realizes the same improvement with increased return difference. 


2.2 Gain equations 


The transducer gain for the circuit in Fig. 2 can be calculated when 
it is assumed that Z, = a/bZ{: 


ere es 


2— Feces pseu are 
Sail? = Topy [i = TP 


(15) 


Again, T is the return ratio and is given by eq. (12). 
In eq. (15), 7 is proportional to k. Thus, for large k, |S2:|? goes to 
1/|ab|. Therefore, eq. (15) can be rewritten as 


[7 /? 


| Soil? = [Seo]? pp (16) 
where 
1 
| Soro]? = Tab] (17) 


2.3 Design procedure 


In the derivations given thus far, a definite effort has been made to 
separate the dependence of k. This was done for two reasons: (2) to 
allow an initial design to be effected with k not a variable, and (iz) 
to allow definitive statements to be easily made concerning the effects 
of k. The former can easily be implemented by assuming k = ~. 

In this case of k = ~, eqs. (8), (9), and (17) are relevant. These 
equations are repeated for convenience: , 


a\* a 
(5) 3 . 
4,0, == (9) 
| Se1e|? = aT (17) 


It should be noted that eq. (8) implies that a/b is real, but a and b 
can be complex. |S2i.|? in eq. (17) is the maximum available gain 
since it is obtained with the input and output matched. 

In summary, the design procedure given below could be used when 
the desired gain g and impedance matches are known. 
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(2) 
(22) 


(272) 


(2v) 


Choose an arbitrary starting Z, and Z, such that Z,Z;, is real. 
Substitute a from eq. (9) into eq. (17), yielding 


1 


= Of ea 
g | So120| ZZ \b |? (18) 
Synthesize b such that 
1 
| = ——— (19) 
Z.Z1\9| 


There are no constraints upon the phase of b except those that 
may result from stability considerations. 
Synthesize a such that 

a = bZ,Z}. (20) 


The value of k is now obtained by considering the practical 
active devices used to simulate the ideal amplifying element. 
With & known, the return ratio T [eq. (12) ] can be calculated ; 
this yields the obtainable impedance match and gain devia- 
tion, eqs. (13), (14), and (16). If the design objectives are not 
met, the previous calculations should make the necessary 


changes evident, e.g., lower Z, or a higher value of k. 


Ill. CURRENT TRANSFER SHUNT; VOLTAGE TRANSFER SERIES FEEDBACK 


The last multiloop feedback circuit to be considered is shown in 
Fig. 3. In this case, the series feedback voltage source is dependent 
upon the output voltage and thus represents a voltage transfer feed- 
back. Similarly, the shunt current source is a current transfer feedback. 
The voltage source is given by aV» and the current source by —blo, 


otherwise Figs. 2 and 3 are identical. 


3.1 Input and output impedance 


The input and output impedances, when evaluated, are given by 


= ee akZ zy 

aaa a oa 
_ Set es — kZb 
Lout —_ Te Seepage > (22) 
For large k, eqs. (21) and (22) become 

Zin = 5 Zi, (23) 
Fey (24) 

out —_ a Se 
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I, = kix “| 2. 





Pout « Zo Po 


Fig. 3—Second feedback form. 


The conditions of input and output matches yield 


OP. a »_4(0\" > 
Bi = Zin = § 2a = $ Cou)* = § (2) Zh (25) 
Thus, 
* 
(5) = (é): 26) 
and 
i a 
ay a me 


Since a/b is real and the real part of Z, and Z, are nonnegative, then 
eq. (27) implies that if Z, is capacitive (inductive) then Zz, must be 
inductive (capacitive). 

Input and output reflection coefficients can be evaluated along with 
the return ratio and return difference. The results are shown below 
for Z, and Zz, satisfying eq. (27). 


2kab Re (Z1)_ 


ie bz, + aZz, 2 
F=1-T. (29) 
= 1 ., ete bz. — aZy 
Pin = Ping FE’ Ping — bz, a aZ* (30) 
1 1 
Pout = Pout, F ae PF’ Pouty = 1. (31) 


Thus, as in the case of the first circuit form, the reflection coefficients 
are improved by the return difference. 
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3.2 Gain equation 
When the load and source impedances satisfy eq. (27), the trans- 
ducer gain for the circuit in Fig. 3 is given by 


[Ele 


| Soil[? = [Sate]? 7G Tl)?’ (32) 
1 
2— —__.. 
| S21. | [ab | (33) 


This is the same form as was given in eqs. (16) and (17) ; thus, the same 
statements apply to the above equations concerning improvement 
with feedback. 


3.3 Design procedure 


Initial circuit design can proceed in a manner similar to the first 
case. The term k again is assumed equal to infinity; this yields the 
germane equations summarized below. 


* 
(2) =5 a» 
Zs 
7 . (27) 
1 
| Soteo |? = [ab| (33) 


The four design steps outlined previously apply except as follows. 


(i) Choose Z, and Z; such that Z5/Zxz is real. 
(22) Substitute eq. (27) into (33) so that 


1 


= 2— —_—. 
g = |Ser0| Ca (34) 
Z1 
Synthesize b such that 
1 
0 =e 35 
|| Zs | (35) 
Da? 
(212) Synthesize a such that 
Zs 


(iv) With & known, the return ratio is obtained from eq. (28). 
Equations (80), (31), and (82) then yield the obtainable impedance 
matches and gain deviation, respectively. 
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IV. DESIGN EXAMPLES 


Results obtained in the last two sections will now be applied to a 
basic amplifier consisting of a cascade of N common emitter transistor 
stages. Transistors will be assumed to be used in a frequency range 
well below cutoff. The first case to be treated is for N odd. 


4.1 N odd 


Consider the circuit given in Fig. 4a. In this circuit, the transistor 
will be modeled by the circuit given in Fig. 4b. The circuit given in 
Fig. 4a will now be converted to the form given in Fig. 2. Z, and Zz, 
have their obvious counterparts. z, is given by the impedance from 
base to ground with al» equal to zero; this is obtained when J» = 0, 
which can be obtained by setting ai (first stage a) to zero. From the 





Fig. 4—Design example N odd. 
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transistor model, this is evidently 
22 = tre t Zu. (37) 


I, is given by the collector current on the Vth transistor when the base 
current on the first is unity. With the cascade of N transistors, this 
yields 


ay 





I, = Bi, 6B dy? 6, = eae k = fi, eS BN (38) 


1 — a,’ 


The term a is given as the value of open circuit input voltage (Z, and 
Zr removed) that exists when J» equals unity. This is given by 


Letre 


a= fy+ 7 ei, N>1 (39) 
a=Zetfe, N=1. (40) 

For |G, ---, Bv|>>|Ze + rel, eq. (39) can be approximated by 
a= Zz, N>1. (41) 


The last remaining parameter b can be obtained by evaluating the 
Yi. parameter of the network Y, yielding 


b= -Wd=-(F)-F (42) 


Loading effects of the Y networks, i.e., y11, y22, can be ignored if they 
are sufficiently small. 

As a numerical example, the value of k, Zr, and Zz are calculated 
to yield an input and output reflection coefficient of 0.18 (return loss 
of 15 dB) and a gain to within a factor of 1.26 (1 dB) of f?/400, f in 
MHz, in the band from 80 to 140 MHz. The remaining parameters of 
the transistors are 7, = 1, andr, = 0.173. 

Following the four-step design procedure yields 


(i) Let Z, = Zr = 20, ZZ, = 400. 


ae re | 
(1) 9 = 309 = 40/0. 
ae ie! 

poe ian 2 








Z;| = f, fin MHz. 

If Zs is chosen as an inductor, then 
Z;| = 2rfL = f; thus, 

L = 1/27 wH, and 

Be oi a alte, 

— OefLj ff) 





b 
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1 _ 400 





(27) a = Z,Z,b = 400b = 400 Z, > nL fi 
400 : 
a= 5! thus, a can be realized when 
Zn(a) is a capacitor of value 
1 


C = 2 (400) pF. 
(iv) Using eq. (12) and a & Zz, the following is obtained: 


Pus —2ab Re (Xx), _ (-2 50 7.20)a / 


ZadZy +a fj fj 
400\ 1 400 
1.173 + 0) =. 20 4 =F | 
( ihe ae Ss) fj 
T= —16000k — {—0.46e-*"7"k at 80 MHz 
~ 8000 + f7423.4 =| —0.27e-**k at 140 MHz. 


For k = 20, |1— T| = |F| at 140 MHz (the worst case point) is 
given by |1 + 5.4e-7°| = [1.75 — 75.85| = 5.63. This reduces the 
reflection coefficient by 1/5.03 = 0.18. Thus, the input and output 
reflection coefficient specification is initially satisfied. 

The gain deviation at 140 MHz is calculated from eq. (15) and is 





This implies a gain deviation from nominal of 0.37 dB, and initially 
satisfies the design requirements. 





Eh (2m) (400) pf 


Fig. 5—Numerical example N odd. 
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The resulting circuit for N odd is given in Fig. 5. The actual trans- 
ducer gain, input impedance, and output impedance values for this 
circuit were obtained by a computer-aided design program.!® A sum- 
mary of the results is given in Fig. 6. It can be seen from this figure 
that this procedure yields a practical first iteration in the design 
procedure. 


TRANSDUCER GAIN IN DECIBELS 


22 


20 


80 


ig COMPUTED INPUT RETURN LOSS 


a“ 
7 


ig COMERTED OUTPUT RETURN LOSS 


DESIRED GAIN —~ 


“N 
“~ COMPUTED GAIN 


DESIRED RETURN LOSS 





100 120 140 
FREQUENCY IN MEGAHERTZ 


Fig. 6—N odd results. 
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4.2 N even 


For an even number of transistors, the circuit in Fig. 7 is used. The 
transistor model of Fig. 4b is again used. The evaluation of the param- 
eters of Fig. 3 follows on the same basis as in the N odd case. Using 
Fig. 7, the calculations are summarized below. 


tere tre + Zu, (43) 
I, = Bi, +29, Byte; k= — Pi; C837 BN, (44) 
ie Zr seen (200Z z)/(Zz + 200) 
200 + Ze Z 182, meee ON 
~ 500° B.> 1, t= 2, pay \Zn| <200 (45) 
2 2 


b=7 45 ™zZ," |Zr| > 2. (46) 
A numerical example is given to show the initial design steps for 
obtaining a maximum input and output reflection coefficient of 0.18 
and a gain to within 1.26 of f? (f in MHz), from 80 to 140 MHz. The 
transistor parameters are again 7, = 1 ohm and r, = 0.173. 
The four design steps yield 





(1) Let Z, = Zp = 20, 25 = 1, 
Zi 
(i) g= f= 
oe ile 
1 y 
b| =2=—a 
ame mae 
|Zr| = 2f 
Let Zp = j2f = j2nfL, Lin pH 
bea 
Tv 
eet a, is aan a A 
(at) a=b= jie fL ~ 300 ~ FF 
200 
Thus, Zz = —>: 
ae: 
This implies that Zz is a capacitor of value 


1 
ee 2m (200) Be: 
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Fig. 7—Design example N even. 


(zv) Using eq. (32), 


p = 2abRe (Zz), _ A/G A/G N20 
bzz + aZy, (1/jf) (ry + re + 200/jf) + 20/7 f 
T = 40 a 0.023e-*"k at 80 MHz 
~ 200 + 21.1737f° ~—-| 0.015e~88°k at 140 MHz. 


For k = —400, |1—T| =|F| at 140 MHz is given by |1+6e—®*"| =6.1. 

The reflection coefficient is reduced by a factor of 1/6.1 (15.7 dB). 
Gain deviation can be calculated and is equal to 0.95 (0.22 dB) ; thus, 
the initial specifications are satisfied. 

Figure 8 gives the resulting circuit. The results of the computer 
analysis of this circuit are given in Fig. 9. Again the data show that 
the approach yields good results. 

It can be seen in Fig. 9 that the difference in gain is greater than the 
computed 0.22 dB. This is due to the fact that a was taken as Zz/200, 
rather than the term given in eq. (45). A more accurate evaluation 
(denoted by the hatted variables) of a is given as 





e= Ze 4. te + (200Z2)/(Zz + 200) _ Ze Le 
~ 200 + Zz Z1B2 “200 ' Z1B2 

a mw LE Ze _ 4E es 

aX 50 + (20) (80) ~ 200 (1 + 0.125) = 1.1254, 


where a was the numerical value previously obtained. Using eq. (33) 
yields 
| ape ena ee 
[ab| -|1.125ab| «1.125 jab] 1.125 |"! * 
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0.833 iey 0.988 ie2 





Fig. 8—Numerical example N even. 


Again the unhatted quantities were the ones used in the five design 
steps. The factor of 1.125 accounts for an additional 0.51 dB of the 
gain difference. 

In this example, the gain difference can easily be reduced by in- 
creasing the 8 of the second-stage transistor. This was not done since 
it was desired to keep the 6162 product at 400. Since 6: is 80, this 
forces 8; to be 5; any high value of 82 results in unrealistic values of §;. 
Nonetheless it is evident that a high 2, ---, 8, product is needed for 
an even number of cascade stages. 


V. CONCLUSIONS 


In this paper, the basic characteristics of two forms of major multi- 
loop feedback have been investigated. The design characteristics 
treated have been input and output impedance and frequency-de- 
pendent power gain. It has been shown that, with sufficient open loop 
gain, the equations that describe the gain and impedance quantities 
are very simple in nature. An initial circuit-design iteration can easily 
be performed since many complicating variables are eliminated. 

This initial circuit-design concept would be extremely useful in 
a computer circuit analysis-optimization program. Well known is the 
major practical limitation of optimization programs: the obtaining of 
a convergent starting point. For dual-loop amplifiers, this paper offers 
the designer a method of easily finding a good starting point. 

Although not reported here, several frequency-shaped amplifiers 
were actually built using multiloop feedback. The excellent perform- 
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_—~ COMPUTED INPUT RETURN LOSS 
a 


DESIRED GAIN ~_ 
~ 


COMPUTED GAIN 


RETURN LOSS IN DECIBELS 


_7 COMPUTED OUTPUT RETURN LOSS 
a 


TRANSDUCER GAIN IN DECIBELS 


DESIRED RETURN LOSS 





80 100 120 140 
FREQUENCY IN MEGAHERTZ 


Fig. 9—N even results. 


ance of these amplifiers, with respect to input and output matching 
and gain shaping, has precipitated the work reported in this paper. We 
anticipate that future papers will discuss more complicated active 
devices, feedback network loading effects, and feedback network 
synthesis. 
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APPENDIX A 
Shunt Transadmittance: Series Transimpedance Feedback Calculations 


A.1 Calculation of Zin 


Vin = aki, 
t= -—*S:?YS 
ez 
tz(22 + ak) = Vin; 
Vo = ca kizZ y. 


Thus, 
Ts: = Le = bVo = 4s + kbt.Z1; 


= -_ UA+kbsZ1) yy . 
(1 + kbZ 1)tx = “2,4 ak V ins 
Z. = Vie: 22 + ka i 
7 Tin 1 + kbZ1 


Aa 
5 
| 





A.2 Calculation of Zou: 
: Z aly 
z= obV : — : 
‘ "zz +Z, @2+2Z,’ 


St ed = Zs _ alo F 
kis = Ty = k (bVo zt ) K(--): 
Io(@z + Z, + ka) a kbZ Lo; 


_ Vo_ Get Z,+ka) 
ey ie kbZ, 














Zo 


A.3 Input reflection coefficient calculation 
 _4in-Z, 7 _ tthe » _ 4 
Pin Gin t+ Za’ in T+ kbZ 1’ ae Are 


ss ee) / Zz + ka a 

pm = (ee bZ 1 (Ae +e) 
_ @2bZ1r + abkZ1, — a — abkZy 
= 7.bZt + abkZ, +a + abkZy’ 








wf eae OA IO 2 ee 
Pin" z2bZ7 +a + abk[2 Re (Zz) ]’ 
en teh = 
Ping ~ Piny x29 — Z20Z, +a 
Therefore, 


— _ &2bZi — a 2ab Re (Zz) y 
pin = EG X (1 zbZt +a” 
2ab Re (Zr) ky 


ars (1 a zxbZ7, + a 
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A.4 Return ratio calculation 

Assuming the output current source J, is disconnected and replaced 
by a l-ampere current source, 7’ is given by the current that flows 
through the disconnected current source. 








ee Ga Z; _ aly ; _ = 
TON AT Ga: UNO, MES AR 
ee ee ae en 2 Ss ty 
le = bau 4g. ae ee 2.+Zs ) Z, = bZ*? 
T=ki = pp 041 (@)/ (027) — a _ _ ,, abZ1 + abZi, 
: 2z + (a)/(bZ7) b2,0%, +a 
> —2ab Re (Zz) i 
bz.Z;, +a 


A.5 Output reflection coefficient 


— Lout — Zh. Z . Oe Zi he. Zr 
ae out — ) 


a 
Post Zout + Bx’ kbZ. ma 


a (tthe a / tet Z.+ ha, a 
peu kbZ, bZ, kbZ,, bZ* )? 
z,+2Z,+ ak — ak / 7 OT ako a at.) 
kbZ, kbZ 2" , 
a zxZ,+ ZZ 
Re Fr A aka sen) 
Zn4,+Z.L5 


= tf ula pe pk 1. 
Pouty Pou | k=0 ZL. + VAT Se 


ll 


Pout 


Therefore, 
= (1 4 u(Z, + 25) a ee 
PS OE aaa eee 
Sica: wen aD) 
Pout = Pouty 22(a/bZ1) + (a/bZ;) (a/bZ 1) 
= kabZy + kabZz,\~* 
= Pouty (1 1 2.0L". cag ) ’ 


2ab Re (Z = 


A.6 Transducer gain calculation 


Assume a voltage source of value V, is inserted in series with the 
source impedance Z, in Fig. 2. Let Z, = a/bZz. Pas will denote the 
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available power from the source and P,y: the real power delivered to 
the load Zr. 


Wales: Vel 1A LOL: 
4Re(Z,)  4Re(Zz)\al’ 


ViFia zz + ka Z. + ka 
ey a Tes (v. i+ ae )/ (aes ee a) 


(2, + ka)bZ;, 


Pas = 





= ghee ees Rabe e 
V.. = (zr + ka)bZiV; _ (2 +ka)bZ, 1 Vv 
"a+ 2.bZ, + 2kabRe (Zt) aa +2202, 1-T ° 
Vin — aly _ : _ kis _ Lo 
Re ee kok 
Thus, 
k 
Io ake im 
= II" _ [kKPRe (Zi) |p, 
Pout = Iolo Re (Z1) = [22 = ka|? Vg ’ 
pp _ [k|? Re (21) |ze + hal?|b|2|Zc|? _|Vo|?_ 
a |z. + kal? la + 2.bZ7,|? [1 — T |? 
Therefore, 
Sa jt = Pow — Lbl? Re x) |b|*] Z2|2] V.|*4 Re Zs) La 
a Pas ret y AGES dULAGLAL |b 
_ kp Re (Ri)labl4 TP? 
las 22bZ;(?|1— 2? lab] (1 7? 
APPENDIX B 


Current Transfer Shunt: Voltage Transfer Series Feedback Calculations 
B.1 Calculation of Zin 


Vin = aV 9 = Vin -+- akizZ 1 


tz = 
ex ea 
i. = Veen : 
* Be — akZ y’ 
ayer - _ at Vin =— bk) . 
lin = blo + tz = ( bk + 1)t. = 2Z — akZ,’ 


= 


= in _ ee akZ iy 
SE a aac AST 
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B.2 Calculation of Zou: 
7 bloZs _ aVo _ Z dL —_ aVo. 
Te, Uae. Leh ay * 
= kZ.bIo — kaV o. 








Ze2> .* 

I)X(Z, + 22 — kZ.b) = — kaVo; 

Fe Vo _ 4st 22 — kbd. 

out I, —ka 

B.3 Return ratio calculation 

Zs 1 _ ae ; 
Uz Oe arg Ne I, = 1, Yo= ZL; 
est bZ; -aZy : _ a px, 
te gee Ze + Z,’ a= 5 4 


ge ee kb(a/b)Zz, + aZzk _ 2kab Re (Zr) 
a + (a/b)Z, ite + TZ, 


B.4 Input reflection coefficient calculation 


_ Zin = Z3, _ 2, — akZ, _ : . 
pe Fie Zin = - T= He Z, = (a/b)Zz, (a/b) real; 


_ (22 —akZ, aby 22 — akZ1 , @ 
pin = (FS) / (ee HG) 


= bz, = abkZ 1, —= aZr, + abkZ 1 
bzz — abkZ, + aZ, — abkZ},’ 





Sena bz, _ aZt 
Pin ~ bzz + aZt — abk Re (Zz)’ 
bz, — aZr 


Ping = Pin) k=0 — bz, + aZ;,’ 


therefore, 
ore ee 2abk Re (Z1) \~! 
Pin Ping bz2 + aZy ° 


B.5 Output reflection coefficient calculation 
—Four- 2p, 7 _ Zi t%:— kbd 
Pout Vie ae Zr’ out —ka ) 
Z, = (a/b)Zj, (a/b) real; 
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Z,+ 22 — kZ.b ‘ Z, t+ 22 —kZb 
Pout (Ste ae _ gy) [( Ate? + 21) 
_ 2,422 —kZ.b + kaZy _ Ze t 22 
~ Ze+24,—kZ,.b —kaZ, Z, +22 — ka2 Re (Zr)’ 


Pouto = Pout) 0 = 1, 


( _ 2ka Re (Z1) y" 7 ( _ Qkab Re (Z1) 5 
POM GID e oe 2s 7 <Dexte aZ, + bz; 


Pout = 


B.6 Transducer gain calculation 

Assume a voltage source of value V, is inserted in series with the 
source impedance Z, in Fig. 3. Let Z, = (a/b)Zj, a/b real. Pag will 
denote the available power from the source and Pou the real power 
delivered to the load Zr. 


2 2 
Brody A 1 


4 Re (Z.) — FRe (Zz) \al’ 


oe Vale 22— akZ;, / Q px , 22 — akZy 
Vin = gee = (V, 1 — bk ) ($2 + 1 — bk ) 


= V (bz, — abkZz) 
~ aZi, — abkZ, + bz, — abkZ yz’ 


: bz, — abkZy bee — abkZn 1 
eS Ve Oo + ba, — 2abk Re (Z1) Vs aZy, + bz, x 1-—T’ 
Vin —aVo_. _ Riz — In | Vin talodr, 
Ze RS eo, ge Ze : 
Ih = Vin (% - az) = - A+ Vin; 
—s 2 Ele oe 
Pec = Iolo Re (Zz) = lz, _ ae (Zr), 
P= |k|? Re (Zz) |bzz — abkZ1|? LAL 
out lze2 — akZz|? |aZ% + bz.(2 [1 — TI? aaaE : 
p _ lel?|b|? Re (Z1)|V." 
out" |aZ, + bz2(2{1 — TI?’ 
therefore, 
Sul? = Ze = |k|?|b|? Re (Zz) |a|4 Re (Zz) | V.|? 
Pas |aZr + bz,|?|1 — T|?|V.|?|d| 
ape eee ae 
[ab| |1 — TI? 
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Variance of Load Measurements in 
Markovian Service Systems 
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The load carried by a queuing system under equilibrium conditions 1s 
the average amount of server usage per unit of time. In telephony, this 
parameter is often evaluated by recording the number of busy servers at 
regular tume intervals; these readings are then cumulated and their sum, 
after division by the number of observations, is an unbiased estimate of 
the carried load. The purpose of this paper ts to derive exact formulas 
for the computation of the variance of this measurement in systems with 
arbitrary input and departure rates. The results obtained here thus apply 
to a wide class of teletraffic models which includes, in particular, the delay- 
and-loss systems with finite- or infinite-source inputs, exponential service 
temes, and arbitrary defection rates from the queue. Problems related to 
computations are also considered, special attention being paid to the 
reduction of both computer time and storage when the number of states is 
large. 


I. INTRODUCTION 


Analysis of the stochastic behavior of traffic measurements is of 
considerable practical relevance, as it provides means for appraising 
field data as well as guidelines for selecting performance standards. 
Load measurements play a central role in this effort, and determina- 
tion of their accuracy is therefore of particular interest. The present 
investigation yields an answer to this problem for a broad class of 
teletraffic models. 

Whenever statistical equilibrium prevails (and it is assumed to 
throughout this paper), the load carried by a service system is the 
average amount of server usage per unit of time or, equivalently, the 
average number of busy servers at an arbitrary instant. In telephony, 
an estimate of this parameter is often obtained by “switch-counting.”’ ! 
This statistic, which is determined by recording the number of busy 
servers at regular intervals and then by taking the arithmetic mean 
of these discrete observations, is an unbiased estimate of the carried 
load. 
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The variance of this measurement, called hereafter the switch-count 
load to distinguish it from the estimate obtained by continuous observa- 
tion, was first determined approximately by Palm? and Hayward! in 
the case of an infinite server group with Poisson input and exponential 
holding times. This result was later extended by Bene§,’? who obtained 
the exact variance of the switch-count load for groups of finite sizes 
without waiting positions (loss systems). A further generalization to 
loss systems with recurrent input and exponential service is due to 
Neal and Kuczura.! Their formal analysis stops, however, with a 
derivation of the Laplace transform of the covariance function of the 
underlying carried-load process. From this point on, they proceed 
numerically, since explicit inversion of the transform appears to be 
difficult in general. 

In this paper we are concerned with derivations of exact formulas 
for the variance of the switch-count load in finite systems with arbitrary 
state-dependent input and departure rates. The results presented here, 
therefore, fill a rather large gap, since they apply to a broad class of 
teletraffic models that includes, in particular, the (finite) delay systems 
with exponential holding-time distributions, arbitrary defection rates 
from the queue (if one is allowed to form) and either Poisson or 
quasi-random input (in the latter case, the traffic is generated by a 
finite number of sources that place demands for service at the same 
constant probability rate when free but that do not submit requests 
while being either served or waiting). 

Let N (t), the state of the system at time t, be defined as the number 
of busy devices at that instant (by device, we mean here either a server 
or a waiting position). Let c and d be, respectively, the number of 
servers and the number of devices. 

Unless stated otherwise, we make the following assumptions: 


(.) When N(t) = n and 0 Sn <d, the probability that a re- 
quest originates during (¢, ! +h), h > 0, is of the form \,A 
+o(h), with \, > 0. 

(71) The requests which are submitted when all the devices are 
occupied are dismissed and, accordingly, Aq is set equal to zero. 

(77) When N(t) = n and 0 < n Sc, the probability that a service 
time terminates during (t, t+ A) is of the form ph + o(h), 
where pp > 0. 

(7v) When N(t) = n > c, the probability that either a service time 
terminates or a waiting request defects from the queue is of 
the form wh + o(h) where up,» > O and» S d. 

(v) When a server becomes free, it is immediately reseized by one 
of the waiting requests if any are present in the system at that 
time. 
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Let N.(t) be the number of busy servers at time ¢ and let c A n be 
the smaller of the two integers c and n. Then 


N(t) if N(@) Sc, 
N.(t) =c a N(t) = 
c if N(t) >, 


and the switch-count load, L,(T), based on n observations (scans) 
made over [0, 7] at times 7, 27, ---, nz, is, by definition, equal to 


nO N.(jr), 
j=l 


where 7 = T/n. 

Let Cov [N.(t:), Ne(tz)] be the covariance between N,(t:) and 
N (tz). Under equilibrium conditions, this covariance depends only on 
| ty —= to | so that 


Cov CN e(t1), N.(te) | = Cov [N.(0), N.({t1 a, t»|)]. 


Hence, the variance of L,(7), cast in a form that will be convenient 
later, is given by the formula (Ref. 3, p. 137): 


Var L,(T) = n= = (n — |k|)Re(kr), (1) 


where 
R.(kr) = Cov [N.(0), Ne(kr) ] 
= Cov [N.(0), N-({k| 7) ]. 


It is clear from (1) that the variance of the switch-count load is 
completely determined by the covariance function R,(-) of the carried- 
load process {N,(t), —2<t<}, and therefore much of what 
follows is concerned with expressing R,(-) in the most convenient form. 

The covariance function can be stated at first in terms of the transi- 
tion probabilities, and the resulting expression can then be reduced 
by taking the structural properties of the process into account. But 
alternate forms can also be obtained by making use of the fact that 
the conditional expectations, E{N.(t)|N(0) = m}, m = 0, 1, ---, d, 
satisfy simple linear differential equations. The covariance formulas 
obtained by these diverse procedures exhibit distinct features that may 
be exploited in the computations. In all cases, however, R.(¢) is ex- 
pressed as a diagonal, positive-definite quadratic form which reveals 
that R.(-) is completely monotonic.® 

Expressions for the transition probabilities, the covariance function, 
and the variance of the switch-count load are derived in Sections IT, 
ITI, and IV, respectively. The variance of load measurements based on 
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continuous observations can be found in Section IV. Extensions of 
the results of Sections II to IV to reversible Markov processes are also 
considered. Questions of a computational nature are dealt with in 
Section V, while Section VI is concerned with some asymptotic 
properties (d large) of the spectrum of the underlying process. 

The formulas presented here have been programmed and used to 
explore the effects of parameter changes on the variance of the switch- 
count load. The result of that investigation will appear in another 
paper. 


Il. TRANSITION PROBABILITIES 


In this section, we express the transition-probability matrix as a 
symmetric product of vectors and matrices. As becomes apparent later, 
this representation makes it possible to write the variance of the switch- 
count load in a way that greatly simplifies its evaluation. 

Let pmn(t) be the probability of a transition from state m to state n 
in time ¢: 


Pmn(t) = Pr LN () = n|N (0) = m], m,n = 0, 1, es d. 
These transition probabilities satisfy the following system of 
differential equations: 
a ee 
dt P mo im, Mipm1 0P mo, 


d 
aim = bhn+1Pm,n+1 + A\n—1P m,n—1 aad (Ra + Ral Din, 1 s n< d, 


d 
= Did = Nd—-1Pm,d—1 — MPdDma- (2) 
dt 
Let 
—Xo No 
My — Ai+e1) AL 
he Me —(Ae+H2) Ae 
Md-1 —(Ag-1+ma-1) Aa-1 
bed —pa 
and 
Pa (t) = [Pmn(t) |, m,n = 0, 1, Sig d, 


be the transition-probability matrix. Capital and lower-case bold-face 
letters are used exclusively to designate matrices and vectors, re- 
spectively. 
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With this notation, the system of differential equations (2) becomes 
4 Pa(t) = Pu(i)-Ay 20, (3) 
so that, fork = 1, 2, ---, 


© Pili) = = P.(t): Ae t 2 0. (4) 


= 

It follows from our assumptions that if the system is in state m at 
time zero [NV (0) = m], then limao Dam(t) = 1 and limo pmn(t) = 0 
for n # m. Hence, with Iz the identity matrix of order d+ 1, the 
initial conditions take the following form: 


P.(0) => lim P,(t) = Ta, 
t40 
and by (8) and (4) we therefore have 
li op (i) = A (5) 
le 


The initial conditions state that Pa(-) is right-continuous at ¢ = 0 and 
imply that Pa(-) is continuous for all ¢ > 0. By (8) and (4), all the de- 
rivatives of Pa(-) exist for t > 0, and by (5) they are also right- 
continuous at t = 0. An application of Taylor’s theorem then yields 
(Ref. 6, pp. 240 ff.) 

ae | 


Pa(t) = exp (At) = py il ; Aa’, i2 0. (6) 


The elements of Ag situated immediately either above or below the 
diagonal are all strictly positive and so Ag can be symmetrized. In- 
deed, let 


Dz = diag [ 50, 61, ees ba | 
with 
ade 4 
= _ Bip2 Km a —4 = sae 
60 c and 6m ¢ (ae cin. ) EDm ) m 1, ) d, 


where (2) ¢ is a nonvanishing but otherwise arbitrary constant, (72) 
the p, are the equilibrium state probabilities, and en) £= (nj. 
Without loss of generality, we can—and shall—set ¢ = Pb so that 
&€=1 and 


Dz' = diag [pi, pi, ---, pil. (7) 
It is easy to verify that 
Sa = Dz'-Aa:Da (8) 
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is symmetric, its nonvanishing elements being 


Sam = — Oe sd Lm), m= 0, I; Rs d, (Xa = 0), 
Smymt1 = Smpiym = (Nmbm4i)}, m=0,1,-°>-,d—1. 


Hence, by (8) we have 
At = D,-S*-D7z!, k = 0,1, 2, <--, (9) 
and, by (6), 
Pi(i) = De-exp (Sal) 'Dz? = Y 4 Da Si-Dz). (10) 


k=0 

The representation of Az in terms of the symmetric tridiagonal matrix 
S.z entails substantial formal simplification of the final results. And it 
is also particularly convenient computationally, since the determina- 
tion of the characteristic values of Aa (which are needed for an exact 
solution) is best carried out after symmetrization. 

The matrices Az and Sz clearly have the same characteristic values, 
ro, T1, °°*, Ta. But Sa is symmetric and is therefore unitarily similar to 
the diagonal matrix 


Ca = diag [1ro, m1, +°-; Ta]. 
This means that an orthogonal matrix Bz exists such that 
Sa = By:Ca:Ba, Bz: Bi = Bi: Ba = Iu, (11) 


where B; is the transpose of Bu. 

But Sz is also tridiagonal, and its off-diagonal! elements never vanish. 
Hence, Sz is nonderogatory and its characteristic values are necessarily 
distinct (Ref. 7, p. 26). The elements in the nth column of B, are then 
the components of the (uniquely defined) normalized characteristic 
vector associated with the nth characteristic value r,(n = 0,1, ---, d). 

We now substitute (11) into (10). This yields 


Pi) = + A (Dy By-CH-By-Dr Dt", 
kao k! 
so that 
P(t) = Da-Ba-exp (Ca-t)-Ba-Dz' 
= D,-Bi- diag [e™*, e', ---, e*4¢]-Ba-Dzl. (12) 


We note now that all the row sums of Aq vanish and one of the 
characteristic roots, 7o, say, must therefore be equal to zero. Further- 
more, known extremal properties of the characteristic values can be 
used to show that ri, re, ---, 7a are negative. It is also readily seen that 


ps” = (pa, pi, a) pi)’ 


is the characteristic vector of Sq that corresponds to the vanishing 
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characteristic root 7». Indeed, let eg and 04 be the (d + 1) dimensional 
(column) vectors whose components are all equal to 1 and 0, respec- 
tively. Then, since Ag-@z = 0a, we have, by (8) and (7), 


Dz:Sa:Dz'-ea = Dz-Sa:ps? = Og. 
But none of the diagonal elements of Dag vanishes and the relation 
Da:Sa:ps” = 04 


can hold if and only if Sz-p{ = og. Thus, p{? is the characteristic 
vector associated with ro(=0), a fact that may be of relevance in the 
computations, as a comparison of p{? with Dy! provides an accuracy 
check for the method used to determine the characteristic vectors. 

In the derivation of formula (12), advantage was taken of the fact 
that the transition-rate matrix Ag is symmetrizable. It is worth 
noting that this relatively simple expression for Pg is a consequence of 
this property, and therefore holds for all (and actually only for) 
reversible Markovian processes with finite state spaces. Indeed, by 
definition, the class of these processes—which includes those of the 
birth-and-death type—is fully characterized by the following condi- 
tions (Refs. 8 and 9): 


DmP mnt) a DaPam(t), mn = 0, 1, Der 69 d, (13) 
or, equivalently, by the single relation: 
D;?-P, = P,: D7’. (14) 


Hence (12), written in terms of Sq, implies that 


Dz?-Pa = Djz'-exp (Sat)-Dz' 
(Dz'-exp (Sat)-Dz')’ = Pa-Da’, 


and (14) is therefore satisfied. 

Conversely, we show next that (14) is a sufficient condition for 
(12) to hold. 

Pre- and post-multiplication of (14) by Da yield 


D;!-Py:Dzg = Dz: P;:D7!. (15) 


I 


I 


Substituting the expansion of Pz as given by (6) into (15), and perform- 
ing the multiplications by Da and Dz’ under the summation sign 
(which is clearly legitimate), we obtain: 

St wp AEDot = Yt (Dy (A -DrDtH, t=O. 

K=0 k! Kook! 
However, this relation cannot be satisfied unless 


D;'!-Aa: Da = D,-Az:D;z', 
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so that, by transposition, 
(Dz!-Ag: Daz)’ = Dz!-Aa: Da. 


This means that Ag is symmetrizable by pre- and post-multiplication 
by Dz? and Da, respectively, and (12) then follows as shown earlier. 

Under the assumption that the process is reversible and that all the 
states communicate with each other” (1.e., pma(t) > 0, m,n = 0, ---, 
d, t > 0), the characteristic roots of Aa are necessarily simple. (Note 
that Ag, and hence Sg = Dz!-Aa-Du, need no longer be tridiagonal.) 
This can be proved as follows. 

The matrix Sz is symmetric and can therefore be tridiagonalized by 
a method from Householder (Ref. 7, pp. 152, 153, 290-2938, and 348). Ac- 
cording to this procedure, the tridiagonalization of Sa is achieved by 
successive right and left multiplications by symmetric orthogonal 
matrices, U1, U2, ---, Ua-1, of the form 


U, = Iq _ 2w,' Wy, 


where w, is a suitably chosen d + 1 dimensional (column) vector 
whose first r components are zero. (All the U, are of order (d + 1) 
and U? = Iu, r = 1, 2, ---, d — 1.) A direct application of the results 
derived in Ref. 7, above, shows that Sa admits of the following repre- 
sentation : 


SS Ui Usa Ugg Ps Ug ee UU 


where T, is a symmetric tridiagonal matrix of order d + 1. 

Let 6; be the elements of Ta. 

We are now faced with two possibilities. Either 6;,i41 = 0i41,; ¥ 0 
for 2 = 0, 1, ---, d — 1, or there is an index j(<d) such that 6,541 
= 03:1,; = 0 so that 

T; 0 
Ta = (16) 
O Taj-1 


In the first instance, all the characteristic roots of Sa, and hence of 
Aga, are distinct (Ref. 7, p. 26). To complete the proof, it is therefore 
sufficient to show that the second contingency cannot occur when all 
states communicate with each other. To this end, we proceed in- 
directly. We assume that (16) is satisfied for some 7 < d and show 
that some states then do not communicate with others. 

When (16) holds for 7 < d, we have, for any k = 0, 


A‘ = D,:Sk- Dz! 

T§ 0 
Die Ue ek Uy: MUgage seat sOqeDe ls 117) 
O° Taj 
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The first row in Dg is (po?, 0, ---, 0) and the elements in the first 
column and first row of each of the U,’s are zero except for their first 
component, which is always equal to 1. Hence, 


(1,0, «++, 0)-Dg-U1- --- -Ua-a = (po#, 0, «++, 0). (18) 
Similarly, since the U,’s are symmetric, we have 

Ua_1: +++ -U;-Dz!: (1, 0, ---, 0)’ = (pi, 0, ---, 0)’. (19) 
Hence, by (17) to (19), 


Poot) = (1, 0, +++, 0) (Spi k 4) (1,0, +++, 0) 
= = 79 Ot, (20) 


where 6§” is the element belonging to the first row and first column 
of T%. 

Let Qmn and Smn, m,n = 0,1, ---, d, be the elements of Ag and Sa, 
respectively. Under the present assumptions, Qmn‘Qnm 2O and 
Smn = (Amn*Anm)?, m, n = 0, 1, ---, d, m # n. The elements in the 
first r rows and columns of Sq are therefore uniquely determined by the 
elements in the first r rows and columns of Ag. Similarly, the vector w, 
depends only on the components, smn, of Sa for which either m S r — 1 
and n =r-—1, ---, d or, by symmetry, n S$r—1andm=r-—1, 

-, d (Ref. 7, pp. 290 ff). Consequently, the elements of T§ (which are 
all obtained after 7 — 1 steps) depend only on the elements of the 
first 7 + 1 rows and columns of Ag. This implies that the transition 
probability Poo(t), as given by (20), is independent of the rates @mn, ™m, 
n > j. However, the process being Markovian, this can only be true 
if Pom(t) = 0 for m > 7 which means (since, by assumption, 7 < d) 
that state 0 does not communicate with states 7 + 1, ---, d, as was to 
be proved. 


Ill. COVARIANCE FUNCTION 
3.1 First version 
The covariance function of the carried-load process is, by definition, 


Ret) = 5 (€ A n):(¢ A m)papam(t) — MY 


m,n=0 


I 


(CA n)-(c A m)palPan(t) — Pm 


m, n=0 
where 


Mi EN OSS eee 


n=0 
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However, if y is an arbitrary constant, the covariance of the process 
{N.(t) + y, t 2 0} is also R.. Hence, with the notation 


pn = (can) +7, nm =0;1, 5d, 
we also have 
d 
R.(t) _ ae PnPmPnL_Pnm(t) = eel (21) 
Let 
Do pi ares Da 
Pp, = sh Pi Pa 
Po Pi bere Da 
and 


Galt) = [pan(t) — Dm]. 
The matrix Pz can be obtained by letting {© in (13). Hence, 
P, = D,-B,-diag [1, 0, 0, ---, 0]-Ba-Dz! 
and 


Ga(t) = Pa(t) — Py 
= D,-B,-diag [0, e!, ---, e7¢]-By-Dz?. 


We now introduce two auxiliary row vectors: 
Ta = (po, pi, ***; pa), Sa = (Popo, P1p1, ***, Papa)- 
Then the coefficient of e7** in the linear form 
Si-Ga(t) Ta = Sq-Da:By- diag [0, e™', ---, e7#*]-By-Dy rq 
is the same as the coefficient of e** in (21), and we may conclude that 
R.(t) = sq: Da: Bz: diag [0, e"', ---, e74¢]-By-Dz?-ra. 


With the notation : 
da = (popd, ae | paPa)s 
we have 
qa = Sa-Da and qa = Dz? -ra, 
so that 
R.(t) = qa-Ba- diag [0, e"', ---, e"4#]-Ba-qa (22) 
or, alternatively, 


d 
R.(t) = >> bier (23) 
i=] 


with b; the ith component of the row vector qz:By. This last expression 
shows that the coefficient of e’** in either (22) or (28) is necessarily 
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nonnegative. Furthermore, since all the r;’s (¢ = 1, ---, d) are nega- 
tive, we also have 
d* 
(—D* aa R(t) 2 0, t20, k = 0, 1, uty 
and &, is therefore completely monotonic over [0, ~).° 

If we now set y = — c, the last d — c components of qa are equal 
to zero, and so the determination of R, by means of (22) necessitates 
only the computation of the first c components of the characteristic 
vectors. Formula (22) is therefore often well-suited to the case of delay 
systems. But unless the number of waiting positions exceeds the 
number of servers, greater reduction of the computations can be 
achieved by means of the formulas derived below. 

In the preceding derivation, the p’s are independent of the arrival 
and departure rates, and the formulas of this subsection therefore hold 
for arbitrary, reversible, Markov processes with finite state spaces. 
In contrast, the results of the next subsection are restricted to birth- 
and-death processes. 


3.2 Alternative forms 


Multiplying the nth equation in (2) by (c a n) and then summing 
with respect to n(0 Sn Sd), we obtain, after rearranging and 
canceling terms, 


d d c—l1 c 
> (c A N)*> Dmn = >» NxDmn —_ ys MnDmn.- (24) 
n=0 dt n=0 n=l 
But 
d 
d. (c A 2)P mn (tb) = E{N.(t)|N (0) = m}, 
so that, by (24), 
d c—l € 
qe tN Cd) | NO) =m} = 2 MnP mn _ 2 HnPmn- (25) 


Adding and subtracting cE {N.(t)|N(0) = m} on the right-hand 
side of (25), we obtain 


d oo 
qo iN |W (0) = m} = cE{N.(t)|N(O) = m} + 2 Pa(K)Dmn; 
m=0,1,---,d, (26) 
where 
(An — Bn — Kn) if tees Ope tee, (uo = 0), 
pn(k) =s = (ROP yi) ifn =, 
— KC ifm=ct+1,---,d. (27) 
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In the preceding formulas, « is an arbitrary real number that may 
be positive, negative, or null. We see later that the covariance formulas 
can occasionally be simplified by appropriate choices of x. 

Taking the initial conditions, 


E{N.(0)|N(O) = m} = (c a m), O<med, 
into account, the solution of (26) is 
E{N.(t)|N(O) = m} 
= (c A m)e* + ex(t—u) oS pi()Pma(u) | -du (28) 
so that 


R.(t) 


Il 


(6A m)pnB(Ne(®)|N(O) = m} — Ma 
M .ekt — M?, 
t d d 
+ fi em & (4 mpm |X ok? ma(w) [du 


where M,. = EN?2(0). By means of (13), the preceding relation can 
be expressed in a much more convenient form: 


R.(t) = Mset — M2, 
t d Ps d 
+ fi eo & plsdpn |S (A m)pam(w) fod 
0 n=0 m=0 
= M oe" = M2, 
t d 
rs ii exw) ¥ pa(k)pall{Ne(u)|N(O) = nj du. (29) 


Next, substituting (28) into (29), we obtain 
d 

R(t) = Meret — Mi + tet D0 (c A 2): pn(k) Dn 
n=0 


t d Me 
+ f eae pa(s) Dn f ED nm(v)-du-du. 
0 n,m =0 0 


Let R¢ be the Laplace transform of R, and pam that of Dam. The pre- 
ceding relation then yields 

2 
Ri(s) = aie _ Ma 


— K § 





1 d ms 
+ Coa: z, (c A N)-pnl(K)* Pa 
+ 5 eile): em(0) ps BEEK. (80) 
n,m =0 § kK) 


We know, however, that pan(t) is of the form 


d 
Dnm(t) = Pm + 2s Yami, 
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so that 
d d 
Ret) = Le (ec an)-(c A m)pn DY Yume. 
m,n=0 t=1 


This implies that the only poles of Rj are r:,7 = 1, --:, d, and that 
lim:.. &.(t) = 0. Taking these two facts into consideration, we see 
at once that (30) reduces to 


* Yami 1 
R.(s) = = 5 Pale) Pm(k) * Pn x (rr; — (7; — )? $s — 1; 
provided x ¥ 7;,7 = 1, 2, ---, d. And referring back to the derivation 
of (22), it is readily seen that 


ert erat 


R(t) = [ai(e) Y-Be- diag | 0  — (i — 6)?’ ey op | Beatle, 
KAT, a=1,---,d, (81) 


where 
Caa(k)]’ = Loo(x)- pd, «++, pa(x) pd). 


The modifications needed when x is equal to one of the characteristic 
roots are immediate. Let 


diag [ao, a1, ---, aa] 


be the diagonal matrix obtained by setting the jth diagonal element of 
diag [do, a1, --+, aa] equal to zero. Then if x = r; we must have, with 
some as-yet-undetermined constant a and o2 the variance of N.(0), 


R.(t) = (0% + aje* 
+LaWY Br-diag® [0, "i= 1, +a] Beate. 
But &.(0) = o2, so that 
a = — [aa(x)]’-Ba-diag® [0, (r: — x), 2 = 1, «++, d]-Ba-qa(x). 
Hence, 
R.(t) = o%e%t + [ai(oY'-Bi-diag® | 0 ee ein DP doy a] 
‘Ba-qa(x). (32) 


It should be noted that (32) is valid even if « #7r;,7= 1, ---, d, 
and that it should be used in the computations [rather than (31) ] 
whenever «x is “close’”’ to one of the characteristic roots so as to avoid 
accuracy losses (see below). Since « is arbitrary, one could always 
choose it so that it is not ‘‘close’”’ to any of the characteristic roots. 
But, as shown next, it is often preferable to select it in such a way as 
to reduce the amount of computation, and this, in turn, may dictate 
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the use of (82). As we have seen, 


ert erat 


R(t) = [aiWOY Be-diag |0, — =, eS | Ba-ailo. 


But, as noted previously, p{? is the characteristic vector associated 
with the vanishing root (ro), so that 


erit erat 
Consequently, (31) remains valid for all [qa(«)]’ of the form 
{ [eo (x) oF vy]: pi, ace Coa(k) = y]-p}}, 


where y is an arbitrary constant. [The same remark, of course, also 
applies to (82). ] 

We are therefore always at liberty to add the same constant to all 
the p;(«)’s. Under some circumstances, this degree of freedom, together 
with the one provided by the introduction of x, can be used to reduce 
the dimension of B’ and B: entire rows in B’ and the corresponding 
columns in B can be set equal to zero without affecting the computation 
either of (31) or (82), or of the variance of the switch-count load. It 
is relevant to note here that this reduction would be largely illusory 
were it not for the fact that the normalized components of any of the 
characteristic vectors can be obtained without having to compute 
other components of that vector (see below). 

According to the result of Section 3.1, the covariance can always be 
cast in a form that involves only the first c components of the character- 
istic vectors. But when the input and departure rates for 0 Sn <c 
are linear in n, the covariance can also be expressed in terms of the 
last d — c + 1 components of these vectors. Indeed, the rates are 
then of the form 


An = AN+N, 
pn + p’, m=O, 1). s-¢ 1, 


(0) -Bi-diag [ 0 ey | Bee = 0. (33) 


Ln 


so that 
An — Mn = (A — wn + QO! — p’). 


Hence, with « = \ — wand y = wp’ — NX’, (27) yields 


0 P= 0 ieee 
pnA— Bw) = 4(¢—Ae—wety’ —N ifn=e, 
(u—AJetp’ —N ifn=ct+]1,---,d. 


For the random (Poisson) and the quasi-random inputs, the p’s take the 
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following simple form whenever the service time is exponential with 
mean 1. 


(2) Random input (A, = a, n = 0, 1, ---): 


0 ifn = 0,1, ---,ce-—1, 
px(—1) = 4 —a ifn =, 
c—a ifn =ct+1,---,d. 
(22) Quasi-random input [N sources, \, = (V — n))\, 
n=0,1,---,N]: 
0 ifm = 0,1, ---,c¢—1, 
pal —-(1 + r)]=4 (ec —- N)vA in = € 


c+ (ec—N)r ifn=ct+1,---,d. 


I 


From the preceding developments, we see that the p’s can be chosen 
in such a way that the number of components of the characteristic 
vectors needed to express R, is the smaller of the two integers c and 
d —c-+ 1. In particular, in the case of loss systems, only the (¢ +1)st 
component of each vector is needed. 

The parameters x and y can also be chosen so that only the first 
ec + 1 components of the characteristic vectors actually enter in the 
expression of R.. This will be the case if we set « = we-1 — Ac—-1 and 
Y = €(Me-1 — Ac-1). 

In Ref. 3, the derivation of the covariance function for loss systems 
[d =c, N.(t) = N(t)] with Poisson input and exponential service 
time makes use of the differential equations 


BAN (®|N(O) se es BON) IMO) Son alt Sy, 


m=0,1, ---,¢. 


These equations appear here as that particular instance of (26) for 
which « = — 1,4, =a,n=0,1,---,c—Il,anduza=n,n=1,---,¢. 
Note also that now 3°{=§ pmn(t) = 1 — pme(t). But we stress that, in 
Ref. 3, the determination of the covariance relies on known recurrence 
relations between the so-called “‘sigma”’ functions (Ref. 3, pp. 129 and 
143 ff.) ; the more general problem considered in the present paper is 
not as readily amenable to such a treatment because of the greater 
complexity of the expressions that would now have to be used instead 
of the sigma functions. As we have seen, however, relatively simple 
formulas for R. can be obtained without extensive algebraic develop- 
ments as long as the underlying process is reversible. 
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IV. VARIANCE OF THE SWITCH-COUNT LOAD 


The variance of the switch-count load is now readily obtained. 
Depending on which expression we select for R,, we have either 


(i) Var L.(T) 


ne? = (n — |k|)Re(br) 


mtg By-diag [0,3 (n= |elyeni =a, a] 
‘Ba-da(x), (34) 

with 

q = (pops, a pc-1p}-1, 0, saa 0), or 


(it) Var L,(T) = n~[qa(«) |’ -Ba 


erilklt . 
diag |0,, me (n — |%|) ——— (r; @,— xp?’ a= 1, al 
‘Ba-qa(k), (38) 
with « # r;,7 = 1, ---, d, or 
(vit) Var L,(T) = n{o% — [aa(x)]'-Ba 
‘diag LO, (ri i 2 me Pea La d]-Ba-qa} 


-y (n — |k|)etl#lt + nL qa(k)]'-Ba 


k=—n 


ene it erilkle 
‘ (3) = ae 
diag“ oo (n — |k|) = Gon: 1, .a| 


‘Ba-qa(x), (36) 
where x ¥ r; for 2 ¥ 7. 
We now make use of the following identity (Ref. 3, p. 137): 


k |= e72nu 
>~ (n— {kl jee = n-cothu — ——3 —— " esch? u 
n=—k 


By means of this relation, (84) to (86) can also be written as 





= n—-a,-B.-di Sia ee 
Var L,(T) = n7-qq: Ba: diag |0, coth ( 9 ) on 





-esch? ( =) ,t=1,---, a| Base; (34a) 
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Var L,(T) = n--[qa(x) ]’-Ba 


: any Lees 
diag |0, ~—+, ay {coth ( 5 jy 
-esch? ( +) a on a a|-Bs- ails, 
KE Ti, 7=1,---,d; (85a) 





and 


Var L,(T) 
=n {o%—[qa(k) ]’-By-diag™ [0, (ri—«)~?, 1=1, ---, d]-Ba-qa(x)} 


-{eoth ( —) ae Fy esole ( F itn + n-[qa(«)] Ba 


ae 1 = Thi Lent 
: (3) essere soe oR (ears eee ol 
diag“ lo, ( — 6) {eoth ( 9 ) on 


-esch? ( =) ,t=1,°-°-, a| Bs ails (36a) 

















where « ¥ r; for? # j. 

Let Var L,,(T) = lim... Var Ln(T) be the variance of the load 
measurement obtained by continuous observation of the number of 
busy servers. If we replace 7 by 7'/n in (34a) to (36a) and then let n 
tend to infinity while keeping 7’ fixed, we obtain the following formulas: 








Var L..(T) = :. Qa: Bi 
-diag ° (1 qt na ),i- ide. a | Beas, (34b) 
Var L(T) = — % aC ONE ‘Bi 
. 1 1 — ev? 5 * 
- diag ae Tr, ),i- 1, +d | Bo-ails 
fees me ee, eee eee) 
and 


2 * , 
Var L,.(T) = 7 {oz — [qa(x) ]’-Ba 
-diag® [0, (r; — x), 7 = 1, ---, d] 
* I 1 —? «T 2 * ? , 
‘Ba-qa(x)} =( i ) ~ 7 Caa(x) |’-Ba 
‘ : 1 1 — en? , * 
‘ qs) eee = UNG j . 
diag“ |0, Fie = (1 + Tr, ) ,t=1, ; a| Ba-qa(k), 
K #1; for tj. (6b) 
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We note that the formula for the variance of sums of dependent 
random variables makes it possible to compute the covariance between 
load measurements performed over distinct time intervals. Indeed, 
consider for instance a sequence of load measurements over the 
intervals (0, T], (T, 2T], (2T, 3T], ---. Let L@(T) be the switch- 
count load over the 7th interval (¢ = 0, 1, ---), S.() = n?L,(t), and 
TY = Cov [LO (T), LE (T)]. Then 


Var Snceryl(k + 1)T } _ (k + 1) Var S(T) 


425 (+1 —dwTO(D) 
1=1 


and 


eto et 


Tre) = Var Lauspl + 1)T | — —~— Var LCT) 


~ ¥ recy. 


i=1 


The preceding formulas may be used to determine the T (7) re- 
currently. But the results of such computations shall be exact only if, 
for some choice of the time origin, all the scanning instants are multiples 
of r. 

We conclude this section with the remark that the variance formulas 
(34), (84a), and (84b) are valid for arbitrary reversible Markov pro- 
cesses with finite state spaces. 


V. NUMERICAL CONSIDERATIONS 


The exact variance formulas of the preceding section are very well 
suited to electronic computation and are easily programmed since, 
apart from straightforward evaluations of hyperbolic functions and 
simple products of matrices and vectors, they only involve the deter- 
mination of characteristic values and vectors for which powerful sub- 
routines are readily available. The fact that Sa is symmetric and tri- 
diagonal (or reducible to tridiagonal form by an orthogonal similarity 
transformation) allows us to use the subprogram TQL2, which is par- 
ticularly efficient under the present circumstances (Ref. 11, pp. 227- 
240). Without going into details, we mention here only that this sub- 
program is based on the so-called QR-algorithm and relies on the con- 
struction of a sequence of symmetric tridiagonal matrices, S{”, 


n = 1, 2, ---, unitarily similar to Sz, which converges to diag [0, 71, 
Te, +++, Taj]. At the nth iteration Sf” is expressed as a product of an 
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orthogonal matrix Q§” and a lower-triangular matrix L”: 
si = Qs? -L?”. 


In Ref. 11, this decomposition is carried out by the Givens’ triangular- 
ization in which the Q§”’s are expressed as products of simple plane 
rotations (Ref. 7, pp. 239-240). 

The (n + 1)st iterate of Sz is then given by 


Spt? = LP-Qp, 


whose unitary similarity to Sf” (and hence to Sa) follows from the 
relation 


LY OF” = (Q%?)’ -S¥? -Q””. 


This method avoids the numerical difficulties frequently associated 
with the computation of the zeros of the characteristic polynomial. 
As shown in Ref. 11, p. 228, 


Qs” = RP. --- -RP, (37) 
where R is of the form: 
1 
C. 8: 
—S; C; <— row 1. (38) 


1 


According to Ref. 11, p. 231, the matrix Bz; (whose columns are the 
characteristic vectors of Sq) is given “almost to working accuracy”’ by 


Qa -- - Qa, (39) 


where n ($380) is the number of iterations needed for the (numerical) 
symmetrization of Sa. Taking (37) to (39) into account, it is then 
readily seen that the elements of B; can be determined row by row 
which, as we have remarked earlier, is a desirable feature in the present 
context. 

Computations have been carried out for systems having as many as 
400 devices (and hence transition-rate matrices of order 401) to deter- 
mine the numerical accuracy of the approach described in the preceding 
sections. Checks were performed by comparing the value of Var Li(T) 
obtained by means of (35a) or (36a) with the corresponding o7, which 
can be calculated directly and independently from the equilibrium 
state probabilities. These two quantities, which are theoretically equal, 
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turned out in all cases to agree to at least 10 decimal places with the 
greatest difference occurring when d was largest. Hence, our procedure 
indeed yields very accurate results for the type of systems that are 
likely to occur in practice. But when d is large, the storage require- 
ments and the amount of computations become critical. It is therefore 
always important to select «x and y in such a way as to minimize the 
number of B’ rows that actually enter into the computations. (It 
follows from earlier remarks that this number, for proper choice of «x 
and y, never exceeds the integral part of (d + 1)/2.) Further reduction 
can also be achieved by excluding the states whose probabilities of 
occurrences are so small that neglecting them will not materially affect 
the final results. In this connection, we make the following remarks. 

The variance of the switch-count load is perturbed by at most 


[oi(k) P-pj-o% 


if p; is set equal to zero in the particular formula used to evaluate 
Var L,(T). Hence, since 


Var L,(T) = o2/n, 


we always have the following upper bound for the relative error, «,;, 
induced by setting p; equal to zero: 


ej < Cox(k) 2+ pin, J = 0, L, ae d. 


For a given relative accuracy of Var L,(T), these inequalities make it 
possible to determine ahead of time whether some components of the 
characteristic vectors can be “safely” eliminated from the computa- 
tions. In large systems, the gains achieved by such a reduction may be 
quite substantial, as either low occupation states [N (t) small] and/or 
high occupation states [NV (t) large] have then frequently very small 
probabilities of occurrences. 

Computations could be arranged to determine only those character- 
istic roots that-are required to reach a given degree of accuracy 
[plus those needed to compute Ba-qz(x)]. This is rather readily 
achieved in loss systems with Poisson input and exponential service 
times since, in this case, the coefficients b? of 


mri\ 1 — entri of Tri 
coin ( HZ 5 ‘) “So eso! ( 2 ) 


in the variance formulas of Section 3.1 are then monotonically de- 
creasing as |7;| increases: 


bj < bj if [rs| > Irsl, 2,9 = 1, “nny d, 
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RELATIVE AMPLITUDE 





0.01 
0.001 
0.0001 
0 10 20 30 40 50 60 70 80 
ROOT INDEX 


Fig. 1—Spectral measure of the carried load process. 


But, in general, one would encounter an additional difficulty, namely, 
that the 6?’s do not have the monotonicity property alluded to above 
and may actually fluctuate widely. This is illustrated in Fig. 1, where 
the roots are assumed to be indexed in order of increasing magnitude 
and the ordinates are the corresponding 6?’s, normalized in such a way 
that max; b? = 1. 

The computations should be based on (35a)—(36a) or on (35b)—(86b) 
in the case of continuous measurements—as these formulas provide 
us with all the flexibility needed to cut down both storage space and 
computation time. When choosing between (35a) and (86a) or between 
(85b) and (86b), one should keep in mind that, for x close to 7:, the 
difference r; — x may not be determinable with enough precision to 
allow accurate computation of Var Z,(T). This is shown in Table I 
where x = — 1 and 1, is the root of smallest positive absolute value. 
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Table |— Loss system, 80 servers, Poisson input, 
exponential service 


Offered Var Li(T) 

Load in l+n oe ce 
Erlangs Formula (36a) | Formula (35a) 
10 1.09 x 10-® 10.000000 10.000000 0.034614 
20 —9.90 X 107% 20.000000 20.000000 0.013173 
30 —1.14 X 10 30.000000 30.000000 20.891365 

40 —3.38 X 


10-6 39.999986 39.999986 39.999926 


(Note that the last two columns of this table should be equal and that 
errors of the same magnitude would arise if one were to use (18) of 
Ref. 3.) In all our computations, we have made use.of (35a) and (36a) 
whenever |x — r;| < 10~‘ for some 7. This bound for |x — r;| is both 
large enough to ward off appreciable accuracy losses and small enough, 
under prevailing conditions, to be satisfied by only one root. | 


VI. REMARKS ON INFINITE SYSTEMS 


It is known that infinite systems can be regarded as limits of finite 
ones,” and it is therefore of practical interest to have information 
concerning the spacing of the characteristic values as the dimension, 
d, of these approximating systems becomes large. Indeed, as d increases, 
computational difficulties may arise because of-a lack of separation 
between these roots. Such problems would ‘certainly come up sooner 
or later if the spectrum of A = limz.,,.Aa happens to be dense over 
some interval as, for instance, in the case of a single-server queue with 
Poisson input, exponential service timé, and unlimited waiting room 
(Ref. 12, pp. 365-366). Infinite systems with well-separated roots do, 
of course, also occur. As an example of this type, we mention the 
systems with an infinite number of servers, Poisson input, and ex- 
ponential service which often provide useful idealizations. (In this 
case, as is well known, the nonvanishing characteristic roots are the 
negative integers, —1, —2, —3, ---.) Other examples of systems with 
discrete spectra are given in Ref. 12, where sufficient conditions for 
this to occur are discussed in some details; but in all these instances 
the A, and yu, both increase as n’ for some v > 0. This condition is un- 
likely to be satisfied in queuing systems; generally, in this particular 
area, the arrival and the departure rates remain bounded: 


05r.SA<~, O0OSnr.5M <a, n=0,1,---. (40) 
As briefly described below, these inequalities imply the existence of 


definite bounds for the spectrum of A. 
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Consider an infinite system, and let A be its (infinite) transition-rate 
matrix. Let Az be the matrix obtained by retaining only the elements 
belonging to the first (d + 1) rows and columns of A and then setting 
Aa equal to 0. Let rao(=0) > rai > --+ be the characteristic roots of 
Aq. Then, under conditions (40) it can be shown that, for any k 2 0: 


(2) lraz] <A+M for d sufficiently large, 
(it) |raa-z] < 2(A + M) for d=k. 


Hither of these two inequalities implies that the characteristic roots 
do not remain separated as d > whenever (40) is satisfied. Under 
the more stringent requirements that (40) holds and that 

lim \, = A, lim pa = M, 


more precise statements can be made, namely, that, for all k’s and d’s, 
lraz| < (VA + VM) 
and that the spectrum of A always comprises a closed interval, viz., 
Q = [—(VA + VM), —(VA — VM)?}. 


(In addition to ©, the spectrum of A may also include a finite number of 
roots in [—(VA — VM)2, 0].) But it turns out (as will be shown 
elsewhere) that, as d increases, the characteristic roots of Aga fill Q 
rather ‘“‘evenly’’; furthermore, for practical accuracy levels, large 
values of d are needed only when the length of 2 tends to be relatively 
large (a circumstance corroborated by extensive computations). 
Hence, within the present framework, it appears that root-spacing is 
not likely to be critical except in the improbable event that extreme 
precision is required. 
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T1 Carbon Transmitter Model for Use in 
Computer-Aided Analysis of Telephone 
Set Transmission Characteristics 


By D. R. MEANS 
(Manuscript received March 5, 1975) 


A carbon transmitter model is presented, the purpose of which is to serve 
as a tool for computer-aided analysis of telephone set transmission 
characteristics. The derivation of the model is based upon the physical 
theory of the device. The parameters in the model are evaluated by com- 
paring the analytically derived expressions for the device characteristics 
to the measured characteristics of a typical device. Because these param- 
eters are related to the physical theory, the model not only serves its de- 
sired practical end, but also serves as a vehicle whereby an understanding 
is obtained of the relationship between device characteristics and physical 
theory. 


1. INTRODUCTION 


Computer-aided optimization of the transmission characteristics 
of telephone sets requires that accurate models be obtained for all 
transmission-related telephone set components. A carbon transmitter 
model has been derived for this purpose. This model has been used 
in a telephone set transmission analysis computer program, and good 
agreement between computed and measured transmission charac- 
teristics was obtained. 

The de V-I characteristic of the carbon transmitter is nonlinear. 
This nonlinearity must be taken into account in the de model so that, 
in the transmission analysis program, the operating point of the 
transmitter can be determined, as well as the operating points of any 
nonlinear telephone speech network components, e.g., silicon carbide 
varistors. Thus, the de model is a voltage-dependent resistance. 

The ac model is similar to that of a vacuum-tube triode, consisting 
of a Thevenin-equivalent resistance and voltage source. The voltage 
source is dependent on the amplitude of the force acting on the carbon 
granules because of the acoustic excitation of the transmitter. This 
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force is in turn the output of a filter which represents the transmitter 
structure itself. The input to the filter is the acoustic sound pressure. 
This filter can be represented by an electrical equivalent circuit,! the 
derivation of which is straightforward. However, it is the charac- 
terization of the effects taking place within the carbon chamber which 
is of primary interest here. The filter is represented simply by its 
measured frequency response. 

Values for the various parameters of the model were determined 
by comparing the expressions derived analytically for the charac- 
teristics of the transmitter to the measured characteristics of a typical 
device. All measurements were made with the transmitter face in a 
vertical plane, in a telephone handset, and in a position relative to the 
artificial mouth as specified in IEEE Standard 269-1971.? Also, the 
transmitter was mechanically and acoustically conditioned prior to 
the measurements. The acoustic conditioning signal was swept be- 
tween 300 and 3300 Hz at a rate of six sweeps per second and fre- 
quency-weighted corresponding to the average sound pressure spec- 
trum of continuous speech, and had an average sound pressure level 
of 94 dB (re 0.0002 dyn/cm?). The conditioning signal was applied 
for 3 s. 


ll. DC MODEL 


As shown in Fig. 1, the de V-I characteristic of the carbon trans- 
mitter is nonlinear. Goucher’ attributed the nonlinearity to the effect 
of joule heating on the contact resistance between carbon granules. 
Later, Mol,* disputing Goucher, attributed the nonlinearity to the 
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Fig. 1—Direct-current V-I characteristic. 
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effect of electrostatic forces between carbon granules. On this basis, 
he derived an expression for the dc transmitter resistance as a function 
of voltage and found it to agree fairly well with measured transmitter 
resistance. However, Mol’s expression for the de voltage dependence 
of the resistance variation when an acoustic signal acts on the trans- 
mitter agrees rather poorly with measurement. This casts doubt on 
the electrostatic force theory. In fact, experimental results indicate 
that the effect of electrostatic forces between carbon granules is 
negligible. These experiments are described in Appendix A. On the 
other hand, a more recent study of the theory of electric contacts 
tends to support Goucher’s theory. Holm‘ treats the subject of the 
effect of joule heating on contact resistance extensively, and his 
results will be applied to the derivation of the transmitter model. 

As will be seen, the nonlinearity of the V-I characteristic cannot 
be accounted for entirely by the effect of joule heating on contact 
resistance. The effect of the thermal expansion of the carbon chamber 
due to joule heating must also be considered. This effect is readily 
demonstrated experimentally because of the relatively long time 
constant involved. If the transmitter current is changed abruptly, a 
slowly decaying voltage transient is observed, owing to the hysteresis 
associated with the expansion of the carbon chamber. The time 
constant is approximately 1 s. 

A cross section of the T1 transmitter is shown in Fig. 2. The carbon 
chamber consists of a movable dome electrode connected to a fixed 
conical back electrode by a flexible, nonconducting chamber closure. 
According to the results of Fritsch’s analysis® of the thermal response 
of the transmitter structure, the transient effect is due primarily to 
the expansion of the dome electrode. As the dome electrode expands, 
it compresses the carbon granules, lowering their resistance. Fritsch 
called this effect ‘‘thermal packing.”’ Of course, after the transient 
has decayed, the transmitter can be reconditioned to unpack the 
granules. However, following the reconditioning, a new thermal 
equilibrium will be established so that some degree of thermal packing 
will still occur. This effect, as well as the effect of joule heating on 
contact resistance, must be included in the model. An understanding 
of these effects is based on a consideration of the factors affecting 
the resistance of a single carbon contact. 


2.1 Effect of contact force on contact resistance 


The contact resistance between two carbon granules is related to 
the magnitude of the force pressing the granules together by the 
expression : 

r, = KP oa (1) 
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Fig. 2—Cross section of T1 transmitter. 


where P is the contact force and K and y are constants. Equation (1), 
which is generally valid for electric contacts, was shown by Goucher’ 
to be valid for carbon granules, based on direct measurement of 
contact resistance and force. Later, this equation was indirectly shown 
by Joscheck’ to be valid based on measurements of the bulk resistance 
of carbon granules as a function of the filling height of the granules 
in the measurement chamber. 

In the carbon chamber in the T1 transmitter, the significant forces 
acting on the carbon granules are the gravitational force resulting 
from the weight of the granules themselves, the force resulting from 
the thermal expansion of the dome electrode, and, when the trans- 
mitter is acoustically excited, the force resulting from the acoustic 
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pressure acting on the diaphragm and dome electrode. Because of the 
random orientation of the contacts throughout the chamber with 
respect to the directions of the applied forces, the contact force and, 
hence, the contact resistance will also vary randomly. However, 
because of the large number of contacts involved, the random variation 
in contact resistance can be ignored. Only the larger overall gradients 
in bulk resistivity will be considered. 

For simplicity, an approximate chamber geometry is assumed in 
which the two electrodes are concentric hemispheres. A spherical 
coordinate system is defined such that the hemispherical chamber 
walls lie at constant radial distances from the origin, the inner and 
outer radii being designated a and b, respectively. The component of 
contact force resulting from the weight of the granules themselves is 
represented by P,(r,6,¢). The component resulting from the 
thermal expansion of the dome electrode, being proportional to the 
power dissipated by the transmitter, is represented by (V?/R) 
Pa(r, 6, @), where V and RF are the transmitter de voltage and resistance, 
respectively. The component resulting from the acoustic excitation 
of the transmitter, being a function of time as well as of the spatial 
coordinates, is represented by AP(r, 8, @, ¢). 

Equation (1), therefore, becomes 


ri (7, 6, g, t) = K[P, Ge 6, ¢) + (V?/R)Pa(r, 0, ¢) 
+ AP(r, 0, 4,0)" (2) 


It will subsequently be seen that, for normal speech signal levels, 
AP is small enough compared to the static components of contact 
force that its effect on the de component of contact resistance is 
negligible. Thus, the relationship between the de contact resistance 
and the contact force is 


re(r, 8,6) = KLP,(r, 6, 6) + (V?/R)Pa(r, 8, 6) 7. (3) 


Although a change in de transmitter resistance is observed when 
acoustic excitation is applied to the transmitter, this is judged to be 
due to the effect of the acoustic excitation on the state of compactness 
of the carbon granules rather than to the effect of the nonlinearity 
of the contact resistance-contact force characteristic. The effect of 
acoustic excitation on the state of compactness of the granules will be 
discussed further. 


2.2 Effect of joule heating on contact resistance 


According to Holm’s analysis’ of the effect of joule heating on 
contact resistance, if certain assumptions regarding the temperature 
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dependence of the electrical resistivity and the thermal conductivity 
of the contact members are satisfied, then the effect of joule heating 
can be accounted for by multiplying the contact resistance by the 
factor 

n(Vi) = [B + (1 — B)(Ax/Vi) tan (Vi/ Ax) 5, (4) 


where V; is the contact voltage and A, and B are constants. Since 
the random variation in contact voltage is of no concern in the model, 
V, is considered to be the average contact voltage, which is propor- 
tional to the total transmitter voltage. Then 


n(V) = LB + (1 — B)(A/V) tan? (V/A) }", (5) 


where A is also a constant. 
The assumptions upon which the derivation of eq. (4) is based 
are that the thermal conductivity satisfies 


A = Ao(1 + BAT), (6) 
and the electrical resistivity satisfies 
p = po(l + eAT)/(1 + BAT), (7) 


where po, Ao, ¢, and 8 are constants, and AT7' is the change in tempera- 
ture because of joule heating. Apparently, these assumptions are 
valid in this case because of the excellent agreement between eq. (5), 
using the values for A and B, listed subsequently, and the measured 
data presented by Hufstutler and Kerns? for the resistivity of granular 
carbon contained in a quartz test chamber having a negligible thermal 
expansion coefficient. This implies that, aside from the effect of 
chamber expansion, the effect of joule heating on contact resistance is 
alone sufficient to account for the nonlinearity of the V-I characteristic. 

Now, if eqs. (5) and (3) are combined, the expression for the de 
contact resistance becomes 


re(r, 8, 6) = Kn(V)[P,(r, 8, 6) + (V?/R)Pa(r, 8, 6) 17. (8) 


2.3 Total de resistance 


If r,(r, 6, 6) is the contact resistance and there are n contacts per 
unit length, then the resistivity of the carbon granules is r:(r, 0, 6)/n, 
where r;(7, 6, @) is given by eq. (8). Then, for the approximate chamber 
geometry which has been assumed, the total de transmitter resist- 


ance is 
b 
R= i a (9) 


° [” [Cn sin o/re(r, 8, 6)]aodo 
0 0 
The mean value theorem can be applied to perform the integrations 


1306 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1975 


with the result 


R= Ro(V)LL +aV?/RY, (10) 

where 
Ro = K(b — a)/[nz’?* sin $P,"(7, 8, $) ], (11) 
a= Pal, 6, 6)/P,(F, 6, $), (12) 


and where 7, 6, and ¢ are constants, being the coordinates of some 
point within the chamber. Note that Ro is the limiting value of R as V 
approaches zero, and that aV?/R is the average ratio of the component 
of contact force due to thermal expansion of the carbon chamber to 
the component of contact force due to the weight of the granules 
themselves. Although FR is not expressed as an explicit function of V, 
a solution to eq. (10) can be obtained using iterative techniques easily 
implemented on the computer. Values for the parameters Ro, a, A, B, 
and y will be determined by fitting eq. (10) simultaneously with 
equations for the transmitter ac resistance and open circuit output 
voltage to measured data. The ac resistance and open circuit output 
voltage will now be considered. 


Ill. SMALL-SIGNAL AC RESISTANCE 


Over the range of frequencies of interest for speech transmission, 
the transmitter ac impedance is purely resistive. However, as is 
obvious from Fig. 3, the ac resistance is not the slope of the de V-I 
characteristic except in the limit as V approaches zero. This is ex- 
plained by the fact that, because of the large hysteresis effect, the 
thermal expansion of the dome electrode cannot follow the ac signal, 
at least not at frequencies above a few hertz. Thus, the difference 
between the ac resistance, which is not affected by the thermal expan- 
sion of the dome electrode, and the slope of the de V-I characteristic, 
which is affected, increases as the power dissipated by the transmitter 
increases. 

Because the thermal expansion of the dome electrode has no effect 
on the ac resistance, the ac resistance is the slope, not of the actual 
V-I characteristic, but of the V-I curve defined by setting the term 
equal to zero which accounts for the expansion of the dome electrode. 
This is the curve defined by 


T= V/[Ro(V)J, (13) 


where Ro and 7(V) are defined by eqs. (11) and (5), respectively. Thus, 
the ac resistance is given by 


rae = Rutt (V) /[a(V) — V ap a(¥) | 
= Ry(A? + V2)/(A? + BV). (14) 
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Fig. 3—Alternating-current resistance characteristic and slope of de V-I characteristic. 


It has been assumed that the ac signal level is low enough so that 
there is no significant joule heating effect due to the ac signal. This 
assumption is valid for signal levels typical of speech transmission. 
For higher signal levels, the ac joule heating effect will cause the ac 
as well as the dc resistance to decrease, as is easily verified 
experimentally. 

At frequencies low enough that the period of the ac signal becomes 
significant compared to the time constant associated with the thermal 
response of the dome electrode, the transmitter ac impedance exhibits 
a reactive component due to the effect of the thermal hysteresis. 
Figure 4, drawn from a photograph of a storage oscilloscope trace, shows 
the effect of the thermal hysteresis in the response to a sinusoidal 
driving voltage at a frequency of 0.2 Hz for four different operating 
points. The effect becomes greater as the dc bias increases, as would 
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Fig. 4—Thermal hysteresis effect at 0.2 Hz. 


be expected since the expansion of the dome electrode is proportional 
to the power dissipated by the transmitter. An expression for the low- 
frequency transmitter impedance that accounts for this effect is 
derived in Appendix B. For frequencies greater than approximately 
10 Hz, the reactive component of the transmitter impedance becomes 
negligible, and the expression derived in Appendix B reduces to the 
expression given by eq. (14). 


IV. OPEN-CIRCUIT OUTPUT VOLTAGE 


When the transmitter is acoustically excited, the contact force will 
vary owing to the effect of the acoustic pressure acting on the dia- 
phragm and dome electrode. The variation in the contact force at the 
point (7,6, ) is designated AP(r, 6, ¢,t). Then the resistance will 
vary by an amount AR(¢) such that 


R+ AR() = Ro (Vl + aV2/R + APO}, (15) 


where 


AP(t) = AP(;, 6, 6, t)/Py(F, 4, 6). (16) 


T1 CARBON TRANSMITTER MODEL ~~ 1309 


The transmitter voltage will change by an amount AV (é) where, if 
the transmitter current J is held constant, 


V+ AV() = I[R + ARQ) ]. (17) 
Then 
AV (t) = IAR(t) 
= VAR(t)/R, (18) 
which is the ac open-circuit output voltage. From eq. (15), 
AR()/R = (1 + AP(@)/( + aV?/R)T7 — 1, (19) 
so that 
AV (t) = Vill + AP(@)/( + aV?/R) 7 — 1}. (20) 


Although this is a nonlinear relationship, AP(t) will be found to be 
small enough compared to 1+ aV?/R at normal speech levels so 
that a linear approximation is justified. Thus, 


AV (t) & — yVAP(t)/(1 + aV?2/R). (21) 


Note that eq. (21) is —yV multiplied by the ratio of the dynamic 
to the static forces acting on the carbon granules. 


V. EVALUATION OF MODEL PARAMETERS 


The parameters Ro, a, y, A, B, and AP(t) were evaluated using 
an iterative optimization computer program to fit eqs. (10), (14), 
and (21) simultaneously to measured dc resistance, ac resistance, and 
open-circuit output voltage, respectively, the latter being measured 
at a frequency of 1 kHz with a sound pressure level of 94 dB (re 
0.0002 dyn/cm?). The measurements were performed on a T1 trans- 
mitter considered to be a typical unit. The resulting parameter values 
are listed in Table I. Of course, transmitter characteristics are subject 
to such factors as aging, temperature, conditioning, orientation, and 
manufacturing variations. The values of the parameters in eqs. (10), 
(14), and (21) will vary accordingly. 


Table | — Parameter values 


Parameter Value 
Ro 111.0 
a 0.94 
¥ 0.43 
A 7.12 
B 6.65 
AP (rms) 0.20 
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Fig. 5—Direct-current resistance characteristic. 


Equations (10), (14), and (21), with the parameter values listed in 
Table I, are plotted in Figs. 5, 6, and 7, respectively, along with 
measured data points. The agreement between calculated and measured 
data is judged to be within the limits of measurement error. 


Vi. INPUT-OUTPUT AND FREQUENCY RESPONSE CHARACTERISTICS 


To complete the model, the input-output and frequency response 
characteristics of the transmitter must be specified. The input-output 
characteristic is nonlinear owing to the effect of the acoustic excitation 
on the compactness of the carbon granules. As the acoustic signal 
level increases, the carbon granules are agitated into a less compact 
state and the mechanical impedance of the granules decreases. There- 
fore, the transmitter efficiency increases as the acoustic signal level 
increases, resulting in an input-output characteristic having a greater- 
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Fig. 6—Alternating-current resistance characteristic. 


than-unity slope. Because of this effect, the transmitter is able to 
discriminate against distant sounds, and thereby reduce interference 
resulting from background noise. However, the weak components in a 
composite signal such as speech are not discriminated against, since 
the compactness of the carbon granules, which is controlled by the 
strong components, is the same for all components of the signal. 
Therefore, the nonlinearity of the input-output characteristic does not 
affect the components of individual speech sounds, and the compact- 
ness of the carbon granules varies only as the overall energy content 
of the speech signal varies. 

As discussed by Bryant,® the frequency response of the transmitter 
is related to the nonlinearity of the input-output characteristic, since 
it also depends on the mechanical impedance of the carbon granules. 
This implies that the frequency response depends on the nature of the 
input signal. The response to a swept frequency sinusoidal signal is 
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Fig. 7—Alternating-current open-circuit output voltage characteristic at 1 kHz 
for a 94-dB (re 0.0002 dyn/cm?) acoustic input level. 


somewhat different from the response to a speech signal. (For further 
discussion, see Ref. 9.) 

For the purposes of the model, a continuous, random, speech input 
signal is assumed. Accordingly, the nonlinearity of the input-output 
characteristic is ignored, and the frequency response is measured as 
suggested by Bryant. The response is shown in Fig. 8 plotted relative 
to the 1-kHz output level. This is the response of the transmitter 
structure itself to the acoustic input signal, or, in the model, the 
response of the input filter. In the computer program, the filter response 
was stored as a table of values at discrete values of frequency. 
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Fig. 8—Transmitter gain (relative) vs frequency. 


Vil. CONCLUSION 


A carbon transmitter model has been presented. The physical theory 
upon which the model is based is supported by the close agreement 
between the calculated characteristics of the model and the measured 
characteristics of an actual device. Thus, it can be concluded that the 
nonlinearity of the de V-I characteristic is due primarily to the effect 
of joule heating on contact resistance and to the effect of the thermal 
expansion of the dome electrode due to joule heating. The effect of 
electrostatic forces is negligible. Furthermore, the difference between 
the ac resistance and the slope of the de V-I characteristic is due to 
the hysteresis associated with the thermal expansion of the dome 
electrode. Finally, the relative resistance change due to the acoustic 
excitation of the transmitter decreases as the de voltage increases 
due also to the effect of the thermal expansion of the dome electrode. 
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APPENDIX A 


In this appendix, the question is considered of whether electrostatic 
forces between carbon granules have a significant effect on transmitter 
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Fig. 9—Phase of diaphragm displacement resulting when transmitter output is 
driven by sinusoidal ac voltage. 


characteristics. Two experiments, the results of which indicate that 
they do not, are described. In both experiments, the transmitter 
output is driven with a sinusoidal voltage, i.e., as if the transmitter 
were a receiver. It is observed that, if the transmitter is driven hard 
enough, an audible signal is generated. Approximately 2-V rms is 
required for the signal to be audible at 1 kHz. This phenomenon could 
be the result of attractive forces between carbon granules owing to 
electric fields or the result of thermal expansion, presumably of the 
carbon granules, since the thermal inertia of the dome electrode and 
diaphragm is too large for their thermal response to follow the instan- 
taneous voltage at frequencies above a few hertz. The results of the 
two experiments which are now described indicate that the forces are 
due to thermal expansion of the carbon granules. This effect is insignifi- 
cant compared to the predominant thermal effects that are accounted 
for in the model. 

In the first experiment, the phase of the dome electrode displacement 
was measured. An outward displacement, in phase with the square 
of the driving voltage at frequencies far enough below resonance so 
that the mass of the system can be ignored, would indicate that the 
force was due to thermal expansion, while an inward displacement 
would indicate that the force was due to electric fields. The displace- 
ment was measured using an optical proximity detector. A small mirror 
was mounted on the dome electrode to provide a flat reflecting surface 
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Fig. 10—Output sound pressure—input power characteristic resulting when trans- 
mitter output is driven by sinusoidal ac voltage. 


for the detector. This was necessary to obtain measurable detector 
output, the displacement being very small. The phase of the outward 
displacement relative to that of the square of the driving voltage is 
plotted as a function of frequency in Fig. 9. The phase angle ap- 
proaches zero at low frequencies, indicating that the force acting on 


0.03 


0.02 


SOUND PRESSURE IN DYN/CM2 
° 
fan) 


0 5 10 15 20 25 30 
VOLTAGE SQUARED IN VOLTS? 


Fig. 11—Output sound pressure—input voltage squared characteristic when trans- 
mitter output is driven by sinusoidal ac voltage. 
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the dome electrode is due to the thermal expansion of the carbon 
granules. 

Since this is the case, the acoustic sound pressure level should be 
proportional to the power dissipated by the transmitter rather than 
the square of the voltage, as would be the case if the forces were due 
to electric fields. (Power is not proportional to the square of the 
voltage, since the resistance is a function of voltage.) This was verified 
by the second experiment, the results of which are plotted in Figs. 10 
and 11. 


APPENDIX B 


In this appendix, an expression is derived for the low-frequency 
impedance of the transmitter, taking into account the effect of the 
thermal expansion of the carbon chamber. 

If the transmitter current is changed abruptly by an incremental 
amount AJ, the displacement of the dome electrode due to the addi- 
tional power dissipation will lag behind the change in current due to 
the thermal hysteresis. According to Fritsch’s analysis,® the transient 
can be expressed as an infinite sum of decaying exponentials. Thus, 


AV (t) = (m+ D kee") ATu(l), (22) 


where m is the slope of the de V-I characteristic and where k; and r; 
are constants. Because the initial change in voltage must be 1r,.A/, 


> Ke = Tee MN. (23) 
i=1 


In the frequency domain, 


AV(s) = [m+ y kis/(s + 1/1.) JAI/s, (24) 
from which 
Z(s) = m+ ps kys/(s + 1/7;). (25) 


As pointed out by Fritsch,° the infinite series solution converges too 
slowly for practical evaluation. A practical expression for ac impedance 
can be obtained by assuming a single time-constant approximation 
for the transient response. Then 


Z(s) Bm + ks/(s + 1/7), (26) 
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where 7 is the effective time constant and 


k = tee — M. (27) 
If eq. (27) is substituted into eq. (26), then 
Z(8) FS Lac(s + m/Tact)/(s + 1/7). (28) 


REFERENCES 


1. A. H. Inglis and W. L. Tuffnell, “An Improved Telephone Set,”’ B.S.T.J., 30, No. 
2 (April 1951), pp. 239-270. 

. IEEE Standard 269-1971, “IEEE Standard Method for Measuring Transmission 
Performance of Telephone Sets.’’ 

. F. S. Goucher, “The Carbon Microphone: An Account of Some Researches 
Bearing on Its Action,’’ B.S.T.J., 13, No. 2 (April 1934), pp. 163-194. 

. H. Mol, “Theorie van de Koolmicrofoon,’’ PTT-Bedriff, 3, No. 4 (June 1951), 
pp. 128-134. 

. R. Holm, Electric Contacts—Theory and Applications, 4th ed., New York: Springer- 
Verlag, 1967. 

. C. A. Fritsch, “Téffects Associated With the Thermal Response of the T1 Tele- 
phone Transmitter,’ B.S.T.J., 47, No. 8 (October 1968), pp. 1615-1636. 

. R. Joscheck, ‘‘Electrishe and Mechanishe Eigenschaften des Kohlengrieszes von 
oe Wiss. Veroff, a.d. Siemens-Werken, 16, No. 1 (1937), pp. 

D-119. 

. M. C. Huffstutler Jr., and B. T. Kerns, unpublished work. 

. H. W. Bryant, ‘Comparable Tests on Linear- and Carbon-Type Microphones,” 
J. Acoust. Soc. Amer., 3, No. 3 (March 1973), pp. 695-698. 


N Oo oO FP WD WW 


Noles) 


1318 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1975 


Copyright © 1975 American Telephone and Telegraph Company 
Tuer Beutu SystEM TECHNICAL JOURNAL 
Vol. 54, No. 7, September 1975 
Printed in U.S.A. 


Some Far-Field Studies of an Offset Launcher 


By M. J. GANS and R. A. SEMPLAK 
(Manuscript received January 24, 1975) 


An offset paraboloidal reflector illuminated by a balanced feed horn 
constitutes an efficient launcher for coupling microwaves into quasi- 
optical beams. Measurements on a launcher with low blockage show low 
cross polarization. The amplitude, phase, and polarization characteristics 
are predicted by two gaussian beam modes, and the resulting formulas 
are found to agree well with measurements at 19 and 28 GHz. For example, 
with increasing offset angles, the ratio of the maximum cross-polarized 
signal in the radiation pattern to the on-axis co-polarized signal is ob- 
served to vary from —44 to —87 dB, within 1 dB of the predicted variation. 


l. INTRODUCTION 


At millimeter wavelengths, normal waveguide losses become too 
large in many applications; for example, long lengths of waveguide 
are required in satellite earth stations between the transceiver and 
the reflector antenna focus. To reduce these losses, one may use quasi- 
optical beams that employ reflectors or lenses for refocussing at 
various intervals, thereby confining the beam within a geometric 
tube with no (lossy) guiding walls. To couple the circuit components 
to these beams, it is desirable to provide a beam launcher that has 
quasi-gaussian amplitude over the aperture, low loss, good polarization 
purity, and high return loss. 

Offset reflectors provide high return loss, i.e., they are well matched, 
because the radiation field of the illuminated aperture bypasses and, 
therefore, does not reenter the feed horn.! If the reflector is made large 
so that the level of the feed-horn illumination at the edge of the 
reflector is low, spillover and diffraction losses are small. Also, since 
the feed horn does not block the aperture of the reflector, blockage 
losses are negligible, and the radiation patterns are unaffected. But 
cross polarization can be serious with offset reflectors, as pointed out 
in Ref. 2; however, by choosing a small offset angle the cross polariza- 
tion can be held to acceptable levels throughout the beam. 

Here, the far-field properties of an offset paraboloidal reflector 
(Fig. 1) are investigated. Section II describes the radiation charac- 
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Fig. 1—Offset launcher. The 76-cm (30-in.) diameter reflector is a numerically 
machined section of a paraboloid. Reflector focal length is 115.7 cm and it is fed, in 
this instance, by a 28.5-GHz dual-mode horn with a polarizer. 


teristics of the dual-mode feed horns, the experimental setups, and 
the measurements of the far-field properties of the complete launcher. 
Section III provides the theoretical formulas showing that the radiation 
performance of offset launchers can be well characterized in terms of 
gaussian modes. 

Specific applications of this type of launcher are in feeding Casse- 
grainian antennas of the type discussed in Ref. 1, and in launching 
and collecting beams on Hertzian cable transmission lines.’ 


il. MEASUREMENTS 
2.1 Dual-mode feed horns 


The dual-mode feed horn designed to feed the offset launcher of 
Fig. 1 is shown in Fig. 2a. The input section generates the TEi; and 
TMi modes in a circular waveguide by means of a conical step in the 
waveguide.‘ This section slides in the input waveguide so that the 
length from the step to the horn aperture, the “‘drift space,’’ may be 
adjusted to co-phase the TE,; and TM, modes to provide zero current 
at the edge of the aperture (minimizing the side lobes and symmetrizing 
the pattern). The small horn taper angle of 7.121 degrees was chosen 
to prevent disturbing the TE,; and TMy modes, and to provide a 
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FREQUENCY, F DIAMETER, D LENGTH, L 
IN GHz IN cm IN cm 
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Fig. 2—(a) Cross-sectional view of the dual-mode feed-horn design used for the 
offset launcher of Fig. 1. (b) Profile of measuring range. (c) Launcher schematic. 


small aperture phase error [e = (22/\)(D?/8L) = 7/2 radians | without 
making the horn too long. 

The azimuth radiation patterns of 28.5-GHz dual-mode feed horns 
for horizontal, vertical, and 45-degree polarizations, along with the 
associated cross polarizations for each case, as measured in an anechoic 
chamber, are shown in Fig. 3. Figure 3a shows the results obtained 
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Fig. 3—The 28.5-GHz dual-mode feed-horn radiation patterns for the principal 
and 45-degree polarizations. Included are the corresponding cross polarizations: (a) 
without polarizer on dual-mode feed; (b) with polarizer on dual-mode feed. 


without a polarizer on a horn. Note that the cross polarization for the 
45-degree polarization condition has peaks of about —28 dB at angles 
of about +8 degrees. (This behavior of dual-mode horns is predictable ; 
superior cross-polarization performance is obtainable from hybrid- 
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Fig. 3 (continued). 


mode horns.) One can also see that, even though the quarter-wave- 
length aperture error fills in the first null in the feed-horn pattern, a 
shoulder appears at about —20 dB relative to the on-axis value for 
the case of horizontal polarization. At least 30 degrees of the feed 
pattern illuminates the reflector of the launcher (resulting in about 
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a 20-dB taper at the edge of the reflector in the principal polarization). 
Thus, the cross-polarization maxima in Fig. 3a also illuminate the 
reflector. However, as shown in Fig. 3b, when a grid polarizer is used 
on the dual-mode feed horn, the cross polarization for the 45-degree 
condition is reduced to a very acceptable —45-dB level. Also, with a 
polarizer on the horn, cross polarization for polarization in the principal 
planes is essentially nonexistent. Since one of the main purposes of the 
experiment is to measure the cross polarization generated by the 
offset reflector, as discussed in the next section, the feed horn per se 
must therefore be devoid of cross-polarized components. For that 
reason, the feed horn was equipped with a polarizer (patterns of Fig. 
3b) for all ensuing measurements. 

The dot-dash curves in Fig. 3b show that the feed-horn pattern 
is well approximated down to about —20 dB by a gaussian beam 
[eqs. (8) and (9) ]. The 10-dB half angle, 6., used in the gaussian beam 
approximation, is 9 degrees. 

The radiation patterns for the 19-GHz dual-mode feed horn are 
essentially the same as in Fig. 3. Significant cross-polarization levels 
were also observed at this frequency, but, using the polarizer, the cross 
polarization is reduced to a very acceptable level (S — 50 dB). 


2.2 Antenna measuring range 


A profile of the antenna range used for measuring the offset launcher 
is shown in Fig. 2b. To determine the cross-polarization characteristics 
at the range, a gently tapered pyramidal horn with a 15- by 15-cm 
aperture equipped with a wire grid polarizer (to eliminate any cross 
polarization caused by the horn itself) was used as the source. A 
standard gain horn with a similar polarizer was used as the receiver 
on the antenna azimuth-elevation positioner. These measurements 
show the cross polarization introduced by the range to be very small; 
in the range of interest, i.e., within +3 degrees of the axial direction, 
it is of the order —47 dB. 


2.3 Offset launcher measurements 


Figure 1 is a photograph of one of the offset reflectors, along with its 
supporting structure. The reflector is illuminated by a 28.5-GHz 
dual-mode feed fitted with an etched grid polarizer as discussed in 
Section 2.1. 

At 28.5 GHz, the far-field patterns in the principal and 45-degree 
polarizations, along with the associated cross polarizations for the 
offset launcher fed by a dual-mode horn with polarizer, are as shown 
in Fig. 4. Note the shoulders rather than sidelobe structure for the 
co-polar patterns. Although the shoulders in the launcher pattern are 
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Fig. 4—Far-field radiation patterns of the principal and 45-degree planes for the 
offset launcher at 28.5 GHz with polarizer on horn. Theoretical (gaussian-beam ) 
calculations are also included. 


also apparent in the feed-horn patterns (compare the vertical polari- 
zation patterns of Figs. 3b and 4), this is not always the case. For 
example, the launcher pattern for horizontal polarization exhibits 
prominent shoulders (Fig. 4), whereas the corresponding feed-horn 
pattern (Fig. 3b) does not. Furthermore, the shoulders in the launcher 
patterns are at about the —24 dB level, while those in the feed-horn 
pattern are at about —20 dB. Most likely, the shoulders in the launcher 
patterns are due to the phase errors in the illumination which are 
caused by the finite taper length of the conical feed horn. 
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At 19.04 GHz, the co-polar and cross-polarization radiation patterns 
for the offset launcher are as shown in Fig. 5. Eixcellent symmetry is 
observed in the cross-polarization patterns even at this —40-dB 
level. Both Figs. 4 and 5 include the theoretical curves discussed in 
the next section. 

One can see from the configuration of Fig. 1 that there is a 
possibility of a small amount of blockage by the feed and its mount 
(and subsequent cross-polarization effects) when the launcher is 
scanned upward in elevation. To examine this, a set of azimuth scans 
for various elevation settings was made at both 19.05 GHz and 28.5 
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Fig. 5—Far-field radiation patterns of the principal and 45-degree planes for the 
offset launcher at 19.05 GHz. Dual-mode feed is equipped with polarizer on horn. 
Theoretical (gaussian-beam) calculations are also included. 
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GHz; the increase in cross polarization due to blockage was found 
to be negligible in both cases. 

One may with justification raise the question as to why the cross 
polarization of the offset launcher, shown in Figs. 4 and 5, evidences 
values of the order of 38 dB even with a feed that has negligible cross 
polarization. There is an inherent depolarization introduced by an 
offset reflection surface,? which increases if 60, the angle between the 
feed axis and the reflector axis, is increased. Figure 6 shows the experi- 
mental results obtained by varying 6> between 12 and 26 degrees; 
the ratio of the maximum cross-polarized signal in the radiation 
pattern to the on-axis co-polarized signal correspondingly varies from 
—44 to —37 dB at both frequencies. Note that the ordinate of Fig. 6 
is the average of the peaks of cross polarization obtained from an 
azimuth scan of the launcher; they should not be misinterpreted as 
on-axis values which, of course, are much lower. In the following 
section, we show that calculations based on gaussian-mode theory 
provide good agreement with the measured data; the theoretical 
result is shown by a solid curve in Fig. 6. 


Ill. THEORY 
3.1 Cross polarization in the aperture 


An approximate method for computing the cross polarization due 
to the offset angle 00 consists of applying geometric optics to compute 


© “~GAUSSIAN MODE 
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Fig. 6—Plot of cross-polarization peaks at 19.05 and 28.5 GHz introduced by the 
reflector itself as a function of offset angle, @. Data are obtained by seanning the 
launcher in azimuth. Gaussian-mode theory is shown by the solid line. Incident 
polarization is horizontal. 
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the reflected field in the aperture from that radiated by the feed horn. 
The aperture field is then decomposed into two gaussian-beam modes 
to predict the far field of the offset launcher; this is a logical procedure 
because dual-mode horns produce an illumination that is approxi- 
mated by a gaussian beam (see Fig. 3b). 

The geometry of the offset launcher is shown in Fig. 7a. The dual- 
mode feed horn approximately provides “balanced feed’’ polarization‘ 
with respect to the coordinates (2’, y’, and 2’) aligned with the horn 
axis which is tilted at angle 4. from the reflector axis. The field radiated 
by the feed horn is given by 
(0', > ) 


E, = (6 sind’ + ¢’ cos g') L028) (1) 


where f’(6’, ¢’) is an arbitrary function of 6’ and ¢’, and 7’, 6’, and ¢’ 
are the usual spherical coordinates associated with the feed (see 
Fig. 7b). The caret indicates a unit vector. In eq. (1), the expressions 
corresponding to ‘vertical’ polarization are used; identical results 
are obtained for “horizontal” polarization. [‘‘Vertical’’ and “‘hori- 
zontal’ are used in the sense that the polarization of the field in the 
aperture of an axisymmetric paraboloidal reflector coaxial with the 
feed-horn axis would be vertical or horizontal when the feed-horn 
polarization is as given in eq. (1). ] 

The axis of the paraboloidal reflector shown in Fig. 7a is co-linear 
with the 2 axis. An aperture field with no cross polarization would, 
therefore, result if the feed illumination were given by 

E, = (6sing + ¢cos ¢) 16,9) 6) (2) 
where r, 6, and ¢ are the usual spherical coordinates associated with 
the x, y, 2 coordinates of Fig. 7a, i.e., the coordinates of a feed horn 
whose axis is aligned with the reflector axis. Theoretically, it is possible 
to hypothesize an asymmetric ‘‘balanced’”’ feed whose axis is aligned 
with the reflector axis, giving the polarization of eq. (2), but whose 
amplitude distribution is offset to illuminate the reflector as would 
the amplitude distribution of a tilted symmetric ‘‘balanced’’ feed. 
Simple means (excluding multiple reflectors, etc.) are not known for 
the construction of such an asymmetric “balanced” feed. Therefore, 
In applications, one must approximate an asymmetric “balanced”’ 
feed horn with a tilted ‘balanced’ feed horn; the cross polarization 
thereby introduced is calculated below. 

Using geometric optics, we assert that if the polarization of a ray 
incident on the reflector from the feed is rotated by a given angle 


1328 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1975 







OFFSET 
~ -PARABOLOIDAL 
REFLECTOR 





AXIS_OF 


FEED HORN 





Fig. 7—(a) Launcher geometry. (b) Spherical coordinates. (c) Spherical triangle. 


around the ray vector, then the polarization of the field of the corre- 
sponding ray in the reflector aperture is rotated by that same angle. 
By geometric optics, the intensity of the field in the aperture along a 
ray is the same as the intensity of the field incident on the reflector 
from the feed at its focus. Thus, the cross polarization in the aperture 
relative to the peak in-line polarization in the aperture can be com- 
puted by projecting the field of a tilted ‘“‘balanced”’ feed horn, incident 
at any point on the reflector, onto the cross-polarized asymmetric 
“balanced”’ field at the same point and dividing by the peak in-line 
field of the tilted “balanced” feed horn; 1.e., 





E,-[6(— cos ¢) + ¢sin o]. (3) 


C(O’, ¢’) = E;-[6 sin @ + 6 c0s $]| peak 
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Performing the scalar products indicated in eq. (3) yields 


C(O’, 6’) = sin (@ — ¢’ + nies) 


| cos (@ = ¢’ + ie) ¢’) 





call (4) 


where 7 is the angle between the primed and unprimed spherical 
angle coordinates shown in Fig. 7c. The angle, ¢’ — ¢ — 7, is equal 
to the area, A (called spherical excess), of the spherical triangle shown 
in Fig. 7c and is related to the offset angle, 40, and the primed (feed 
horn) coordinates by the formula® 





— @ 
cos ( 
oi nbeats Se 2 oy m 
A=? @— 7 = 2 arctan (aewyn(3 +4) 
2 


The feed-horn amplitude pattern, f’, is approximately uniform in 
¢’ and maximum on axis (as for a corrugated or a dual-mode feed 
horn), and r’, the distance from the feed to the reflector, is relatively 
constant over the region in which the cross-polarized field is a signifi- 
cant fraction of the peak in-line field. In this case, the cross-polarization 
amplitude can be approximated by 


ee ot 1) 


/ 0’) 
C(6’, ¢’ — sin AL) ( 6 
6, o) = — sin AGT. (6) 
For a given 6’, the ¢’ which maximizes A and the cross polarization is 
dm = arcsin (tan i“ tan 5 ): (7) 


For offset angles, 60, less than or equal to 90 degrees and 6’ less than 
30 degrees, which covers the case of interest, dn is less than 16 degrees. 
This leads us to approximate ¢, by 0 degrees, which results in a 
particularly simple yet accurate formula for the peak cross-polarization 
amplitude ratio, C, as will be shown. C is the ratio of the maximum 
cross-polarized amplitude to the maximum in-line polarized amplitude. 

For dual-mode and corrugated horns, the pattern is approximated 
by a gaussian beam,’ 





f'(') = ea”, (8) 
where a is related to the 10-dB half-angle beamwidth, 6, by 
In 10 
a= 3. (9) 
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Approximating ¢, by 0 degrees, we have the following relation® for a 
right spherical triangle (Fig. 7c): 


sin 69 sin 6” 





a= 1 + cos 8 cos 6” (10) 
Since we are interested in small 6’, 
‘ _ 6’ sin Ao 
sin A= 1 + cos by } (11) 
and eq. (6) becomes 
, ee sin Ao 1a '2 
C6 ) dm) 1 + cos Oo . (12) 


By differentiating eq. (12) to locate the angle which maximizes 
C'(6’, bm), We find 


= 55 = 6,/vIn 10, (13) 


and 


Ho 
0, tan Dy) 


ve ln 10. 


As mentioned above, C is the ratio of the maximum cross-polarized 
amplitude to the maximum in-line polarized amplitude. The subscript 
a indicates that the ratio is of the maxima found in the reflector 
aperture. If we denote by 67 the half-angle of the gaussian beam 
approximation to the feed-horn pattern, where the power is 7 dB 
below that on axis, then eq. (14) becomes 


. a | 10 
C, = — 6p tan 5 eT In 10 (15) 


By comparing eq. (14) with the exact formula for gaussian beams, 
eqs. (6), (7), and (8), eq. (14) is found accurate to within 0.1 dB for 
all offset angles, #0, less than or equal to 90 degrees and all 10-dB 
half-beamwidth angles, 6@,, less than or equal to 45 degrees. The 
maximum value of C in the exact formula was found by trial and error 
with the aid of a calculator. 

The above calculations give the cross polarization in the reflector 
aperture; as will be shown in the next section, eq. (14) is also a good 
approximation for the far field in most cases. 


C,+=- (14) 


3.2 Two-mode approximation to the aperture field 


After reflection from the offset reflector, the gaussian beam from 
the horn is converted into two gaussian beam modes in the aperture: 
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a fundamental mode with the in-line polarization (denoted Eoo) and 

a higher-order gaussian beam mode which includes the cross polariza- 

tion (denoted Epi). Depending on the polarization of the feed horn, 

the polarization of the fundamental and higher-order modes vary 

as shown in Fig. 8. For an arbitrary balanced feed polarization, a 

superposition of the two polarization cases shown in Fig. 8 can be made. 
The expressions for the gaussian beam modes are? 


(1) Fundamental mode: 


2 
Eoo = (Hoot + Voo gy exp {- dke" — - 
Wio 
Del! kp? 
+ j arctan (=, ) _~ a || , (16) 
(iz) Higher-order mode: 
Eo. = |Vo(é" cosa — 9’ sina) — Hou (é”’ sna + Y” cosa) own 


Peds kp 
9) cutee : 
ry De ae il2 arctan (= )- IR. (17) 


nt + an - \\/) 


“exp |- gk2"" — 


Eo0 + Eo, — EToT 
FUNDAMENTAL HIGHER ORDER TOTAL 
(a) 

—— ie Men an 

oO > i a 
ae | = eae 

—_—_ ee 

—_ a. —-P 

Ea0 + Foy — Eror 

FUNDAMENTAL HIGHER ORDER TOTAL 


(b) 


Fig. 8—Two-mode decomposition of aperture field (polarization looking in positive 
z direction, 1.e., looking at the reflector). (a) Feed horn vertically polarized. (b) Feed 
horn horizontally polarized. 
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Hoo, Voo and Hoi, Voi are the phasor coefficients for horizontally and 
vertically polarized feeds. The subscripts refer to the standard TEM oo 
and TEM, mode notations of Ref. 8; p, a, and 2’’ are cylindrical 
coordinates, with 2’’ denoting the distance along the beam axis from 
the beam waist. At the beam waist, the radius of curvature of the 
phase front, R, is infinite, and the field varies with increasing distance, 
p, from the axis at a rate determined by w. For the fundamental mode, 
the field decreases to the 1/e value at p = Woo. For the higher-order 
mode the field is maximum at p = Woi/v2 and decreases to V2/e of 
that value at p=Wo1. Away from the beam waist, 2’’ ~ 0, the field 
amplitude varies with p at a rate determined by w instead of w, and 
the phase front has a finite radius of curvature, R. w and FR are given 
by® 








Dat N2 
w=w 1+ (4), (18) 
and 
. kw? \?) 
k=2 [1+ ($e) | (19) 


Both modes have a characteristic exponential attenuation with 
distance from axis, e~°’/”’, and a spherical wave front near the axis at 
constant 2’’, denoted by the term, e~-?/2®, Passing through a beam 
waist, the on-axis phase advances by z for the fundamental mode and 
2m for the higher-order mode (relative to the plane-wave retardation, 
e—Jk’’). Thus, if the cross polarization and the in-line polarization 
are in phase at the beam waist (normally near the reflector aperture), 
they will be in phase quadrature at large distances from the beam 
waist (the far field). This phase quadrature relation gives rise to a 
beam shift with circular polarization as described in Ref. 2. 

The choice of eq. (17) as the appropriate higher-order mode is based 
on its ability to approximate simultaneously both the cross polariza- 
tion and the “space” taper (amplitude asymmetry from top to bottom 
of dish) of offset reflectors. 

The in-line and cross-polarized fields in the aperture of the offset 
launcher of Fig. 1 were computed exactly by means of eq. (3); the 
resulting field-amplitude contours are shown in Fig. 9. In the aperture 
plane, the in-line and cross-polarized fields are in phase. Thus, the 
corresponding gaussian beam modes have their beam waists at the 
aperture and are in phase. This implies that the total field is linearly 
polarized everywhere in the aperture, and the direction of polarization 
varies in a manner determined by the ratio of the in-line and cross- 
polarized fields. 


OFFSET LAUNCHER = 1333 


KEY: SOLID CURVES 
DASHED CURVES 


GEOMETRIC OPTICS WITH FEED-HORN PATTERN OF eg. (8) 
GAUSSIAN — BEAM MODES 


10cm 





7 “EDGE OF 


Y REFLECTOR 
a 
oe 
er 
IN-LINE POLARIZATION Se CROSS — POLARIZATION 
CONTOURS —- —— CONTOURS 
F=115.7cm 


76-cm DIAMETER APERTURE 


Tig. 9—Amplitude contours for the example reflector shown in Fig. 1. 


The gaussian beam fields required to match the exactly computed ap- 
erture fields are found by first choosing a fundamental mode centered 
on the aperture with beam waist radius woo such that it decreases 
8.686 dB in power at the same radius as does the computed in-line 
polarized field, both being normalized to unit amplitude on axis. 

With these criteria and the approximation that r’ is nearly constant 
over the reflector, 


== Fseot(Q), (20) 


the parameters of the fundamental mode are determined as follows. 
The radius, p1;e, at which the field drops to 1/e times its on-axis value 
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is found in the direction a = 0. (Note that p;;, is largerin thea = —7/2 
direction and smaller in the a = 7/2 direction due to the space taper 
in these directions. Thus, p1;-, as determined from the a = 0 direction, 
approximates the average p1;- over all directions.) 


Woo A prye = 270 Sin ($61/c), (21) 


where, from the gaussian approximation for the feed-horn pattern, 


2 
6. = O. mene 29 
1/ In 10 ( ) 


For unit amplitude on-axis, eq. (16) requires (we restrict our dis- 
cussion here to the vertically polarized case; the horizontally polarized 
case yields identical expressions) 


Voo0 = 1. (23) 


The higher-order mode parameters are found from the cross-polari- 
zation characteristics. Since the cross polarization is maximum at 6, 
we have, from eq. (13), 

' : i] 
Wo = 2v2ro sin (36,,) = 2v2r sin ( : ): 24 
DS arn Aa gee AD Gina i 
By comparing eqs. (24) and (21), it is seen that, for small 6., Wo1 = Woo. 
The amplitude of the higher-order gaussian beam mode is given 


by the maximum cross-polarization amplitude ratio (C,). From eqs. 
(14) and (17): 





The phase of the higher-order mode follows from the fact that at the 
beam waist the in-line and cross-polarized fields are in phase. 

Using eqs. (21) through (25), the gaussian beam mode approxi- 
mations to the aperture fields, plotted as dashed contours in Fig. 9, 
compare favorably with those obtained by geometrical optics (solid 
contours). 


3.3 The far fields 
The parameters of the gaussian beam modes being thus determined, 
it is possible to compute the in-line and cross-polarized fields at any 
position in the main beam of the field radiated from the reflector. 
The far field is of particular interest, 
kw? 


2" >>, (26) 
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allowing eqs. (18) and (19) to be approximated by 


tt 
w= a and R= 2”. (27) 


Also, if the beamwidth is small, 


=V@P FR Haw + SE; HF 


22 5,773 2 


(28) 


From eqs. (16) and (17), the in-line field is therefore given by 


jke- srs 67 
Ly = Z V vow exp { — a 
; 2 

21s Yoo 


2 
— jv2V owes ue (sin a) exp (- we. ) | (29) 
Yo Vor 
and the cross-polarized field by 





_ kee jkr 


62 
Ey = V2V ewes — -, (eos a) exp (~ i) : (30) 
o1 





2r; 
for the ‘‘vertically” polarized feed. We have defined 





2 
Yoo A EWoo and You A ’ (31) 


the angular beam radii in the far field. 
The C,; from eqs. (14), (23), (25), (29), and (80), where the subscript, 
f, indicates the far-field maximum cross-polarization amplitude 


ratio, 1s 
. Vou (= Wo1 , (= ) 
- =c,( 2"). 2 
1 ’ Voove Woo Woo oe 


Thus, the far-field C is 7(wWoi/Woo)? times that in the aperture and 
occurs in the azimuth plane at an angle 





vou = ~ radians. (33) 
01 


As mentioned above, for small feed beamwidths, 4, 
Wo1 — Woo, (34) 


so the far-field C is approximately equal to that of the aperture, and 
the peak cross-polarization lobe occurs at approximately the —4.34- 
dB level of the main beam. 

A comparison of the experimental and theoretical far-field patterns 
for in-line and cross-polarized fields from the offset launcher of Fig. 1 
at 28.56 GHz and 19.04 GHz are shown in Figs. 4 and 5. The main 
discrepancy between theory and experiment are the shoulders on the 
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sides of the experimental pattern at a level of about —24 dB not 
present on the sides of the theoretical main beam. As mentioned in 
Section 2.3, the shoulders are probably due to phase error in the horn 
aperture which is not in the gaussian feed pattern assumed in the 
theory. 


3.4 Truncation effects 

The effect of truncating the aperture of the launcher at various 
circular contours of the fundamental mode can be computed. Let the 
radius of the aperture at the truncation be c. Then the taper at the 


truncation is 
20 Ee 
P= (0) ea) (8) 


The radiation integrals in the 2’’, z’’ plane (where the cross-polarization 
is largest) for the fundamental mode and higher-order mode are 


Qn 


c 2 
Kou = | pdp ‘ da V 9 exp |- (£) — jkp cosa sin |, (36) 


Woo 
and 
_ c Qr v2p 
Koy i I pdp { daVorT 
2 
‘exp |- (= ) — jkp cosa sin | cosa cos 6s. (37) 
Wo1 
From eq. (36), the in-line polarization far-field on-axis is 
2 
E oo, (84 = 0) => rVoow3o| 1 — exp (- <)| < (38) 
00 


The cross-polarized far-field pattern from (37) is, letting cos 6; = 1, 


. 2 
. JtV Wor 


Koy = “Ae. I(6;), (39) 
where 
c/wo1 . 
1(6;) & 2V2e iF Pe-"°J (twos sin 6,)dt, (40) 
0 


and J1(2) is the first-order Bessel function of the first kind. From eqs. 
(84) and (85), choosing the taper, 7, determines /(@;). Numerically 
integrating eq. (40), the peak value, Jp, of J(6,) and the location, 
6;,, at which J(@;) is maximum were determined for various tapers, 7’. 
The resulting far field C; relative to that of an infinite aperture, 


C I 
— 20 logio (s5) = — 20 logio ( we) ? (41) 
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is plotted in Fig. 10 along with the angular position, 6;,, of the peak 
in the cross-polarization lobe relative to that, @sp., of an infinite 
aperture. Notice that the effect of truncation is to reduce the cross- 
polarized signal relative to the co-polarized signal (higher C,;) and to 
move the cross-polarization lobes out to larger angles off-axis. Further- 
more, these truncation effects depend only on the truncation taper 
and are essentially independent of the offset geometry. 

Figure 11 shows Cy as determined from the gaussian beam formulas 
for a wide range of offset reflector geometries indicating an infinite 
aperture (no truncation) and the aperture truncated at a 10-dB taper. 
With the infinite aperture, the C; from the gaussian beam theory are 
within 0.2 dB of the cross polarization in the aperture, C, (not shown), 
for all geometries on the figure. When the aperture is truncated at the 
10-dB level, Cy is, from Fig. 10, 2.2 dB smaller than that for the 
infinite aperture, as is also seen from the dashed curve in Tig. 11. 
These dashed curves agree with cross polarization obtained by nu- 
merically computing? the radiation integral over the aperture field 
found by geometrical optics projection of the radiation pattern of a 
dual-mode feed horn. 

Using eqs. (14) and (32), the maximum cross-polarization levels 
as a function of offset angle were computed for the precise geometry 
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Fig. 10—Effect of aperture truncation on cross-polarization pattern. 
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INFINITE APERTURE (NO TRUNCATION) 
APERTURE TRUNCATED AT 10-dB TAPER 


10—dB 
HALF 
BEAMWIDTH 


20 LOGi9 (Cy) IN dB 





0 15 30 45 60 75 90 
OFFSET ANGLE IN DEGREES 


Fig. 11—Maximum cross-polarization amplitude ratio for a range of offset. geometries. 


of the offset launcher of Fig. 1. The factor (Wo1/Woo)? which converts 
aperture cross polarization to far-field cross polarization is only 1.0004 
(i.e., 0.008 dB). The 10-dB half angle, 6., of the feed-horn gaussian- 
beam approximation (discussed in Section 2.1) is 9 degrees. The 
calculated cross polarization is compared with the measured cross 
polarization as a function of offset angle in Fig. 6. The theory appears 
to be in good agreement with the measurements at both frequencies. 


IV. CONCLUSIONS 


It has been demonstrated that an offset launcher can provide low 
cross polarization and a low-sidelobe symmetrical beam when fed 
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with a suitable horn at a small enough offset angle. Simple formulas 
for the far-field performance of the launcher are derived in terms of 
two gaussian modes; comparison with measurements at 19 and 28 
GHz shows good agreement. The maximum cross-polarization ampli- 
tude ratio is found to change little from aperture to far field. Offset 
reflector geometries have also been useful for multiple-beam 
applications.?: 
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Upper Bound on Error Probability for Detection 
With Unbounded Intersymbol Interference 


By A. D. WYNER 
(Manuscript received February 27, 1975) 


Forney’s asymptotic upper bound for per-bit error probability in the 
detection of pulse-amplitude-modulated digital data in the presence of 
additive white gaussian noise was obtained for the case where the dura- 
tion of the intersymbol interference is bounded. In this paper, we show the 
validity of Forney’s bound under much weaker assumptions that allow 
unbounded intersymbol interference. 


]. INTRODUCTION 


We consider the situation where a data sequence do, ---, @v-10f +1’s 
is transmitted via pulse amplitude modulation as }(?=o h(t — kT)ax 
and received in the presence of additive white gaussian noise with one- 
sided spectral density o?. In a recent series of papers, Forney,! 
Foschini,? and Mazo® developed an asymptotic (as o?—0) upper 
bound on the error probability per data bit P.: 


- ft +oal, (1 


where d(h) is the minimum £, distance between distinct modulated 
pulse sequences. This bound holds under the strong assumption that 
the pulse A(¢) is supported on finite interval. 

In this paper, we show that (1) is valid for a considerably wider class 
of h(t). Roughly speaking, our assumptions are little more than that 
h(t) is in £L1(— ©, ~) and L2(— ~, ~), and that H(f), the Fourier 
transform of h(t), does not vanish on an interval. The precise conditions 
on A(t) under which (1) holds are given below. In particular, (1) is 
valid when H(f) is a rational function. 

In Section II we give a precise statement of our results, and the proof 
follows in Section IIT. 


P, S exp 





Il. FORMAL STATEMENT OF PROBLEM AND RESULTS 


In this section, we give a precise statement of the problem and the 
results that were stated informally in Section I. 

We begin with some definitions. We denote N vectors by boldface 
superscripted letters, and components by subscripted letters, e.g. 
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u* = (uo, +++, Uv—1). When the dimension N is clear from the context, 
we omit the superscript. Define the sets @y, @ix, Gy, by 
@v = {wWiuj= +105 785 N —-1}, 
Qt, = {uY:u = +1uj=+1,j3# k}, (2) 
@vx = {UX up = —lu= 41,74 k}. 


I 


ll 


Of course, @vy = Qi, U @y,. Again, when N is clear from the context, 
we write @y = @, @ix = Qi, Qh: = Qi - 

Next, let f(), yg), and —»<t<o be real-valued measurable 
functions. The inner product of f and g is denoted by 


(9) =f Og@at, (3a) 


and the norm of f is 
l= i = (Cf Poa) 


For a vector u* € @y, and f(t), -»<t <, a real-valued function, 
let the function h#u = s be defined by 


1 
2 


(3b) 


si) = D3 f(t — kT)ux, 


where 7’ > 0 is a fixed parameter. 

We are concerned with the following modulation scheme. Let 
aN = (do, --+, @v-1) € @y denote the data to be transmitted. Assume 
that all the 2% vectors in @y are equally likely. The transmitted signal 
is the function hxa%, where the pulse h(¢) is a fixed function for which 
\|h||< co. The received signal is 


y(t) = (ha) (t) + 2), —a<lti<coa, (4) 


where z(t) is a sample from a white gaussian noise process with zero 
mean and one-sided spectral density o”. 

The decoder associates with the received signal y, a vector 
D(y) = a” © @y. Corresponding to a given decoder function D, let 
the bit error probability be 


1 N=1 . 
Pin(D) = FH a Pr {d, 4 ay}. (B) 
k=0 
Also, define the optimum error probability 


Pty = Pty(h, 6”). = inf Pey(D). (6) 
D 


We are concerned here with the asymptotics of Pey(h, 0”), as «2 — 0, 
i.e., as the signal-to-noise ratio approaches infinity. Accordingly, define 


Ey(h) = — lim inf 6? log Pty(h, 0), (7a) 
o?—0 
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so that, as o? > 0, 





Bl) Tat ety hs (7b) 


og 


Pin (h, a”) S exp )— 


Next, consider a particular decoder that is of special interest here— 
the maximum-likelihood decoder, denoted D,. In the present problem, 
Diy) can be taken to be that i © @w such that for all u € Gy, 
u ~ G, 

Ily1 — hxtil]<|lys — haul], (8) 
where yi is the projection of y onto the subspace of £.(— ~~, «) 
spanned by the signals h#u, u © @y. With probability 1, (8) will be 
satisfied for some i € @y. 

Now, subject to the condition that A(¢) has finite support, i.e., there 
exists a to > O such that 


h(t) = 0, for |t| > to, (9) 


Forney,! Foschini,? and Mazo’ have shown that E(h) = d?(h)/4, where 
the “minimum distance’’ d(h) is defined by 


d(h) = lim inf min |/h*xu — hav]. (10) 
Nos ane 


Thus, as o? — 0, 


Pey(h, 0?) S exp }— 





d?(h) 

pe [1 + 0) (11) 
Inequality (11) is established by showing that the error probability 

for the maximum likelihood decoder, P.w(Dz), is overbounded by the 

right member of (11). This is done by writing (this is not as difficult as 

it looks) 


1x ' 
Pex(Ds) = HD, 2, OP Pr (Day) © @xaly = how + 2) 
=1 u€QGru 
= FE ESO Pr | YU [Daly) = vi ly = how +2] 
NF weat veQGr 
1 
as 9-(N-1) 
a NY eae 
‘Prt Ue (lly: — Asul] 2 lly: — havi|}ly = heu + 2} 
ve Gk 
=5d SFO Pr U te he(v — w)) 
NS uc Qe veEQs 


2 3llh*(v — u) II} 


S px (h, 0°). (12) 
Relation (12) is valid for any h(t). Subject to condition (9), it is then 
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shown that, as o? > 0, 


vy (h, a”) = exp |- a) [1 =e oi(1) ] ’ (13) 


where 01(1) does not depend on N. Thus, since Pty < Pey(Dn), (11) 
holds. Further, the o(1) term in (11) does not depend on N. An 
interesting by-product of these results is that the performance indi- 
cated in (11) is achievable via the decoder D,. This decoder can be 
instrumented (using the Viterbi algorithm) with a complexity which 
remains bounded as N > «. 

We now drop the assumption that A(t) has finite support. Instead, 
we assume that /(t) satisfies the following conditions: 


(t) There exists a nonnegative £; function go(t), 1.e., 


| go(t)dt <<, 


such that 
Ih(t)| S got), —-e<t <a, (14) 
and such that go is monotone in |é]. 
(77) Let 
HGy= / . h(e®@Vdl, —o<f<x (15) 


be the Fourier transform of h(t). By (7), f|h(t)|dt <«, so that H(f) 
is well defined for all f. We assume that there exists a nonnegative £, 
function Gi(f) which is monotone in | f|, such that 


IA(NP Sa, —-#9<f<o. (16) 
(217) Let the “folded spectrum”’ of h be 
sn- % lw(t+%)), Osrsy an 








We show in Appendix A that S(f), 0 S f S$ 1/7, is finite and con- 
tinuous. We assume that S(f) > 0,0 S f S$ 1/T. Let 


m= min S(f)>0, (18) 


Osf $1/T 


where the existence of the minimum follows from the continuity of 
S(f) on the compact interval [0, 1/T']. 


Remarks: 

(1) Condition (2) is just slightly stronger than simply requiring h 
to be in £1(— ~, ~). Condition (14) forces h(t) to go to zero as |f| > 
in a “‘well-behaved”’ manner. Condition (77) imposes a similar condition 
on |H(f)|?. 

(2) For the very important special case where H(f) is a rational 
function, ie., H(f) = P(@2rf)/Q(27f), and P, Q are polynomials 
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with the degree of P < degree of Q, then conditions (2) and (72) are 
satisfied. Since H(f) has only a finite number of zeros, condition (772) 
is also satisfied. 

(3) Suppose that H(f) has no more than a countable number of 
zeros, but that S(f) = 0 for some f € [0, 1/7']. It is easy to see that 
some arbitrarily small change in 7 will cause S(f) to be strictly 
positive for all f € [0, 1/7]. Thus, condition (722) is not especially 
restrictive. 

We now state our main result, the proof of which is in Section III. 


Theorem 1: Let h satisfy conditions (2), (it), and (a1) above. Then, for 
all e > O, there exists a To = Tole) sufficiently large so that, for all 


7 > Fo, 


P.n(Dy,) < WwLh., o7(1 + e)*], 


hd), [t| = 7; 
Ogee (19) 
0, Falecar 


where 


A 


is the truncated version of h(t). The quantity ro does not depend on N. 
Since h, has finite support, we conclude from Theorem 1 and (13) 
that, for all « > 0 and - sufficiently large, 
Pew(h, e)s Pen (Dn,) 
(hz) 





S exp |- 421 + [1 + o2(1)]}, (20a) 
[where 02(1) is independent of N ] so that 
2 
By(h) = — lim inf o? log Pi(h, o#) 2 wey (20b) 
We show in Appendix B that 
d(h,) — adh), as T—> 0, (21) 
so that letting «> 0 and 7 © in (20b) yields 
¥3 
Ey(h) =" ), (22) 


We state this as 


Corollary 2: Let h satisfy conditions (2) to (21) above. Then, as o? — 0, 





h 
Pin(h, o Ss exp - ") Lhe + 03(1) ] 
where 03(1) 1s independent of N. 
We conclude this section with a remark concerning the relationship 


of the bound of Forney et al. (11) with the result of Corollary 2. We 


UNBOUNDED INTERSYMBOL INTERFERENCE 1345 


can rewrite (11) as 


Pen (h, a”) aa Ki (to, re | a es (23) 
and the bound of Corollary 2 as 
Piv(h, 02) S Ko(|lAlla, m, o2)e7@?? 1427, (24) 


Here, we made explicit the dependence of A, on the support interval 
to of h(t) [see (9)], and the dependence of Ky on ||h\|1, the £1 norm 
of h, and on m = min S(f). Both Ky and Ke increase in 1/o? slower 
than e?()/40?, But K1(to, o?) > as ty, and Ke(|\h||1, m, 0?) >, 
as ||h||; 0 or as m—> 0. Thus, although it might seem reasonable 
to assume that all A (é) satisfy (9) for some to, the bound of (23) depends 
on that t and becomes meaningless as tp >. Similarly, although 
it might be reasonable to assume for any h/(é) that ||h\|1 < © and 
m = min S(f) > 0, the bound of (24) depends on these quantities 
and also becomes meaningless as ||h||1;} > © or m — 0. Therefore, both 
bounds have their limitations; the new one, however, is considerably 
less limited. 


Ill. PROOF OF THEOREM 1 
_ Let A satisfy (2) to (wi). Let h,(f) be as defined in (19), and let 
h,(t) = h(t) — h-(0), ie., 
z = 0, |¢| = 7; Or 
iO = Ia, Ue > — 


T. 


Then, if the data sequence is u € @y, the received sequence is y 
= heu + 2 = h,xu + 2, where 


S=24+ hyxu. (26) 


Following the same steps as in (12), we obtain 


Pev(Di,) $2 OY Qe 
N ve at 


Pr U {(,hex(v — u)) = SllAsx(v — u)|]}, (27) 
ve: 


where 2 is given in (26). 
We will show that, for arbitrary « > 0, there exists a to = To(e, h) 
(ro independent of N), such that for r 2 70, the event 


{(2, hex(v — u)) 2 allher(v — u)||?} 


Cl (2, heey — w)) 2 SEO aay — w) Ph, 


(28) 


for all u € @f, v € @;. Substituting (28) into (27) yields, on com- 
parison with (12), 


Pen(Drn,) s VvLh-, o*(1 a e)* |, 
which is Theorem 1. It remains to establish (28). 
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Relation (28) will follow immediately when we show the existence 
of a ro(e, h) such that, for r 2 79 and all u, v, 


| i+u, hee(v — u))| S 5llhee(v — u) I. (29) 


If (29) holds, the event in the left member of (28) 


{ (2, hex(v — u)) 2 3||h-*(v — u)||?} 
= {(z,h.*(v — u)) & 3llh.*(v — u)||? — (h,4u, he*(v — u))} 





C]@ nel) = Eine — wie, Go) 
which is the right member of (28). Thus, it remains to establish (29). 

Let w = (Wo, +++, Wy-1) = v—u. The entries of w are 0, +2. 
Also set gq = 4-*u, andr = h,«w. Then 


lA 


| Give, hee(v —u))| = Harri Sf lal Ir@lae 
C sup lel f- ir@la. Gn 


IIA 


Consider 
20 co |N—-1 
/ [7 (t)|dt = ; ps h,é@ —kT)w: 


00 N-1 
SX [wel fo [be — kT) |dt < [all SO leat, (82) 


dt 





where ||h||1 = /|h(t)|dt <«, by condition (7). 

We obtain an upper bound on © |w;| as follows. Since w, = 0, +2, 
we have > |wz| = 4 > wz. Now, let A,(f) = (2. h-(He®™'dt be 
the Fourier transform of h,(t), and let 

ro) n 2 
s.in= = |HA(r+ g)) 


0Sfsa (33) 





be the corresponding folded spectrum. Then, from Parseval’s theorem, 


In|? 


Weawilt =f [20-2121 wne*! [2a 
—00 k 
1/7 
[SOE wre” Pap 
0 k 


[ inf snipe |S> wpe®t*7s | 2df 
O<f S1/T 0 
= Linf §-(/) J 2 WE = Linf S-(f) J > |wej-2. (34) 


IV 


Therefore, 
DY lwel SaliriPC inf SCA P. (35) 
B O<f <1/T 
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Combining (31), (82), and (35), we have 
lar) $C sup [q@|JC inf S.(f)} Sllalfalirii2 (86) 
—a<t<o O<f S1/T 


Now, we show in Appendix A that 
liminf[ inf S.,(f)|]2m>0, (37) 
Sf S1/T 


TO Oo<fs 


where m = min S(f) > 0 [see condition (777) ]. Further, using condi- 
tion (2) [particularly the monotonicity of go(t) ], we have 


GO| =| DAG bp SD ie - kD)| 
<> |k@—-kD|= ST [ke -kD| 
k=—o k:|t—kT| 2r 


= > golt -kT)s Dz [go(r + 7T) =f go(—T —j7T)] 
k:|t-kT| 27 j=0 
1 se 1 se 

<5 = go(at+= i: " _ga(t)dt +0, as ree, (38) 


Combining (386), (87), and (38), we obtain 


Mary] _ [Gieu, haw — wu)! 9 
InP [how — ul : 


as ro. This is equivalent to (29), so that the proof of Theorem 1 
is complete. 
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APPENDIX A 
The Folded Spectrum S(f) 


We first show that S(f) as given in (17) is always finite, i.e., the 
series in (17) converges for all f € [0,1/7]. From condition (12), 
using the monotonicity of Gi, 


(no, f) |n| 2no a(t - tI) = inn als = 7) 
EL atite)+ Eats), 


nS-n 


n+1 n 7 
e caf ‘a )+ a a(t) - ‘seer mane 


nS—no n2=no 





I 


lA 





+ af Gi(z)dtz—> 0, asm, (89) 
ng—-1/T 


so that the series in (17) converges. 


1348 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1975 


To establish the continuity of S(f), write 


A(t +H) + eo, O<f<1/T. (40) 





si= = 


=no 


For arbitrary 6,0 $f < 1/T, 


isin —su+ ol s|z [art a) | [e+e al] 
+[E(mo, f)| + [€(mo,f + 8)|. (41) 


Now since h(t) © £1(—”, ~), A(f) is continuous. To make the 
right member of ineq. (41) Se, first let no be sufficiently large so that 
the last two terms of the right member of ineq. (41) Se/2; then 
choose |6| sufficiently small so that the first term of the right member 
of inequality (41) S«/2. This establishes the continuity of S(f). 

We next verify (37), which concerns S,(f). Since h, isin £:(— ~, ~), 
H,(f) exists for all f € (— #, «). Thus, S,(f) as defined in (33) is 
meaningful, though perhaps infinite on a set of measure zero. With 


E(no f) as in (39), write 
S()-SN= E_ lH-(r+%) ‘~ E |a(r+ zy) 


= 5 [ert ge) - |e(r+g)/ | - soon. 42 


|n| Sno 











Now let « > 0 be arbitrary. From (39) we can choose no sufficiently 
large such that &(no, f) S ¢/2, for f € [0, 1/7]. With no so chosen, 


reed) hGeay He 
Now let #,(f) be the Fourier transform of 4,. Then 
[H(f)| = |H-(f) + A.(f)| S |H-(/)| + |4-(A/)I. 
Therefore 
| (f) |? 








Crees 


[n| Sno 


I 


Dl? + HANAN) + EO 
H(A)? + 2llhelaiella + Weel 


where || ||; denotes £: norm. Since ||/,||1 20, as r— ~, if 7 is 
sufficiently large, then 


WA IIA 


€ 1 
HDI — HOI 2 -$ oy, [E(-*, ). 4 


Inequalities (43) and (44) imply that, for all « > 0, there exists a 
ro(e) such that for all 7 = ro(e), 


S(N2ZSN-6 OSSSF (45) 
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Thus, for r 2 To(e), 


inf S.(f)2 inf S(f)-—-e=m-—e. (46) 
Osf S1/T OSf S1/T 


Letting 7— « and e— 0 in (46) yields (87). 


APPENDIX B 
Convergence of the Minimum Distance 


In this appendix, we shall verify (21), i.e., 
d(h,) ~d(h), asr—> ©. (47) 


From the definition of d(h,) (10), for arbitrary « > 0, we are assured 
of the existence of a w = u — v such that u, v € @y, and 


|hexw|| S d(h,) + «. (48) 
Repeating the steps in (34), we obtain 


1/T N-1 
[omol? = fo S.NIE wert |2af & Wine $.(] O lwnl. 49) 
0 k k=0 


From (37) we can choose 7 sufficiently large so that 


inf S,(f) 2m. (50) 


Hence, for such a choice of 7, 


< Ld(hs) + ef 


= |wz| S pean roe eer (51) 
Now 
d(h) S ||hew|| = |lhew — heew|| S |[heew|| + || ew] 
S d(h1) + ¢+ ||A*w]]. (52) 
Since 


(h-*w) (t) = LX h(t — kT)-wy,, 


we have, with 7 large enough to satisfy (51), 


[Fevwe S  [ovel Well = We bel s Mel cea + er. 3) 


Combining (52) and (58) yields for 7 sufficiently large (and e > 0 
arbitrary) 
lh, 


| ' 
2 Cath) + B, (54) 


d(h) S d(hiz) tet Well 
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Letting 7 and e — 0 yields 
d(h) S lim inf d(h,). (55) 


The identical argument with A and h, reversed yields for all + > 0, 
e> 0, 


a(h.) sah) + + Ml aay + 
so that (letting 7 ~~, e > 0) 
lim sup d(h,) S d(h). (56) 


Tao 


Inequalities (55) and (56) yield (47) or (21), completing the proof. 
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